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ABSTRACT 


MICROWAVE  RAr'l  OMETRY  AND  ITS 
APPLICATIONS  TO  MARINE 
METEOROLOGY  AND  OCEANOGRAPHY 
January  1969 
Jack  Frederick  Paris 
B.S.,  Terns  ASM  University 
B.S.,  University  of  Washington 
Directed  ’  •  Dr.  George  L.  ihif'-  Jr. 

Past  developments  in  microwave  radioraetry,  a  discipline  of 
remote  sensing,  are  reviewed  comprehensively  to  establish  a 
continuity  between  microwave  physics,  microwave  engineering,  and 
applications  to  marine  meteorology  and  oceanography.  Molecular 
oxygen,  water  vapor,  and  liquid  water  absorb,  emit,  and  scatter 
microwave  radiation.  The  ocean  reflects  or  scatters  a  majority 
of  the  microwave  radiation  incident  upon  its  surface  and  emits 
microwave  energy.  The  functional  relationships  between  these 
interactions  and  the  physical  state  of  the  atmosphere  and  hydro¬ 
sphere  are  described  in  detail. 

The  emission  of  sea  water  is  almost  constant  with  temperature 

and  salinity  for  frequencies  of  8  GHz  to  30  GHz,  is  linearly 

proportional  to  temperature  for  frequencies  near  5.4  GHz,  and 

is  strongly  dependent  upon  salinity  for  frequencies  near  1  GHz. 

Sea  foam  causes  the  microwave  emission  of  the  sea  to  increase 

ii  i 


greatly  for  all  microwave  frequencies  and  angles  of  viewing. 

The  natural  variability  of  water  vapor  in  the  atmosphere 
affects  the  absolute  value  of  upwelling  microwave  radiation 
greatly  for  frequencies  near  22.235  GHz.  The  attenuation  and 
emission  of  the  atmosphere  are  predictable  for  frequencies  less 
than  8  GHz. 

Microwave  measurements  near  5.4  GHz  and  1  GHz  may  be  used 
to  survey,  remotely,  the  temperature  and  salinity  of  coastal 
water.  Microwave  measurements  near  22.235  GHz  may  be  used  to 
measure  the  precipitable  water  in  the  atmosphere. 

It  is  recommended  that  (1)  electric  properties  of  sea 
water  be  measured  precisely,  (2)  controlled  environmental  tests 
of  current  microwave  theory  be  conducted,  (3)  the  problem  of  the 
microwave  emission  and  extinction  of  raining  clouds  over  water 
be  studied  in  great  detail  using  realistic  models  of  clouds  and 
using  the  Mie  theory,  and  (4)  the  emissive  properties  of  sea 
foam,  bubbler.,  and  rough  sea  surfaces  be  determined  through 
empiricism. 
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CHAPTER  I 


INTRODUCTION 

Nature  of  the  Problem 

Meteorological  and  oceanographic  studies  are  based  upon  ob¬ 
servations  of  the  environment.  Although  much  progress  has  been 
made  during  this  century  in  the  art  and  science  of  environmental 
observations,  much  of  the  hydrosphere  and  atmosphere  goes  unsur¬ 
veyed  or  is  inadequately  surveyed.  An  artificial  satellite  in 
Earth  orbit  is  a  platform  from  which  observations  of  phenomena 
on  Earth  can  be  made  with  unprecedented  coverage,  synopticity, 
and  repetitiveness.  These  observations  are  necessarily  made  by 
remote  sensors. 

Although  some  remote  sensors  measure  the  gravitational,  elec¬ 
trical,  or  magnetic  force  field,  the  vast  majority  of  remote  sen¬ 
sors  are  radiation  sensors  which  are  used  to  determine  the  dis¬ 
tribution  of  the  intensity  of  the  radiational  field  in  the  vicin¬ 
ity  of  the  satellite.  Since  this  distribution  is  a  function  of 
interactions  that  have  taken  place  between  radiation  and  matter 
along  the  path  taken  by  the  radiant  energy,  scientists  are  able 
to  relate  the  remotely  sensed  data  to  physical  phenomena  in  the 
field  of  view  of  the  sensor. 

The  citations  on  the  following  pages  follow  the  style  of  the 
Journal  of  Geophysical  Research. 


t 


Remote  sensors  may  be  classified  as  passive,  active,  qualita¬ 
tive,  or  quantitative.  Most  remote  sensors  are  passive  sensors, 
that  is,  they  detect  natural  radiation.  Sensors  which  provide 
their  cwn  illumination,  such  as  radar  sensors,  are  called  active 
sensors.  The  difference  between  qualitative  and  quantitative 
remote  sensors  is  not  distinct.  In  general,  if  the  output  of 
a  sensor  is  recorded  photographically,  it  is  a  qualitative  sensor. 
If  the  output  of  a  sensor  lends  itself  to  calibration  and  numer¬ 
ation,  it  is  a  quantitative  sensor. 

Gases,  water  droplets,  and  other  atmospheric  constituents 
selectively  absorb,  scatter,  and  emit  radiation.  Scientists  who 
desire  to  view  surface  phenomena  from  satellites  select  frequency 
domains  in  the  electromagnetic  spectrum  (Figure  1)  of  maximum 
atmospheric  transparency.  These  domains  are  called  atmospheric 
'windows'  and  are  listed  in  Table  I.  These  'windows'  are  neither 
completely  transparent  nor  sharply  bound.  The  optical  and  in¬ 
frared  windows  become  opaque  when  clouds  are  in  the  field  of  view. 
The  upper  end  of  the  radio  window  may  becc'.v-  opaque  over  raining 
clouds.  Some  scientists  may  desire  to  itudy  certain  atmospheric 
phenomena  using  non-window  regions. 

In  1931,  thermal  radiation  was  first  detected  at  radio  fre¬ 
quencies  from  extra-terrestrial  sources  (for  example,  the  Milky 
Way  Galaxy)  using  a  crude  radio  receiver  operating  at  20.5  MHz^ 

^■MHz  =  10^  cycles  per  second  (one  mega-Hertz) . 


WAVELENGH 


*  3m  30  m  300  m  -3cm  3cm  30cm  3m  30m  300m  3km  30km  300km 

_ t _ I _ t _ I _ < _ I _ I _ I _ » _ I _ » _ 

“w.C.Vi,  U"'  ,A0I°  U  *M'°  « 


ATMOSPHERIC  TRANSMISSION 


VIIRATION  ROTATIONS! ELECTRIC,  MAGNETIC 

FIELD  FLUCTUATIONS 


PHENOMENA 

DETECTED 


*  Tftttttff!  tttftltttt 
•020  >0i8  10‘6  I014  1012  I010  I08  I06  to4  10*  I 

- -  FREQUENCY  (cps) 


NASA  SM65- 15182 
T/8/65 


Figure  1.  The  electromagnetic  spectrum 
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and  30  MHz.  Highly  sensitive  radio  receiveis  and  highly  direc¬ 
tional  radio  antennas  were  being  developed.  This  research, 
coupled  with  intensive  studies  of  the  absorption  and  scattering 
properties  of  atmospheric  gases  and  water  droplets  for  radar 
propagation,  led  to  the  development  of  radio  astronomy.  Micro- 
wave  radiometry  evolved  from  the  science  of  radio  astronomy. 

A  microwave  radiometer  is  a  quantitative,  passive  remote  sen¬ 
sor  i>at  makes  precise  measurements  of  the  intensity  of  thermal 
radiation  at  microwave  frequencies.  Early  receivers  in  microwave 
radiometry  suffered  mostly  from  unsteady  amplification.  Dicke 
(1946 )  proposed  a  technique  to  lessen  errors  in  measurements  of 
microwave  radiance  due  to  amplification  fluctuations.  By  switching 
rapidly  from  the  power  received  by  the  antenna  system  to  the  known 
power  of  an  internal  source,  he  was  able  to  reduce  the  effect  of 
amplification  drift. 

Much  of  the  scientific  community  in  remote  sensing  is  pres¬ 
ently  interested  in  possible  meteorological  and  oceanographic 
applications  of  radiometric  measurements  at  microwave  frequencies 
from  spacecraft  and  aircraft. 

The  pvirpose  of  this  study  is  to  review  comorehensivelv  past 
developments  related  to  mi  rowave  radiometry  and  to  present  them 
in  a  manner  understandable  to  oceanographers  and  meteorologists 


who  do, ice  to  use  this  remote  sensor  as  a  tool  in  their  research. 


Past  Studies 
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During  World  War  II,  many  investigations  were  made  into  prob¬ 
lems  of  propagation  of  radar  waves  in  the  atmosphere.  After  the 
war,  the  results  of  these  investigations  appeared  in  scientific 
publications  for  the  first  time.  Of  concern  to  microwave  radio— 
metry  are  papers  dealing  with  the  absorption  and  scattering  pro¬ 
cesses  in  the  atmosphere,  the  emission  o.  the  sea  surface,  the 
electrical  properties  of  water,  and  developments  in  the  engineer¬ 
ing  of  microwave  radiometers. 

Mie  (1908)  formulated  the  complete  theory  of  scattering  and 
absorption  of  electromagnetic  waves  by  a  homogenous  dielectric 
sphere.  Stratton  (1930)  computed  values  of  attenuation  in  uni¬ 
form  clouds  of  water  droplets  using  the  Mie  theory;  however,  he 
erroneously  assumed  that  the  imaginary  part  of  the  complex  index 
of  refraction  (and  hence  the  absorption)  was  zero.  Franz  (1940) 
correctly  computed  the  attenuation  of  microwaves  in  fog;  but,  he 
did  not  consider  rain  droplets  or  ice  crystals.  Using  the  exper¬ 
imentally  determined  values  of  the  complex  dielectric  constant  of 
water  by  Saxton  and  Lane  (1946)  and  the  complex  dielectric  constant 
of  ice  by  Dunsmuir  and  Lamb  (1945),  Ryde  (1946)  computed  the  ab¬ 
sorption  and  scattering  of  microwaves  by  fog,  cloud,  rain,  hail, 
and  srow  using  the  Mie  theory  as  developed  by  Stratton  (1941). 

He  assumed  the  distribution  of  drop  sizes  given  by  Laws  and  Parsons 


(1943).  Goldstein  (1951)  elaborated  on  Ryde's  work.  Aden  and 
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Kerker  (1951)  extended  the  Mie  theory  to  concentric  spheres  of 
homogeneous  dielectric  materials.  Gunn  and  East  (1954)  presented 
a  very  complete  treatment  of  osorption  of  microwaves  by  fogs, 
rain,  atmospheric  gases,  clouds,  hail,  and  ice.  This  work  sum¬ 
marized  all  previous  investigations  and  updated  Ryde's  calcula¬ 
tions.  Imai  (1957)  determined  the  absorption  and  scattering  co¬ 
efficients  of  rain  drops  to  a  great  degree  of  accuracy.  Using 
Aden's  method  for  Mie  theory,  Stephens  (1961)  computed  the 
attenuation  coefficient  of  water  droplets. 

Buettner  (1963,  1965)  suggested  that  the  microwave  emission  of 
a  cloud  could  be  related  to  its  distribution  of  liquid  water. 
Cummings  and  Hull  (1966)  observed  large  increases  in  sky  radiance 
under  various  cloud  and  rain  conditions  and  confirmed  the  pre¬ 
dictions  of  Buettner.  An  increase  in  the  intensity  of  upwelling 
radiation  was  observed  over  marine  clouds  during  flight  tests 
of  a  scanning  microwave  radiometer  operating  at  19.35  GHz  over  the 
Gulf  of  Mexico  (Catoe  et  al. ,  1967).  Kreiss  (1968)  computed  the¬ 
oretical  values  of  the  intensity  of  upwelling  microwave  radiation 
over  clouds  of  various  thicknesses,  heights, and  concentrations  of 
liquid  water.  Using  microwave  radiometers  from  the  ground.  Decker 
and  Dutton  (1968)  determined  the  distribution  of  liquid  water  in 
thunderstorms  over  Colorado. 

Van  Vleck  (1947a,  1947b,  1951a),  and  Van  Vleck  and  Weisskopf 
(1945)  determined  the  absorption  of  microwaves  by  oolecular 
oxygen  and  water  vapor  from  quantum  theory.  This  outstanding 


work  stands  essentially  unaltered  to  the  present  day  except  for 
a  correction  for  microwave  absorption  by  pressure-broadened  ab¬ 
sorption-lines  in  the  infrared  spectrum  (Figure  1)  of  water  vapor 
(Becker  and  Autler ,  1946) . 

The  absorption  by  molecular  oxygen  for  wavelengths  near  5 
millimeters  has  received  much  attention  lately.  Meeks  and  Lilley 
(1963)  proposed  that  the  microwave  emission  from  molecular  oxygen 
in  the  vicinity  of  5  mm  could  be  used  to  determine  the  distribu¬ 
tion  of  air  temperature  with  height.  Westwater  (1965)  determined 
profiles  of  air  temperature  by  inverting  measurements  of  the  micro- 
wave  radiance  of  the  sky  from  ground-based  microwave  radiometers. 

A  discussion  of  present  efforts  being  made  to  invert  microwave 
measurements  of  oxygen  emission  in  the  atmosphere  to  obtain  temper¬ 
ature  profiles  in  the  atmosphere  will  not  be  included  in  this 
paper.  One  may  determine  the  present  state  of  these  efforts  by 
referring  to  Gille  (1968). 

The  absorption  by  water  vapor  in  the  atmosphere  is  much  greater 
than  that  by  molecular  oxygen  near  wavelengths  of  1.35  centimeters 
(Van  Vleck,  1947b).  Thus,  it  should  be  possible  to  obtain  infor¬ 
mation  about  the  distribution  of  water  vapor  in  the  atmosphere 
from  microwave  measurements.  Barrett  and  Chung  (1962),  and  Staelin 
(1966)  have  studied  these  possibilities. 

The  microwave  emission  and  reflection  of  the  sea  surface  has 
been  studied  by  many  investigators.  The  reflection  and  emission  of 
homogenous,  flat  dielectrics  is  well  known  (Jordan,  1950). 


Peake  (1959)  defined  the  problem  of  scattering  and  emission  from 


homogenous  dielectrics  having  an  irregular  surface  and  used  various 
models  of  roughness  to  compute  the  emission  of  various  land  sur¬ 
faces  (Chen  and  Peake,  1961).  The  theory  for  a  smooth  multi-layer 
dielectric  was  used  by  Pascalar  and  Sakamoto  (1965)  to  compute  the 
emission  and  reflection  of  ice  sheets  over  water.  Barath  (1965) 
and  Mardon  (1965)  suggested  that  measurements  of  ice  distribution, 
ice  thickness,  water  temperature,  and  sea  state  could  be  made  from 
aloft  using  microwave  radiometers.  McAlister  and  McLeish  (1965) 
proposed  a  method  whereby  the  total  vertical  heat  flow  from  the 
ocean's  surface  could  be  determined  by  dual-frequency  microwave 
radiometers.  Conway  and  Sakamoto  (1965)  reported  on  measurements 
of  microwave  radiance  of  various  terrestrial  objects  at  13,  18, 

35,  and  70  GHz?,  Stogryn  (1967)  used  the  Kirchhof f  Approximation 
and  the  distribution  of  wave  slopes  obtained  by  Cox  and  Munk 
(1954a,  1954b)  to  determine  the  effects  of  sea  state  on  microwave 
radiance  of  the  sea  at  19.4  GHz  and  35  GHz. 

In  1967,  the  U.S.  Coast  Guard  began  to  use  microwave  scanners 
over  the  North  Atlantic  Ocean  to  detect  icebergs  through  fog 
(Roeder ,  1967).  Sirounian  (1968)  and  Paris  (1968)  determined  the 
effect  of  salinity  on  microwave  radiance  at  3  GHz.  Williams  (1968) 
showed  that  microwave  emission  of  foam  greatly  exceeds  that  of  a 

2  9 

GHz  =  10  cycles  per  second  (one  giga-Hertz) . 
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smooth  sea  surface. 

Recently,  measurements  of  microwave  radiance  over  the  Gulf 
of  Mexico  over  large  ranges  of  temperatures,  salinities,  and 
concentrations  of  liquid  water  have  been  furnished  to  the  author 
by  Blinn  (1967a),  Nordberg  (personal  communication,  1968),  and 
the  Manned  Spacecraft  Center.  These  data  were  taken  by  tne  MR-62 
and  MR-64  microwave  radiometers  (Blinn,  1967a)  at  9.3,  15.8, 

22.235,  and  34.0  GHz  and  by  a  scanning  19.35  GHz  microwave  radio¬ 
meter  built  by  Aerojet  General  Corportation  (Catoe  et  al.,  1967; 
Tobin,  1967). 

The  only  microwave  radiometer  that  has  been  flown  in  a  space¬ 
craft  is  one  dual-frequency  microwave  radiometer  (model  MR-62) 
carried  by  the  Mariner  II  spacecraft  on  the  Venus  mission  (Barath 
et  al. ,  1964;  Glass tone,  1965). 

Although  many  aspects  of  applied  microwave  radiometry  have 
been  adequately  considered  by  the  above  investigators,  no  presen¬ 
tation  of  microwave  radiometry  exists  in  the  literature  which  com¬ 
prehensively  considers  the  total  effects  of  water  temperature, 
water  salinity,  surface  roughness,  foam,  absorption  and  emission  by 
atmospheric  gases,  fogs,  clouds,  and  rain  on  microwave  radiometric 
measurements . 

In  addition  to  the  physical  theory  of  microwave  emission  and 
extinction  in  the  ocean  and  atmosphere,  one  must  be  aware  of  the 
capabilities  and  limitations  of  the  radiometric  instruments. 

Dicke  (1946),  Bracewell  (1962),  Ko  (1964),  and  Barath  (1967)  fully 
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discuss  these  factors.  Also,  since  microwave  radiometers  are  very 
sensitive  instruments  and  since  man-made  sources  of  microwave  ra¬ 
diation  may  be  many  times  more  intense  than  natural  sources  of 
microwave  radiation,  scientists  must  use  frequencies  void  of  man¬ 
made  interferences.  In  1964,  certain  frequencies  were  set  aside 
for  use  by  radio  astronomers.  These  radio  astronomy  bands  are 
shown  in  Table  II  (Smith-Rose,  1964). 


Table  II.  Microwave  bands  reserved 
Smith-Rose,  1964). 

Frequency  band  (GHz) _ 

1.4  -  1.427  (hydrogen  line) 

1.66  -  1.69  (OH  line) 

2.69  -  2.70 
3.165  -  3.195 

4.8  -  4.81 
4.99  -  5.00 

5.8  -  5.815 

8.68  -  8.70 

10.68  -  10.70 
15.35  -  15.4 

19.3  -  19.4 

31.3  -  31.5 
33  -  33.4 

33.4  -  34.0 

36.5  -  37.5 


for  radio  astronomy  (after 


A1 locations 


Exclusive  world-wide 
secondarv  (.urinary: 
meteorological  satellites 
exclusive  world-wide). 


Europe,  Cuba  and  India 


Objective  and  Scope 

The  objective  of  this  study  is  to  review  comnrehensivelv 
past  and  present  activities  in  microwave  radiometry  as  they  are 
related  to  marine  applications.  The  development  of  certain 
ohysical  and  engineering  concepts  will  be  presented  in  a  rigorous 
manner  so  that  the  nature  of  the  simplifying  assumptions  can  be 


fully  understood.  This  theory  will  be  compared  to  measurements 
when  appropriate  and  available. 

The  purposes  of  this  review  are  to  lay  a  firm  foundation  for 
future  research  in  specific  areas  of  concern  to  microwave  radio- 
metry  and  hopefully,  to  engender  interest  among  other  environmental 
scientists  who  may  be  able  to  use  microwave  radiometers  as  a  tool 
in  their  research. 

Since  the  problem  of  determining  the  vertical  distribution  of 
air  temperature  by  inverting  microwave  measurements  of  oxygen 
emission  is  separate  from  all  other  aspects  of  applied  microwave 
radiometry,  it  will  not  be  included  in  the  scope  of  this  paper. 

The  main  interests  of  this  study,  as  far  as  applications  are 
concerned,  are  techniques  that  are  beine  developed  for  marine 
environmental  measurements. 

In  Chapter  II  the  fundamental  concepts  of  the  electrical  and 
magnetic  properties  of  water,  ice,  and  air  are  defined;  practical 
expressions  are  developed  for  the  complex  dielectric  constant  of 
water  as  a  function  of  its  temperature,  its  salinity,  and  the 
frequency  of  observation. 

The  fundamental  problem  of  radiative  transfer  is  discussed  in 
Chapter  III,  and  certain  aspects  of  radiation  in  general  are  dis¬ 
cussed.  The  concept  of  brightness  temperature  is  introduced. 

In  Chapter  IV,  the  extinction  and  emission  processes  in  the 
atmosphere  are  discussed,  and  the  equation  of  transfer  is  solved 
for  various  meteorological  conditions. 
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The  emission  and  scattering  by  the  hydrosphere  is  determined 
in  Chapter  V  for  various  conditions. 

The  results  of  Chapter  IV  and  V  are  combined  in  Chapter  VI 
to  give  expressions  for  the  intensity  of  the  upwelling  radiation  as 
a  function  of  altitude,  frequency,  polarization,  angle,  and  physi¬ 
cal  state. 

Internal  and  external  aspects  of  microwave  radiometers  are 
discussed  in  Chapter  VII,  and  the  effect  of  instrument  character¬ 
istics  on  the  measured  flux  of  upwelling  microwave  radiation  is 
determined. 

Available  measurements  of  microwave  radiance  are  compared 
to  theoretical  values  in  Chapter  VIII,  and  conclusions  and  recom¬ 
mendations  are  made  in  Chapter  IX. 
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CHAPTER  II 

ELECTRIC  AND  MAGNETIC  PROPERTIES  OF  MATTER 
Maxwell's  Equations 

Maxwell  (1892)  generalized  the  laws  of  Ampere,  Faraday, 
and  Gauss  into  a  set  of  electromagnetic  field  equations  which 
can  be  used  to  describe  the  macroscopic  properties  of  matter. 
Maxwell's  equations  for  a  chargeless  medium  may  be  expressed  in 
the  Giorgl  (1901)  systems  of  units  as  follows: 


V  • 

-* 

E  = 

0 

(volts/m^t . 

(1) 

v  • 

-*■ 

H  = 

0 

(amps/m^) . 

(2) 

V  x 

-* 

E  = 

3H 

~  U°  3t 

(volts/m^) . 

(3) 

V  x 

-► 

H  = 

°l+  -’••If 

(amps/m~) . 

(4) 

E  is  the  electric  field  vector  (volts  per  meter),  H  is  the 
magnetic  field  vector  (amperes  per  meter),  vQ  is  the  permeability 
(Henrys  per  meter)  of  a  vacuum,  t  is  time  (seconds),  o  is  the 
conductivity  (mhos  per  meter),  e'  is  the  real  relative  permittivity 
(dimensionless),  and  zq  is  the  permittivity  (farads  per  meter)  of 
a  vacuum.  By  international  consent  u0  is  4 n  x  10“^  henrys  per 
meter.  It  is  assumed  that  the  permeability  of  all  mediums  dis¬ 
cussed  in  this  paper  is  yQ.  eQ  is  approximately  10~9/3f>v  farads 

^The  Giorgi  system  or  rationalized  MKS  system  of  units  uses 
the  basic  units:  meter,  kilogram,  second,  and  coulomb. 
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per  meter  (Von  Hlppel,  1954) . 

It  is  assumed  that  all  components  of  £  and  1?  are  sinusoidal 
in  nature,  that  is, 

E  =  Em  exp  ( jwt)  (volts/m) ,  (5) 

and 

H  =  exp  (jut)  (amps/m) ,  (6) 

where  h  and  E  are  the  maximum  values  of  H  and  E,  respectively, 
mm 

w  is  the  angular  frequency  (radians  per  second),  and'j  =  /-i . 

According  to  Von  Hippel,  o  represents  an  effective  conductivity 
which  includes  all  dissipative  effects.  The  customary  approach 
has  been  to  define  a  relative  loss  factor,  e",  such  that 


o  «  weoe"  (mhos/ra).  (7) 

Substitution  of  (5)  and  (7)  into  (4)  gives: 


V  x  H  =  (we  e"  +  jwe'e  )E 
o  J  o 

=  jweo(e'  -  je")E 

-V 

=  (e'  -  je")|^  (amps/m) .  (8) 


Complex  Dielectric  Constant 


Theory 


The  parenthetical  expression  in  (8)  is  dimensionless  and  is 
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called  the  complex  dielectric  constant  of  the  medium.  According 
to  Sucher  and  Fox  (1963),  e*  is  associated  with  the  ability  of 
a  material  to  store  electric  energy,  and  e"  is  associated  with 
the  losses  that  occur  in  the  material.  These  losses  include 
viscous  dissipation  of  energy  in  polar  liquids  as  well  as  ohmic 
losses.  The  curl  of  (3)  is: 

VX(7x?)=-p  V  x 

o  at 


»  -  p  3(v  x  5)/3t  (volts/m3).  (9) 

o 

By  the  use  of  a  vector  identity  and  (8),  (9)  becomes: 

2 

V  (v  -f)  -  7  22  -  -  uoeo(e'  -  je">  (volts/m3).  (10) 

Since  7  •  2  is  zero  by  (1), 


^  .  p  E  (£'  -  je”) 

o  o  3t 

Also,  the  curl  of  (8)  is: 

7  x  ( v  x  3)  -  e  (e'  -  je")3( V  x  f)/a t 
0 

By  the  use  of  a  vector  identity  and  (3) , 

2 

7(7  -  5)  -  72S  -  -y  e  (e’  -  je") 

o  0  dt 

Since  7  *3  is  zero  by  (2) , 

72S  -  p  e  (e1  -  je")  -Mf 
00  0 1 


(volts/m^) . 

(amps/m  ) . 
(12)  becomes: 

2 

(amps/m  ) . 

(araps/m3) . 


(ID 


(12) 


(13) 


(14) 


i/ 


Equations  (11)  and  (1A)  are  the  wave  equations  of  the  electro¬ 
magnetic  field  in  a  chargeless  medium  having  permeability  equal 

to  u  . 

c 

Without  loss  of  generality  one  may  assume  that  the  electro¬ 
magnetic  wave  described  by  (11)  and  (1A)  is  traveling  along  an 
x-axis  in  space  and  is  a  plane  wave.  Thus, 


3E  =  3E  _  3H  m  3H 
3y  3z  3y  3z 


(15) 


and  (11)  and  (1A)  reduce  to: 


voeo(e'  "  St"*  Jjl 


(volts/m^),  (16) 


and 


(amps/n)3) ,  (17) 


respectively. 

A  set  of  solutions  to  (16)  and  (17)  according  to  Ramo,  et  al., 
(1965)  is: 

E  =  Eq  exp(-ax)  exp(-jgx)  (volts/m),  (18) 

and 

H  =  H  exp(-ax)  exp(-jgx)  (amns/m).  (19) 

o 

E  and  H  are  initial  values  ol:  E  and  H,  respectively,  where: 
oo 
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a  =  a)  A\iQzQz'  H)  /I  +  (z"/z')?-  1  (nepers)-4,  (20) 

dad 

8  *  a)  /<p0E0e72)  /I  +  (c,7e,)*+l  (nepers).  (21) 

a  is  the  attenuation  factor  (nepers) ,  and  g is  the  phase  fac¬ 
tor  (nepers) . 

If  e"  is  zero  (and  hence  the  medium  is  lossless),  a  is  zero, 

and 

6  =  a)/  \x  z  e'  (nepers).  (22) 

o  o 

Since  the  form  of  (16)  and  (17)  is  that  of  a  wave  traveling 
at  velocity  v  (meter) 

V  =  1 N  M  z  (e'  -  je”)  (m/sec),  (23) 

o  o 

and  since  the  index  of  refraction,  m,  is  defined  as: 

m  =  c/v  (dimensionless),  (24) 

then: 

m^  =  c'  -  je"  (dimensionless).  (25) 

c  is  the  speed  of  propagation  of  electromagnetic  energy  in 
a  vacuum  and  is  given  by:  (26) 

c  =  1/  /  u  t  (ra/'sec).  (26) 

o  o 

^  A  neper  is  an  engineering  unit  which  means  per  meter . 
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m  is  complex  and  may  be  defined  as: 

m  =  n  +  jK  (dimensionless).  (27) 

n  is  the  real  index  of  refraction,  ai  1  k  is  the  index  of 
absorption. 

From  (25)  and  (27)  we  derive  the  relations: 

e'  =  n^  -  (dimensionless),  (28) 

and 

e"  =  2n<  (dimensionless).  (29) 

Since  it  has  been  assumed  that  the  permeability  of  all  matter 
under  study  is  equal  to  uq,  it  is  only  necessary  to  know  c' and 
as  a  function  of  frequency,  temperature, and  physical  state  to 
completely  describe  the  electric  and  magnetic  properties  of  the 
matter  in  question. 

Debye  (1929)  suggested  that  the  complex  dielectric  constant 
of  a  polar  liquid,  such  as  water,  could  be  expressed  as: 

‘  +  <es  -  O/d  + 

=  [ e  +(e_-c)/(l+  w^t^) 

1  00  5  CD 

-  j  [lcs  ”  co#)wi/(1  +  u)2-2)]  (dimensionless).  (30) 

e  is  the  static  relative  permittivity,  t  is  the  relaxation  time 
s 

(seconds)  and  e  is  the  relative  permittivity  at  high  frequencies. 
In  reality,  e  varies  with  frequency  and  physical  state;  however, 

O) 

since  (3)  is  to  be  used  only  for  radio  and  microwave  frequencies, 
a  constant  value  is  assumed  for  c  •  In  general,  z  and  t  are 

oo  g 

functions  of  the  physical  state  but  not  of  the  frequency. 
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Lane  and  Saxton  (1951,  1952;  Saxton  and  Lane,  1946;  Saxton 
and  Lane,  1952),  Saxton  (1946a,  1946b,  1946c,  1951)  and  Collie, 
et  al . ,  (1947)  have  studied  the  electric  properties  of  liquid 
water  with  and  without  sodium  choloride  in  solution.  Measurements 
have  been  made  of  e’  and  e"  by  these  investigators  for  frequencies 
of  3,  9.4,  10,  18,  24,  24.2, and  48.5  GHz,  for  temperatures  ranging 
from  -8°  C  to  50°  C,  and  for  concentrations  of  sodium  chloride 
ranging  up  to  3N.  These  measurements  indicate  that  (30)  is  valid 
for  wafer  having  concentrations  of  sodium  chloride  in  solution 
(Lane  and  Saxton,  1952) . 

The  relative  loss  factor  is  quite  appreciable  in  pure  water 
even  though  the  ionic  conductivity  of  pure  water  is  quite  small- 
This  fact  indicates  the  importance  of  losses  due  to  viscous  dis¬ 
sipation  known  as  Debye  relaxation.  The  addition  of  sodium 
chloride  slightly  affects  losses  due  to  viscous  dissipation  and, 
of  course,  introduces  additional  losses  due  to  the  ionic  conduc¬ 
tivity  as  the  result  of  the  sodium  and  chloride  ions  present.  It 
has  been  customary  to  separate  these  two  mechanisms  of  loss  by 
introducing  an  additional  factor  into  the  expression  for  s". 

From  (30),  it  follows  that: 

e*  “  e  +  (e  -  e  )/(l  +  ui^t^)  (dimensionless),  (31) 

oo  S  00 

and 


e"  =  wT(e  -e  )/(l  +  u^t^) 
s  00 


(dimensionless).  (32) 


*>1 


Now  it  is  assumed  that  (31)  and  (32)  are  valid  for  pure 
water  where  eg  and  i  are  functions  of  temperature,  T,  (degrees 
Kelvin).  In  the  case  of  salt  water,  it  can  be  assumed  (Lane  and 
Saxton.  1952)  that: 

e*  *  e  +  (e  -  e  ) / (3  +  ;j2r2  )  (dimensionless),  (33) 

00  o  co 

and 

e"  =  wt(-  -  e  )/(l  +  w2t2)  +  o./jf  (dimensionless)  (34) 
s  10 

where  is  the  ionic  conductivity  of  the  solution  (mhos  per  meter) . 

In  Lane's  and  Saxton's  pacers,  electrostatic  units  are  used 
for  e  and  o^.  Also,  e'  and  e"  are  the  real  and  imaginary  parts 
of  the  complex  permittivity  which  are  numerically  equal  to  the 
real  and  imaginary  parts  of  the  dielectric  constant  used  in  this 
study. 

In  (33)  and  (34),  e  ,  T,and  o.  are  functions  of  temperature 

s  i 

and  salinity  only. 

There  now  arises  a  problem  concerning  the  salt  concentration 
of  a  solution  of  pure  water  and  sodium  chloride  and  the  salt  con¬ 
centration  of  ocean  water.  All  previous  measurements  of  the  effects 
of  salts  in  solution  on  the  dielectric  properties  of  water  have 
been  made  using  only  sodium  chloride.  Ocean  water,  on  the  other 
hand,  has  many  types  of  salts  in  solution.  Table  III  ^ists  the 
salt  constituents  most  abundantly  found  in  sea  water  and  gives 
their  percent  contribution  by  weight  (McLellan,  1965) . 
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Table  III. 


Principal  io.  Lc  constituents  of  sea  water  (34.4  o/oo 
salinity)  (after  McLella_i,  19b5). 


Amount  (g) 


Percentage 


per  kg  of  Molecular  or  of  total 

Species _ solution _  atomic  weight*  salt 

Cations: 


Sodium 

10.47 

22.997 

30.4 

Magnesium 

1.28 

24.32 

3.7 

Calcium 

0.41 

40.08 

1.2 

Potassium 

0.38 

39.096 

1.1 

Strontium 

0.013 

87.63 

0.05 

Anions : 

Chloride 

18.97 

25.457 

55.2 

Sulfate 

2.65 

32.06** 

7.7 

Bromide 

0.065 

79.916 

0.2 

Bicarbonate 

0.14 

0.4 

Borate 

0.027 

10.82 

0.08 

*1>.  Van  Nostrand  Company,  Inc.  (1958). 
**Value  for  sulfur. 


From  Table  II,  one  sees  that  85.6  percent  of  the  total  dis¬ 
solved  salts  are  sodium  and  chloride.  The  mean  molecular  weight  of 
sodium  chloride  in  sea  water  is  31.0.  The  mean  molecular  weight 
of  all  cations  and  anions  listed  in  Table  III  is  31.2.  Based  on 
this  fact,  it  is  assumed  that  a  solution  of  sodium  chloride  and 
water  approximates  sea  water  as  far  as  its  electrical  properties 
are  concerned. 

According  to  McLellan  (1965)  the  chlorinitv  of  a  solution  is 
the  total  amount  of  chlorine,  bromine,  and  iodine  in  grams  con¬ 
tained  in  one  kilogram  of  sea  water  assuming  the  bromine  and  iodine 
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to  be  replaced  by  chlorine.  Assuming  that  the  relative  proportions 
of  the  major  constituents  of  sea  water  remain  constant,  one  may 

relate  the  salinity  of  sea  water  to  its  chlorinity  by  the  following 
expression  (Cochrane,  personal  communication,  1967): 

salinity  =  1.80655  x  chlorinity  (o/oo).  (35) 

Practical  treatment 

Lane  and  Saxton  (1952)  have  determined  values  of  e  ,  e  , 

-  ®  s 

and  t  at  specific  values  of  temperatures  and  concentrations. 

They  expressed  concentration  in  terms  of  normal  solutions.  A 
solution  that  has  a  normality  of  unity  contains  one  gram-equiva¬ 
lent-molecular-weight  of  dissolved  substance  per  liter  of  solu¬ 
tion  (Chemical  Rubber  Co.,  1959).  Since  the  valences  of  sodium 
and  chloride  are  unity,  one  normal  solution  contains  58.454  grams 
of  sodium  chloride  per  liter.  The  mass  of  one  liter  of  solution 
is  a  function  of  its  temperature,  the  pressure  and  the  percentage 
of  solute  by  weight. 

The  salinity  of  a  solution  may  be  defined  as  the  total  amount 
of  solid  material  in  grams  contained  in  one  kilogram  of  solution 
when  all  carbonate  has  been  converted  to  oxide,  the  bromine  and 
iodine  are  replaced  by  chlorine,  and  all  organic  matter  oxidized 
(McLellan,  1965).  This  definition  was  formulated  for  sea  water; 
it  applies  to  sodium  chloride  and  water,  fhus,  one  may  express 
the  percentage  of  solute  by  weight  as  a  salinity. 


24 


Table  IV  gives  the  specific  gravity  of  solutions  of  sodium 
chloride  and  water  as  a  function  of  temperature  and  salinity  at 
one  atmosphere  of  pressure  (Chemical  Rubber  Co. ,  1959). 

Table  IV.  Specific  gravity  of  solution  of  sodium  chloride  and 
water  as  a  function  of  temperature  and  salinity  at 
one  atmosphere  of  pressure. 


Salinity, 

Temperature,  T(°C) 

S  (o/oo) 

-8 

0 

4  10  20 

30 

0 

0.99869 

0.99987 

1.0000  0.99973  0.99823 

0.99567 

10 

— 

-- 

1.0053 

— 

20 

30 

40 

— 

— 

1.0125 

— 

_ 

_ 

1.0268 

_ 

Since,  in  :.ome  cases,  the  specific  gravity  in  Table  TV  is 
below  unity,  and,  in  some  cases,  it  is  above  unity,  it  is  assumed 
that  the  specific  gravity  of  a  solution  of  sodium  chloride  and 
water  is  approximately  unity.  Thus  there  are  approximately  58.45 
grams  of  sodium  chloride  per  kilogram  of  solution  for  a  one- 
normal  solution,  that  is: 

S  =  58.45N  (o/oo),  (35) 

where  N  is  the  concentration  of  the  solution  expressed  in  terms 
of  normality. 

Table  V  gives  values  of  r  ,  t,  and  as  a  function  of 
temperature  and  salinity  from  measurement  made  by  Saxton  and 


Lane  (1952). 
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Using  the  values  of  es,t,and  given  in  Table  V,  Huebner  and 
McGraw  of  Texas  A&M  University  formed  a  set  of  second  order  re¬ 
gression  equations  in  temperature  (T,  °C)  and  salinity  (S,  o/oo) 
to  express  es>and  t.  These  expressions  are  as  follows: 

Oi(T,S)  =  (O.Q8533905)S  -  (0.000135394)S2  +  (0.004047859111)T 

+  (0. 0000277351188) T2  +  (0.03255122)ST 

-  (0. 000007332967 22) S2T  +  (0.0000000549132477)S2T2 

-  (0 . 00000883166144) ST2  (mhos/m) .  (36) 

es(T,S)  =  87.336049  -  (0.34438426)S  +  (0.00063045524)S2 

-  (0.39006802)T  +  (0.00068027202)T2  +  (0.0023219662)ST 

-  (0.0000098371722)S2T  +  (0. )00000099566511)S2T2 

-  (0 . 000012220872) ST2  (dimensionless).  (37) 

t(T,S)  «  18.494012  -  (0. 04522073) S  +  (0. 000082162302) S2 

-  (0.54708843)T  -r  (0.0058367346)T2  +  (0.0017348405)ST 

-  (0.0000027579928) S2T2  +  (0.000173711)ST2  (sec). (38) 

e  is  assumed  to  be  4.9  (Lane  and  Saxton,  1952). 

Values  of  es  and  j  determined  by  Collie,  et  al.  (1948)  and 
Hasted  (1961)  are  very  close  to  those  given  in  Table  V. 


Electric  and  magnetic  properties  of  water 


The  permeability  of  water  is  y  (4  x  10  2  hsnrys  per  meter) 

o 


(Lane  and  Saxton,  1951) . 


Equations  (36),  (37)  and  (38)  can  be  used  together  with  (33) 
and  (34)  to  generate  values  of  e'  and  e"  as  a  function  of  T,  S»and 
v.  The  results  are  shown  in  Table  VI. 

To  illustrate  the  variations  of  e'  and  e"  with  temperature, 
salinity,  and  frequency.  Figures  2  and  3,  adapted  from  Saxton 
and  Lane  (1952),  are  presented.  The  salinity  of  'sea  water'  is 
assumed  to  be  36  o/oo.  The  salinity  of  'fresh  water'  is  such 
that  its  ionic  conductivity  is  about  0.01  mho  per  meter.  The 
salinity  of  pure  -ater  is  zero.  Values  of  e'  and  e"  are  plotted 
against  frequncy  for  'pure',  'fresh'  and  'sea  water'  at  Ga  C 
and  20°  C.  Note  that  these  values  are  identical  for  frequencies 
greater  than  about  10  GHz. 

It  will  be  shown  later  that  the  emissivity  and  reflectivity 
of  sea  water  depend  on  its  complex  dielectric  constant,  incidence 
angle,  polarization,  and  the  nature  of  the  shape-  of  the  air-sea 
Interface.  JThus,  it  is  useful  at  this  time  to  make  some  gen¬ 
eral  comments  on  the  behavior  of  the  complex  dielectric  constant 
of  sea  water  as  a  function  of  its  temperature,  salinity,  and 
frequency  of  observation. 

For  all  the  frequencies  chosen,  e'  and  s"  are  dependent  on 
the  temperature  and  salinity  of  the  water.  The  dependence  on 
salinity  is  very  slight  except  for  the  3  GHz  radiation,  e' 
appears  to  be  fairly  steady  with  temperature  at  3  GHz  but  increases 
with  temperature  at  the  other  frequencies,  e"  decreases  with 
temperature  at  3  GHz  and  9.3  GHz,  remains  steady  with  temperature 


Table  VI.  Calculated  values  of  e*  and  e"  for  selected  values  of  temperature  (T,  °C) ,  salinity 
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at  IS. 8  GHz  and  19.35  GHz,  and  increases  with  temperature  at 
22  235  GHz  and  34.0  GHz.  The  marked  change  in  e"  with  salinity 
at  the  warmer  temperatures  at  3  GHz  shows  that  ionic  conductivity 
is  a  dominant  factor  in  (34).  Apparently,  Debye  relaxation 
dominates  at  other  frequencies  (Debye ,  1929). 

The  polynomial  expressions  used  above  are  based  on  measure¬ 
ments  made  for  a  wide  range  of  salinity  and  ccver  salinities 
which  are  much  larger  than  those  observed  In  the  World  Ocean. 

These  measurements  should  be  repeat  J  for  actual  sea  water  for 
many  values  of  salinity  between  0  o/co  and  40  o/oo. 

Electric  and  magnetic  properties  of  ice 

The  permeability  of  ice  is  ^(4*  xl0~?  henrys  per  meter). 
According  to  Gunn  and  East  (1954) ,  the  dielectric  constant 
of  ice  is  not  a  function  of  frequency  in  the  microwave  range. 

Table  VII  shows  the  values  of  s’  and  e"  for  ice  for  various 
temperatures. 

Table  VII.  Dielectric  constant  of  ice  (after  Gunn  and  East,  1954) 
Temperature,  T  (°C) _ ej _ _ 


0 

3.165 

0.00855 

-10 

3.165 

0.00281 

-20 

3.165 

0.00196 

The  permittivity  of  ice  is  assumed  to  be  constant  over  the 
range  of  temperatures  and  frequencies  considered  in  this  paper 
and  equal  to  eQ(3.165  -  j  0.003)  farods  per  meter. 

The  large  difference  between  the  complex  permittivities  of 
water  and  Ice  lead  to  large  differences  between  the  emission  of 
ice  and  water.  It  will  be  sho»:u  later  how  this  can  be  used  to 
detect  ice/ocean  boundaries  and  possibly  to  infer  ice  thickness. 


35 


CHAPTER  III 
RADIATIVE  TRANSFER 


Elementary  Radiation 


Intensity 

In  this  study,  the  word  intensity  denotes  the  specifi-;  nter  - 
sity  of  electromagnetic  radiation.  Consider  the  flow  of  radiant 
energy  as  shown  in  Figure  4.  The  intensity,  I  (v,P,a,t)  is  defined 
by: 

^  Ip  A 

I  (v,P,a,t)  *  dh'dA  dv  dt  (watts/m  sferad ian/Hz) ,  (39) 

where  dE  is?  the  amount  of  radiant  energy  (joules)  that  passes 

2 

through  an  element  of  area  (dA,  m  )  in  an  interval  of  time  (dt,sec). 
This  energy  travels  along  a  direction  indicated  by  the  unit  vector 
a,  and  is  contained  with  an  element  of  solid  angle,  dfi  ,  surround¬ 
ing  a.  This  energy  is  carried  by  electromagnetic  waves  having 
frequencies  in  the  range  from  v  to  v  +vd  (Hertz).  I  is  a  function 
of  frequency,  position,  direction,  and  time.  The  use  of  differen¬ 
tial  notation  in  (39)  implies  that  I  exists  only  if  the  limit 
exists.  For  practical  cases,  one  must  use  finite  areas,  solid 
angles,  times,  and  bandwidths.  The  actual  energy  that  flows  is 
the  integrated  intensity.  It  is  convenient  to  retain  the  concept 
of  a  monochromatic  intensity  since  it  will  be  necessary  to  weigh 
the  intensity  field  by  the  gain  pattern  c'  ’■»  antenna  system  used 
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by  a  given  microwave  radiometer.  Thus,  the  intensity  gives  the 
rate  of  flow  of  energy  per  unit  time  per  unit  solid  angle  per  unit 
area  and  per  unit  frequency  band  width. 


Poynting's  theorem 

The  relationship  between  the  time  varying  electromagnetic 
field  and  the  corresponding  electromagnetic  radiational  intensity 
can  be  deduced  from  Maxwell's  equations  in  a  manner  first  shown 
by  Poynting.  This  development  will  te  outlined  here  using  con¬ 
cepts  established  in  Chapter  II. 

Suppose  one  takes  the  dot  product  of  (8)  and  E,  that  is: 

E  *  V  x  H  =  e  (e*  -  je")  £  •  (watts/m3).  (40) 

O  at 

By  a  vector  identity, 

E-V*3=3.vx£-V*  (£xH)  (watts/m3) .  (41) 


Also, 

_  tf.g) ,  illii ,  a  .  >1 

at  v  ’  at  at 


and 


at 


(H*H) 


-  iifii  =  2ff  .  $ 

at  at 


(volts^/m^  sec), (42) 


(amps^/m^  sec) .  (43) 
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Substitution  of  (41),  (42),  (43),  and  (3)  into  (40)  yields  the 
expression: 

2 

V  •  (E  x  H)  =  -  me  e"E  •  E  -  (e  e' |I|  /2 
0  dt  0 

+  u  |h|2/2)  (watts/m^).  (44) 

o 

The  first  term  on  the  right-hand  side  of  (44)  represents 
the  power  lost  in  the  medium;  the  second  term  represents  the  power 
stored  in  the  electric  and  magnetic  field-  If  one  integrates  over 
a  volume  V  and  uses  the  divergence  theorem,  one  obtains  the  ex¬ 
pression: 

•  #(E  x  H)  •  dA  =  -  ///toe  e”E.E  dV 
A  V 

-  ///  (e  c* |e|2/2  +  u  |h|2/2)  dV  (watts),  (45) 

ot  y  0  O 

where  A  encloses  V  and  dA  is  the  outward  normal  vector  to  the 
surface  A.  From  the  form  of  (45)  one  must  conclude  that  E  x  K 
is  the  power  flowing  across  A  since  energy  must  be  conserved. 

The  first  tena  on  the  right-hand  side  of  (45)  represents  all  the 
losses  occurring  in  V, and  the  second  term  represents  the  time  rate 
of  change  of  the  energy  stored  in  the  electric  and  magnetic  fields 
in  V. 

One  defines  Poynting's  vector  ?  such  that: 

(watts/m2) . 


P  =  E  x  H 


(46) 
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Polarization 

Measurements  of  microwave  intensity  are  accomplished  by 
using  a  radiometer  whose  antenna  is  sensitive  to  only  part  of 
the  radiational  field.  Also,  microwave  power  is  received  and 
integrated  over  periods  of  one  second  or  longer  due  to  the  low 
intensities  involved.  The  intensity  field  consists  of  statisti¬ 
cally  independent  waves.  Microwave  radiometers  are  usually 
sensitive  to  one  plane  of  electrical  polarization.  It  is  con¬ 
venient,  therefore,  to  divide  the  intensity  field  into  two 
components  —  those  waves  whose  electric  vectors  lie  in  the  ver¬ 
tical  plane  and  those  waves  whose  electric  vectors  lie  in  the 
plane  perpendicular  to  the  vertical  plane.  Consider  the  radia- 
tion  traveling  along  the  direction  indicated  by  a.  There  is 

an  infinite  number  of  planes  containing  a.  The  plane  which 
->• 

contains  both  a  and  a^  the  vertical  unit  vector;  is  called  the 

vertical  plane  and  is  denoted  by  the  subscript  v.  The  plane 

*♦- 

containing  a  that  is  perpendicular  to  the  vertical  plane  is 
commonly  called  the  'horizontil  plane',  and  is  denoted  by  the 
subscript  h.  The  horizontal  plane  is  rarely  level,  that  is 
tangent  to  the  geoid;  but,  this  is  the  common  reference  made  to 
this  plane  in  the  literature. 

The  intensity  may  be  thought  of  as  being  the  sum  of  tne 
horizontally  polarized  and  vertically  polarized  components  of 
the  intensity,  that  is. 


40 


I  =  Ih  +  Iy  (watts/m^/steradian/Hz ) .  (47) 

1^  and  are,  in  general,  not  equal.  It  is  assumed  that 
Ih  and  Iv  are  mean  values  integrated  over  some  finite  time  period 
which  is  many  times  longer  than  the  time  period  of  the  radiation. 
Also,  1^  and  I  bear  no  phase  relationship  to  each  other  and 
are  independent  of  each  other.  Thus,  the  time  dependences  in 
(39)  are  only  those  slowly  varying  changes  associated  with  over¬ 
all  changes  in  mean  intensity. 

The  resolution  of  the  intensity  field  into  its  horizontally 
polarized  and  vertically  polarized  components  results  in  the 
necessity  of  expressing  many  of  the  following  equations  as  coup¬ 
lets,  that  is,  one  would  need  one  equation  for  horizontal  po¬ 
larization  and  one  equation  for  vertical  polarization.  To 
reduce  the  number  of  equations  in  the  remaining  pages  of  this 

paper,  the  subscript  p  is  used  when  either  horizontal  or  ver- 

© 

tical  polarization  is  applicable,  that  is,  p  =  h  or  v.  In 
cases  where  cross  polarization  effects  are  under  study,  this  short¬ 
hand  notation  will  be  dropped. 

Blackbody  radiation 

Consider  the  intensity  field  within  a  constant  temperature 
enclosure.  A  constant  temperature  enclosure  is  necessarily  in 
a  state  of  thermodynamic  equilibrium.  Consider  a  small  area  dA 
within  this  enclosure  on  the  walls  of  the  enclosure.  Let 


B(v,T)  be  the  intensity  of  the  radiation  emitted  by  an  area  dA 
where 

dE  =  B(v,T)  dA  dv  dt  dft  (joules).  (48) 

One  major  discovery  of  recent  time  was  the  form  of  B(v,T) 
by  Planck  (1906)  where: 

B( v,T)  -  2hv3/c2  [exp  (hv/kT)  -  1] 

2 

(watts/m  /steradian/Hz) .  (49) 

T  is  the  temperature  (degrees  Kelvin,  °K)  of  the  medium, 

h  is  Planck’s  constant  (6.623  x  10-3^  joule-seconds),  c  is  the 

speed  of  light  in  a  vacuum  (2.9979  x  10®  meters  per  second), 

—23 

and  k  is  Boltzman's  constant  (1.380  x  10  joules  per  degree 
Kelvin) . 

In  the  case  of  frequencies  less  than  100  GHz  and  for  tem¬ 
peratures  more  than  200°  K,  (hv/kT)  is  very  much  less  than 
unity.  A  Maclaurin  expansion  of  the  bracketed  quantity  in  (49) 
is: 

exp  (hv/kT)  -  1  =  hv/kT  +  h2v2/2k2T2  +  .  .  . 

(dimensionless).  (30) 

Equation  (49)  reduces  to 

B(v,T)  =  (2kv2/c2)T  (watts/m2/steradian/Hz)  (51) 

when  higher  order  terms  in  (50)  are  neglected.  Equation  (51)  is 
the  Rayleigh  -  Jeans  Radiation  Law  and  is  valid  for  considerations 
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in  this  paper. 

Suppose  one  removes  the  emitting  matter  from  the  enclosure 
without  changing  its  physical  state.  Will  the  emission  function 
still  be  the  same  as  it  was  in  the  enclosure?  Goody  (1964)  says 
that  it  is  doubtful  that  radiation  laws  valid  inside  a  constant 
temperature  enclosure  would  be  valid  outside  the  enclosure  due 
to  the  existence  of  stimulated  emission.  He  goes  on  to  prove  that 
the  form  of  Planck’s  Law  (and  thus  Rayleigh-Jean's  Law)  is  valid 
as  long  as  local  thermodynamic  equilibrium  is  maintained  or  as 
long  as  the  matter  is  excited  predominantly  through  collisions 
(as  opposed  to  photon  excitment) . 

It  is  assumed  that  Rayleigh-Jean’s  Law  is  valid  for  all 
purposes  in  this  paper. 

A  perfect  radiator  is  called  a  blackbody  since  it  absorbs 
completely  all  radiation  incident  upon  its  surface.  Natural  ob¬ 
jects  may  approximate  blackbodies  at  least  for  some  range  of  fre¬ 
quencies.  The  emission  from  a  blackbody  is  completely  random, 
and;  therefore,  the  polarized  intensity  of  emissions  (Bp,  watts/ 
m^/steradian/Hz)  can  be  expressed  as: 

Bp  =  (kv^/c^)T  (watts /m^/steradian/Hz)  (52) 

since  one-half  the  total  intensity  of  emission  would  exist  in 
each  plane  of  polarization. 

A  natural  body  is  less  efficient  than  a  perfect  emitter 
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(blackbody).  It  is  assumed  that  the  form  of  the  emission  is 
the  same  as  that  of  a  blackbody  or  enclosed  body,  and  that  the 
polarized  intensity  of  emission  is  given  by: 

I6  =  £  B  (watts/ra^/steradian/Hz)  (53) 

P  P  P 


where  ? 

P 


is  the  polarized  emissivity. 


Brightness  Temperature 


In  (52)  there  is  a  linear  relationship  between  the  physical 
temperature  of  a  perfect  emitter  and  its  polarized  intensity  of 
emission-  This  convenient  relationship  led  to  the  common  prac¬ 
tice  in  microwave  physics  of  referring  to  any  given  intensity 
as  an  apparent  temperature  or  brightness  temperature.  This 
practice  results  in  a  simplification  of  resulting  transfer 
equations  and  the  units  used  in  microwave  radiometry. 

All  investigators  have  used  (51)  to  relate  an  intensity 
to  a  brightness  temperature.  Since  microwave  radiometers  are 
sensitive  to  only  one  component  of  tne  natural  field  of  radiation, 
it  is  erroneous  to  employ  a  polarized  brightness  temperature  that 
relates  to  the  total  intensity  of  emission  of  a  blackbody. 

This  incorrect  definition  of  polarized  brightness  temperature  is 
corrected  by  these  same  investigators  by  erroneously  assuming 
that  the  natural  radiation  field  is  unpolarized  and  random. 

The  natural  field  of  radiation  is  usually  biased  toward 
one  particular  plane  of  polarization,  that  is,  the  vertical 
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component  may  be  statistically  more  intense  than  che  horizontal 
component.  One  would  expect,  however,  that  the  emission  from 
a  blackbcdy  is  unooiarized  and  completely  random.  Thus,  one 
should  use  (52)  13  define  polarized  brightness  temperature.  Sup¬ 
pose  the  natural  field  of  microwave  radiation  is  defined  by 
Ip, the  polarized  intensity.  Then  T  ,  the  polarized  brightness 
temperature  (degrees  Kelvin),  is  given  by: 

T  =  (c2/ ko2)  1  (°K).  (54) 

P  P 

It  is  convenient  to  refer  to  certain  components  of  the  total 
polarized  field  of  intensity  as  certain  types  of  polarized  bright¬ 
ness  temperature. 


The  polarized  intensity  (IS^)  of  the  microwave  radiation 

P 

incident  upon  the  solid  or  liquid  surface  of  Earth  is  referred 
to  as  TS^,  the  polarized  sky  temperature  which  is  given  as: 


msky 

V 


(c2/kv2>  lsly 


CK). 


(56) 
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The  following  developments  of  the  equations  of  radiative 
transfer  follow  closely  that  given  by  Chandrasekhar  (1950)  and 
Goody  (1964). 
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Extinction  and  emission 

Radiation  traveling  along  a  particular  path  in  space  and 
existing  in  one  particular  plane  of  polarization  is  called  a 
pencil  of  radiation.  A  pencil  of  radiation  may  be  weakened 
through  extinction  processes  and  strengthened  through  emission 
processes  as  it  trave'/s  through  a  medium.  According  to  Goody 
(1964)  the  fundamental  law  of  extinction  is  Lambert’s  Law  which 
states  that  the  extinction  processes  are  linear  and  independent 
of  the  intensity  of  radiation  and  in  the  amount  of  matter  pro¬ 
vided  that  the  physical  state  is  held  constant.  This  law  is 
assumed  to  be  valid  for  the  purposes  of  this  paper. 

Lambert's  Law  may  be  stated  mathematically  as: 

2 

dl  =  el  de  (watts/m  /steradian/Hz) ,  (60) 

P  P 

where  I  is  the  intensity  of  the  polarized  radiation  under  study, 
dfc  is  the  differential  path  length  (meters) ,  and  e  is  the  volume 
extinction-coefficient  (nepers),  e  is  a  function  of  the  physical 
properties  of  the  matter  and  the  frequency.  It  is  assumed  that 
e  is  isoplaner,  that  is,  it  does  not  depend  upon  the  particular 

plane  of  polarization  under  consideration. 

The  optical  path  length  A (£.,£•' )  is  defined  as: 


A (£,*')  -  / 

Z' 


e(u)du 


(dimensionless) . 


(61) 
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Extinction  consists  of  absorption  (conversion  of  energy  to 
internal  thermal  energy)  and  scattering  (redistribution  of  energy 
with  direction  and/or  frequency).  Since  all  extinction  processes 
are  linear  one  may  say  that 

e  =  a  +  s  (nepers) ,  (62) 

where  a  is  the  volume  absorption-coefficient  (nepers)  and  s  is 
the  volume  scattering-coefficient  (nepers). 

It  is  now  assumed  that  extinction  processes  are  isotropic 
for  a  differential  volume  in  the  atmosphere  or  ocean.  On  the 
other  hand,  the  distribution  of  the  intensity  of  the  scattered 
radiation  is  a  function  of  direction. 

Formally  one  may  state  that  dlp,  the  differential  change  in 
polarized  intensity  due  to  emission  processes,  is  given  by: 

2 

dl  =  eJ  d£.  (watts/m‘Vsteradian/Hz) ,  (63) 

P  P 

where  is  the  polarized  source  function  having  units  of  intensity. 
General  equation  of  transfer 

The  total  differential  change  in  polarized  intensity  is  the 
result  of  the  processes  defined  bv  (60)  and  (63)  and  is  given  by: 

2 

dl  =  (I  -  J  )  dA  (watts/m  /steradian/H z) ,  (64) 

P  P  P 

which  is  known  as  Schwarzschild ' s  Equation  of  Radiative  Transfer. 

The  nature  of  the  polarized  source  function  is  crucial  to 


problems  of  radiative  transfer.  Ihe  simplest  case  is  that  of 
local  thermodynamic  equilibrium.  In  this  case,  the  polarized 
source  function  is  Planckian  (Goody ,  1964)  and  is  given  by: 

2  2  ? 

J  (emission)  =  (kv  /c  )  T  (watts/m“/steradian/Hz) .  '65) 
P 

Conversely,  the  source  of  radiation  may  be  due  to  the  total 
amount  of  scattered  radiation,  that  is, 

Jh  (scattering)  =  /[Thv  (a,  a  )  1^  (ar)  +  Fhh  (a,  a J 
4  r 

Ih  (a^)  1  d:i  2 

(watts/m  /steradian/Hz)  (66) 

and 

Jv  (scattering)  -  /[r^  (a,  a„)  Iy  (afl)  +  ryh  (a,  a  ) 

4tt 

I  (an)]dQ 

(watts/m  /steradian/Hz)  (67) 

where  T,  ,  r. ,  ,  T  .and  F  .  are  the  coefficients  of  scattering 
in  the  medium,  Iv  and  1^  are  the  polarized  intensities  of  the 
incident  field  of  radiation,  and  a0  is  a  unit  vector  which  is 
allowed  to  vary  over  the  entire  hemisphere.  The  integrations  are 
performed  over  4n  steradians.  Now  I  (a^)  and  I  (a  )  are  dif- 

V  P 

ficult  to  assess  in  a  practical  situation  since  they  are  also 
the  result  of  extinction  and  emission  in  the  space  around  the 
differential  volume  under  question.  Most  investigators  avoid 
this  problem  by  assuming  that  scattering  is  negligible.  This 
assumption  is  adopted  in  the  present  paper.  When  cases  arise  for 
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which  this  assumption  seems  to  be  significantly  in  error,  the 
effects  of  scattering  will  be  discussed  in  a  qualitative  manner. 
The  author  intends  to  pursue  the  case  of  a  scattering  and  ther¬ 
mally  emitting  medium  i:i  future  research.  Thus,  with  this  as¬ 
sumption, 

e  =  a  (nepers) f  (68) 

and 

Jp  (scattering)  =  0.  (69) 

Transfer  equation  in  a  non-scattering  medium  under  local  thermo¬ 
dynamic  equilibrium  in  terms  of  brightness  temperatures 

Consider  the  transfer  of  microwave  radiation  that  occurs  in 
a  horizontally  uniform  atmosphere  overlying  a  uniform  ocean  as 
shown  in  Figure  5.  The  polarized  intensities  of  the  microwave 
radiation  incident  upon  the  air-sea  interface  are  the  result  of 
two  processes.  Polarized  cosmic  radiation  is  partially  extincted 
throughout  the  atmosphere  and  reaches  the  surface.  The  atmosphere 
emits  radiation  which  is  partially  extincted  by  the  atmosphere 
between  that  level  and  the  surface.  It  follows  from  equation 
(64)  and  (65)  that  the  polarized  intensity  ISky  of  the  sky  radia¬ 
tion  is  given  by: 

I  ^  “  1  exp  (—A  sectl)  +  (kv^/c^)  seci{>/ 

P  P  m  o 

•  T(A')  exp  ( -A'sect|))dA  (watts/m^/steradian/Hz)  ,  (70) 


adiatlve  transfer 
le,  narallel  atroosn 


where: 


A  =  /  a(z)  dz 
m  o 

(.dimensionless)  , 

(71) 

A ’  =  /  a(z)  dz 

0 

(dimensionless) , 

(72) 

dA'  =  a(z)  dz 

(dimensionless) . 

(73) 

z  is  the  distance  (meters)  above  the  air-sea  interface. 

If  the  intensities  are  expressed  as  polarized  brightness 
temperatures  according  to  (55)  and  (56),  then: 

Tsky  =  TC  exp(-A  secili)  +  sec*A»  T(A’)  exp(-A,Recv!»)dA' 
p  p  m  o 

(°K) .  (74) 

In  most  cases  the  contribution  from  cosmic  sources  can  be 
ignored  when  compared  to  emission  from  the  atmosphere.  Solar 
cosmic  radiation  cannot  be  ignored  since  the  sun  radiates  with  an 
effective  temperature  of  11,000°  K  in  the  microwave  region  (Barrett 
and  Chung »  1962).  The  attenuation  of  solar  microwave  radiation  by 
atmospheric  gases  has  been  measured  by  Staelin  (1966) ,  Dicke  (1946) , 
and  other  investigators  to  determine  the  distribution  of  water 
vapor,  oxygen,  and  temperature  in  the  atmosphere. 

Extinction  and  emission  at  the  air-sea  interface 

At  the  air-sea  interface  some  microwave  radiation  is  reflected 
or  scattered  into  the  atmosphere  and  the  rest  enters  the  hydro¬ 
sphere.  In  the  hydrosphere  additional  extinction  takes  place. 
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Unless  the  depth  of  water  is  extremely  small,  all  of  the  microwave 
energy  entering  the  hydrosphere  is  absorbed. 

If  the  air-sea  interface  is  flat,  the  reflected  pencil  of 
radiation  will  emerge  along  a  direction  inclined  at  *  which 
is  equal  to  i|i  and  is  opposite  *o  \j,  (see  Figure  5) .  The  incident 
and  reflected  pencils  lie  in  a  vertical  plane.  If  the  air-sea 
interface  is  smooth  but  is  not  flat  the  reflected  pencil  lies 
along  an  angle  coangular  to  the  incident  angle  with  respect  to 
the  normal  of  the  inclined  smooth  surface.  A  surface  is  defined 
to  be  smooth  if  the  radius  of  curvature  of  the  surface  at  a 
point  is  many  times  larger  than  the  wavelength  of  the  micro- 
wave  radiation  (Peake,  1959).  Since  microwave  wavelengths 
are  typically  one  to  one-hundred  centimeters,  the  assumption 
applies  only  to  long  wavelength  ocean  swells.  In  most  cases,  the 
sea  surface  is  rough  and  must  be  treated  from  the  standpoint  of 
physical  optics. 

Peake  (1959)  has  defined  a  set  of  quantities  known  as  the 

bistatic  scattering  coefficients  per  unit  area  which  give  the 

scattering  properties  of  an  interface  separating  two  homogenous 

dielectrics.  Let  1^  and  I*  be  the  polarized  intensities  of 

radiation  incident  upon  the  scattering  surface  along  a  direction 

s  s 

indicated  by  a^, a  unit  vector.  Let  I  and  be  the  polarized 
intensities  of  radiation  scattered  along  a  direction  indicated 
by  ag  c?om  the  scattering  interface.  Let  D  be  the  radial  dis¬ 
tance  from  the  point  of  observation  and  A  be  the  area  of  the 
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scatterer.  The  scattering  coefficients  per  unit  area  are  defined 
by  Peake  (refer  to  Figure  6)  as: 


Yhh 

lim 
=  A-*“ 
R-x* 

AttD2 

g 

I.  /Acosii'. 
h  i 

(®~2) • 

(75) 

Yhv 

lim 
=  A-x» 
R-ko 

4ttD2 

I^/Acos’J'i 

(m"2) . 

(76) 

Yw 

lim 
*=  A-x» 

R-x» 

AtiD^ 

Is/Acos’4'. 
v  i 

(m‘2) . 

(77) 

Yvh 

lim 
*  A-*oo 
R-x» 

4*D2 

I®/Acos^ 

(m-2) . 

(76) 

where  ij^is  the  angle  between  and  a . 

Consider  the  transfer  of  radiation  through  a  homogenous  slab 
of  dielectric  material.  The  polarized  intensity  of  the  radiation 
entering  the  material  is  less  than  I*, the  polarized  intensity  of 
the  radiation  incident  upon  the  slab.  If  the  surface  of  the  slab 
is  a  smooth  surface,  then  1°,  the  polarized  intensity  of  the  radia¬ 
tion  entering  the  slab* is  given  by: 


1°  =  (1  -  r  )  I1 
P  P  P 


(watts/m  /steradiar.  )  (79) 


where  r^  is  the  polarized  reflectivity.  If  the  surface  is  a  scat¬ 
tering  surface,  1°  is  dependent  upon  the  distribution  of  y’s. 
According  to  Peake,  the  polarized  albedos  of  a  surface  per  unit 
solid  angle  is  given  by: 
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gh  =  (l/4ir)  /[Yhh  (Si,  t)  +  Yhv(at,  t))dh  (m“2),  (80) 

2it 

and 

8  v  =  (l/4n)  /[Yvh  (ti,  t)  +  Y^C^.^JdA  (m'2),  (81) 

2x 

where  g^  and  g^  are  the  polarized  albedos. 

Thus, 

1°  =  (1  -  3^)  1^  (watts/m2/steradian/Hz) .  (82) 

The  entering  radiation  will  suffer  some  attenuation  as  it 
travels  in  the  dielectric  slab  and  will  gain  some  energy  due  to 
emissions  in  the  slab.  The  polarized  intensity  of  the 
resulting  transmitted  energy  is 

I*  -  1°  exp  (-  /Xadx)  +  (kv2/c^)/X  aT 
p  p  o  o 

•exp  (-/Xadx)dx  (watts/m2/steradian/Hz) .  (83) 

o 

therefore, 

I*  *  1°  exp  (-  ax)  +  (kv2/c2)T  [1  -  exp  (-  ax)], 

P  P 

2 

(watts/m  /steradian/Hz) ,  (84) 

where  a  is  the  absorption  coefficient  of  the  slab,  T  is  the  slab 
temperature,  and  x  is  the  slab  thickness.  The  term,  exp  (-ax  )» 
is  the  polarized  transmissivity  of  the  medium.  The  amount  of 
radiation  absorbed  is  1°  [1  -  exp  (-  ax)]  since  energy  must  be 
conserved.  The  second  term  of  the  righthand  side  of  (.84)  is 
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T(1  -  exp  (-  ax)l»  a^d  this  form  suggests  that  the  bracketed 
multiplier  takes  on  the  role  of  a  polarized  emissivity  for  the 
entire  slab;  however,  the  same  term  gives  us  the  expression  for 
the  absorptivity.  This  equivalence  of  absorbtivity  and  emissivity 
is  the  result  of  the  assumption  that  local  thermodynamic  equilib¬ 
rium  is  maintained  as  was  first  stated  by  Kirchhoff.  If x  is 
increased  to  such  an  extent  that  the  transmissivity  is  practi¬ 
cally  zero,  then  one  may  draw  certain  conclusions  about  the 
relationships  between  the  scattering,  absorbing  and  emissive 
properties  of  the  air-sea  interface,  and  these  are  as  follows: 

Tor  a  smooth  surface: 

*  1  -  r  (dimensionless) ,  (85) 

and  for  a  rough  surface: 


(dimensionless) . 


(86) 


Consider  the  intensity  of  the  upwelling  microwave  radiation 
at  some  point  A"  above  the  surface  of  the  sea.  The  expression  for 
Tp, polarized  brightness  temperature, in  a  horizontally  uniform 
atmosphere  is 

T  (A",e)  *  {(1  -  r  )  T  +  r  [secV*  T(A’) 

p  p  bca  P  O 

•  exp  (-sec^A^dA*  +  TC  exp  (-sec^A  )]) 

p  m 

•  exp  (-sec^A”)  +  sec T(A,M)  exp  (-secif>A,M)dA"' 

°  (°K),  (87) 


where: 


dO 

=  I  a(z)dz 

0 

(dimensionless) , 

(88) 

=  /Z  u(z)dz 

0 

(dimensionless) , 

(89) 

=  I**0  a(z)dz 
o 

(dimensionless) , 

(90) 

=  I**0  a(z)dz 

7 

(dimensionless) . 

(91) 

In  the  above  equation,  hQ  is  the  height  of  the  observer 
above  sea  level  in  meters  and  z  is  the  height  of  a  level  under 
consideration.  The  computations  in  (87)  through  (91)  can  be 
performed  by  various  finite  difference  schemes  on  the  electronic 
computer  if  d(z)  and  T(z)  are  known.  a(z)  is  usually  dependent 
upon  temperature,  density,  pressure,  and  type  of  absorbing  medium. 

In  the  microwave  region  the  atmospheric  gases  of  molecular 
oxygen  and  water  vapor  and  the  atmospheric  constituents  of  liquid 
water  droplets  and  ice  are  the  principal  absorbers.  Their  rela¬ 
tive  and  absolute  values  are  the  subject  of  Chapter  IV. 


CHAPTER  IV 


JO 


EXTINCTION  AND  EMISSION  IN  THE  ATMOSPHERE 

General 

Atmospheric  gases  and  other  atmospheric  constituents  selec¬ 
tively  absorb,  scatter,  and  emit  microwave  radiation  as  it  travels 
through  the  atmosphere.  Molecular  •  xygen  and  water  vapor  are  the 
primary  absorbing  gases  in  the  atmosphere  in  the  microwave  region 
of  the  electromagnetic  spectrum  (Van  Vleck,  1951a) .  Droplets  of 
liquid  water  in  clouds  and  rain  strongly  absorb  and  scatter  mi¬ 
crowave  radiation.  As  a  consequence  of  Kirchhoff's  Law,  n.clecu- 
lar  oxygen,  water  vapor,  and  droplets  of  liquid  water  are  the 
primary  sources  of  emitted  microwave  energy  in  the  atmosphere. 

All  investigators  in  microwave  radiometry  have  assumed  that  scat¬ 
tering  is  negligible  as  an  extinction  process  as  well  as  an  emis¬ 
sion  process.  This  assumption  greatly  simplifies  the  solutions 
to  radiative  transfer  since  it  is  necessary  to  know  the  entire 
field  of  intensity  throughout  the  volume  being  considered  if 
one  is  to  treat  a  scattering  atmosphere. 

Since  this  study  is  a  review  of  the  present  state  of  investi¬ 
gations  in  microwave  radiometry,  qualitative  speculation  of  the 
effects  of  scattering  under  actual  atmospheric  conditions  will  be 
offered  along  with  quantitative  results  from  other  investigations 
that  treat  the  atmosphere  as  a  non-scattering  medium  in  local 


thermodynamic  equilibrium. 
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Previous  studies  of  the  absorption  coefficient  and  scattering 
coefficient  of  molecular  oxygen,  water  vapor,  and  liquid  water  were 
often  presented  from  the  standpoint  of  applications  to  radar.  Con¬ 
sequently,  the  total  attenuation  coefficient  or  extinction  coef¬ 
ficient  is  usually  discussed  in  these  papers.  Thermal  emission  of 
rain  clouds  is  of  little  concern  to  radar  studies  since  the  inten¬ 
sities  involved  in  radar  are  several  orders  of  magnitude  larger 
than  the  intensity  of  thermal  emission.  For  the  purposes  of  the 
present  paper,  a  knowledge  of  the  absorption  characteristics  of 
the  atmosphere  is  needed  as  opposed  to  the  total  attenuation 
characteristics. 

One  other  point  of  semantics  needs  clarification  before  going 
into  detailed  discussion  of  absorption  and  emission  in  the  atmo¬ 
sphere.  Most  investigators  have  used  absorption  coefficients 
which  are  expressed  in  terms  of  decibels  per  kilometer"*  (db/km). 

In  Chapter  III,  a  is  expressed  in  nepers.  Let  1°  be  the  polarized 
intensity  of  a  pencil  of  radiation  at  some  point  along  its  path, 
x  =  0,  and  let  1^  be  the  polarized  intensity  of  this  pencil  at  x. 
Let  be  the  absorption  coefficient  along  x  in  terms  of  decibels 
per  kilometers,  and  let  a  be  the  absorption  coefficient  along  x  in 


terms  of  nepers.  Then, 

IX/I°  =  exp  (-ox) 

P  P 


(dimensionless) ,  (92) 


Lengths  may  be  expressed  as  centimeters  (cm),  meters(m)  or 
kilometers  (km)  in  various  expressions  in  the  following  develop¬ 
ments  . 


r 


-10  log10  (I*  /  1°)  -  10-3  adb  X  (db) .  (93) 

Equation  (92)  follows  from  (60),  and  (93)  follows  from  the 
definition  of  a  decibel  where  x  is  in  meters  in  both  cases. 

From  (92)  it  follows  that: 

-in  (IX  /  1°)  =  ax  (dimensionless) ,  (94) 

P  P 

or 

-10  log..  (IX  /  1°)  *  4.3429  ax  (db) .  (95) 

10  p  p 

Thus,  from  (93)  and  (95),  it  follows  that: 

odb  -  4342.9a  (db/km).  (96) 

Conversely, 

a «  0.00023026  a,,  (nepers).  (97) 

ab 

If  Is  expressed  as  decibels  per  meter,  then 
ab 

a*  0.23026a.,  (nepers).  (98) 

ab 

Atmospheric  Gases 


Molecular  oxygen 

The  oxygen  molecule  (0^)  absorbs  microwaves  because  it  pos¬ 
sesses  a  permanent  magnetic  moment.  Van  Vleck  (1947a)  made 


predictions  of  the  absorptive  characteristics  of  molecular  oxygen 
from  quantum  theory.  As  a  result  of  collisions  with  neighboring 
molecules,  the  energy  levels  of  the  oxygen  molecule  are  broadened, 
and  a  statistical  distribution  of  allowable  transitions  results. 
Van  Vleck  used  the  broadening  distribution  developed  by  Van  Vleck 
and  Weisskopf  (1945). 

The  general  quantum  mechanical  expression  for  a  ,  the  ab¬ 
sorption  coefficient  of  molecular  oxygen,  is  quite  complicated. 
Oxygen  has  a  series  of  absorption  lines  near  60  GHz.  Outside  this 
resonance  region,  one  may  assume  that  one  central  resonance  line 
exists  and  may  predict  absorption  based  on  broadening  of  this 

line.  Van  Vleck  (1947a)  gives  a  simplified  expression  for  a 

-  '  ox 

for  a  temperature  of  293°  K  and  for  a  pressure  of  1013.25  milli¬ 
bars  (mb)**  that  is  applicable  for  frequencies  less  than  35  GHz; 
this  expression  is: 

a  =0.34  (v2/c2)  fpr - ,  /tn 

ox  j_[2  -  (v/c)j^  +  (Av/c)4 

_ A\>/ c  _  _ Av/c _ 7|(db/kn,  (99) 

+  [2  +  (v/c)]^  +  (Av/c)^  (v/c)^  +  (Av/c)^J 


where : 

c  =  speed  of  light  (cm/sec), 
v  =  frequency  (GHz) , 

(Av/C)  =  line  breadth  constant  (cm-*). 

^In  the  rationalized  MKS  system  of  units  1  millibar  (mb)  = 
100  newtons  per  square  meter  (nt/mz). 
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The  value  of  (dv/c)  has  to  be  determined  from  experimental 

data  since  the  nature  of  the  cross  section  for  collisions  between 

atmospheric  gases  is  poorly  understood.  Van  Vleck  selected  a 

value  for  (dv/c)  of  0.02  cm“^  and  made  calculations  of  a  based 

ox 

on  equation  (99) .  These  values  are  given  in  Table  VIII  along  with 
similar  values  based  on  the  quantum  expression  for  for  fre¬ 
quencies  above  35  GHz. 

Table  VIII.  Theoretical  values  of  the  absorption  coefficient  of 
molecular  oxygen  (aQX,  db/km)  at  a  pressure  of 
1013.25  mb  and  a  temperature  of  293.0°  K  (after 
Van  Vleck,  1951). 


V 

(GHz) 

a 

ox 

(db/km) 

V 

(GHz) 

aox 

(db/km) 

V 

(GHz) 

aox 

(db/km) 

0.3 

0.0014 

30 

0.014 

64.5 

5.0 

1.0 

0.0050 

45 

0.077 

69 

0.51 

3.0 

0.0066 

51 

0.32 

75 

0.19 

10 

0.0072 

54 

1.99 

150 

0.03 

20 

0.0089 

60 

14.0 

The  values  of  a  given  in  Table  VIII  are  for  standard  condi- 
tions  (p  *  1013.25  mb;  T  ■  293°  K).  These  values  served  a  useful 
purpose  as  far  as  horizontal  propagation  of  microwave  energy  in 
the  atmosplere  is  concerned.  Vertical  and  slant-range  propagation 
is  of  concern  in  this  study.  For  other  pressures  and  other  tem¬ 
peratures,  the  expression  given  in  (99)  is  altered  in  two  ways. 
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First,  the  line-breadth  constant  is  affected  by  temperature  and 
pressure.  This  dependence  is  not  veil  established.  Van  Vleck 
(1947a)  states  that: 

(Av/c)  «  pT‘1/2.  (100) 

Meeks  and  Lilley  (1963)  stated  that  measurements  by  Hill  and 
Gordy  (1954)  indicate  that: 

(Av/c)  *  pT-0*85.  O01) 

Falcone  (1966)  confirmed  (101)  through  measurements  of  sky 
radiance  where  the  value  of  (Av/c)  at  1013.25  mb  and  293°  K  was 
taken  to  be  0.025  cm-'1'.  There  may  be  some  non-linearity  in  (Av/c) 
in  regard  to  its  dep  idence  on  p.  Zimmer  and  Mizushlma  (1961) 
report  that  the  slope  of  Av  versus  p  is  2  MHz  per  mm  of  mercury 
at  very  low  pressures  and  is  0.8  MHz  per  mm  of  mercury  at  one 
atmosphere. 

Broadening  also  results  from  Zeeman  splitting  and  thermal 
agitation  (Doppler  broadening).  Above  40  kilometers,  the  Zeeman 
effect  becomes  a  major  factor;  above  80  kilometers,  Doppler  broad¬ 
ening  is  the  main  cause  of  broadening  (Meeks  and  Lilley ,  1963) . 

The  effect  of  these  causes  can  be  accounted  for  by  adding  a  line 
breadth  constant  for  Doppler  (Ziw^c)  and  Zeeman  (Av^/c)  effects. 

In  addition  to  the  effects  of  pressure  and  temperature  on 
the  line-breadth  constant  in  (99),  the  constant  term  of  (99)  is 
affected  by  pressure  and  temperature. 


. a - -—  . . . M . . ^-u..,^.,-.t....^.v.il,.r.......um«ta»i.  . . . 
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Water  vapor 

The  molecule  of  water  (H-O-H)  in  its  gaseous  state  absorbs 

microwaves  due  to  its  permanent  electric  dipole.  Van  Vleck  (1947L) 

determined  its  absorptive  characteristics  from  quantum  theory  in 

a  manner  similar  to  that  used  to  determine  a  .  In  the  case  of 

ox 

water  vapor,  there  is  a  single  absorption  line  at  22.235  GHz. 

Since  water  vapor  is  electrically  polar,  strong  absorption  occurs 
at  infrared  frequencies.  Broadening  of  infrared  absorption  lines 
results  in  some  additional  absorption  in  the  microwave  region 
which  is  called  residual  absorption.  Becker  and  Autler  (1946) 
have  shown  that  the  contribution  to  a^,  the  absorption  coefficient 
for  water  vapor,  due  to  residual  absorption  was  about  five  times 
larger  than  that  predicted  by  Van  Vleck.  The  quantum  mechanical 
expression  for  a w  developed  by  Van  Vleck  can  be  simplified  and 
corrected  for  the  observed  residual  absorption.  The  resulting 
expression  given  by  Staelin  (1966)  is: 

a  =  140.71  exp(-644/T)v2pp  T-3,125  (1  -  0.0147  p  T/p) 

W  W  w 

-  U/l(v  -  22. 234)2  +  (Av)2]  +  l/[(v  +  22.234)2  +  (Av)2} } 

+  0.01107  p  2AvT-1*5  (db/km),  (102) 

w 

where : 

v  =  frequency  (GHz) 

T  =  temperature  (°K) 


P  * 


65 


p  =  total  pressure  (mb) 

3 

Pw  =  density  of  water  vapor  (g/m  ) 
and 

Av  =  2.58  x  10'3  p(l  +  0.0147p  T/p)  (T/318)'0,625  (cm-1). 

w 

Van  Vleck  (1951)  has  calculated  a  at  1013.25  mb  pressure 
”  w 

and  for  a  temperature  of  293°  K.  These  calculations  are  given  in 
Table  IX. 

Total  gaseous  absorption 

Dicke,  et  al.,  (1946)  determined  a  through  measurements  of 

'  V 

sky  radiance  at  20  GHz,  24  GHz, and  30  GHz.  These  values  are  0.014, 

3 

0.026, and  0.011  (db/km  per  g/m  ),  respectively,  and  they  agree 

closely  with  the  corrected  theoretical  values  given  in  Table  IX. 

The  density  of  water  vapor  depends  somewhat  on  the  temperature 

of  the  air  and  varies  over  a  large  range  throughout  the  year. 

At  sea  level,  saturated  air  at  a  temperature  of  20°  C  contains 

about  17  grams  of  water  vapor  per  cubic  meter  (Chemical  Rubber  Co., 

3 

1959).  Thus,  o  may  be  17  g/m  or  less  under  these  conditions. 

V 

In  model  atmospheres,  it  has  been  popular  to  assume  that  is 
7.5  grams  per  cubic  meter  at  sea  level.  In  Table  X,  values  of 
a__.  and  a  have  been  calculated  for  a  pressure  of  1013.25  mb,  for 

u*  V 

3 

a  temperature  of  293°  K,  and  p  equal  to  3,  7.5,  and  17  g/m  .  It 

w 

can  be  seen  in  Figure  7  that  water  vapor  is  the  principal  absorber 
for  frequencies  between  8  GHz  and  50  GHz.  Water  vapor  is  quite 


*  Values  based  on  theory 
**  Theoretical  values  times  five 
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FREQUENCY,  v(GHz) 

Figure  7.  The  absorption  coefficient  of  molecular  oxygen  and 
water  vapor  at  sea  level  versus  frequency  for  a 
temperature  of  293°  K  and  for  water-vapor  densities 
of  0,  3,  7.5,  and  17  grams  per  cubic  meter. 
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Table  X.  Total  absorption  coefficient  for  molecular  oxygen  and 
water  vapor  (a  ,  db/km)  for  a  pressure  of  1013.25  mb, 
for  a  temperature  of  293°  K  and  for  water-vapor 
densities  of  3,  7.5  and  17  g/m3. 


Frequency 

(GHz) 

(8i>/km) 

p'  =  3 
(g/m3) 

p  =  VV5 
(g/m3) 

- - - - 

1 V  (g/m3) 

3 

0.0066 

0.0068 

0.0072 

0.0079 

10 

0.0072 

0.0075 

0.014 

.0226 

15 

0.0078* 

0.0162 

0.029 

.0556 

20 

0.0089 

0.047 

0.104 

.2248 

21 

0.0092 

0.063 

0.141 

.308 

22.235 

0.0096* 

0.077 

0.178 

.392 

24 

0.0102 

0.072 

0.164 

.359 

30 

0.014 

0.046 

0.125 

.193 

45 

0.077 

0.121 

0.187 

.325 

*  Value  is  estimated  from  plot  of  a 

variable  in  the  atmosphere  in  amount  and  distribution  with  height. 

In  most  cases  the  amount  and  distribution  of  water  vapor  distribution 

tion  in  the  atmosphere  is  not  well  known.  Thus,  one  would  endeavor 

to  use  frequencies  less  than  8  GHz  if  one  wishes  to  view  surface 

phenomena.  However,  if  one  is  interested  in  measuring  remotely 

some  aspect  of  the  field  of  water  vapor,  then  frequencies  near 

22.235  GHz  seem  appropriate. 

Tables  VIII  through  X  and  Figure  7  present  data  for  sea 

level  pressure  and  temperature  only.  Since  the  transmission  of 

microwave  energy  through  the  atmosphere  along  vertical  and  slanting 

paths  is  of  concern  to  this  study,  the  eftects  of  pressure  and 

temperature  on  a  must  be  considered .  Perhaps  the  best  way  to 
8 

proceed  now  is  to  examine  the  results  of  past  calculations  of  sky 
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radiance  and  total  atmospheric  transmission  *or  a  cloudless  atmo¬ 
sphere.  Most  investigators  have  calculated  sky  radiance  or  total 
attenuation  by  assuming  a  model  atmosphere  and  by  using  Van  Vleck's 

formulas  for  a _ and  a  . 

ox  w 

Total  atmospheric  attenuation  for  a  clear  atmosphere 

If  o(z)  is  known,  then  L^,,  the  total  zenith  attenaution  of 
microwave  radiation,  is: 

L  «  /°V(z)dz  (db),  (103) 

g  o  g 

where: 

L  =  total  zenith  attenuation  (db) , 

O 

a  -  attenuation  coefficient  for  molecular  oxygen  and  water 
°  vapor  (db/km), 

and 

z  =  altitude  (km). 

Equation  (103)  may  be  evaluated  by  finite  differences  or  any 
other  approximation  scheme. 

Using  the  distribution  of  pressure  and  temperature  given  by 

the  1960  ARDC  model  atmosphere.  Keeks  and  Lllley  (1963)  calculated 

L  for  oxygen  only  over  the  range  of  frequencies  from  1  GHz  to 
g 

120  GHz.  In  the  range  from  1  GHz  to  10  GHz,  h  was  approximately 

g 

0.04  db. 


Stogryn  (1967)  calculated  values  of  the  fractional  absorption 
of  a  1960  ARDC  model  atmosphere  for  various  values  of  altitude  and 


frequencies  of  x9.4  GHz  and  35  GHz.  He  assumed  that  the  vertical 
distribution  of  water  vapor  was  given  by: 

p  (z)  *  9.35  exp  (-0.5z)  (g/m^) ,  (104) 

w 

where  z  is  the  altitude  (km). 

In  terms  of  decibel  loss,  the  attenuation  due  to  absorption 
by  oxygen  and  water  vapor  in  a  layer  of  atmosphere  extending  from 
the  surface  to  altitude  z(km)  is: 

L  (z)  «  4.3429  /Z  o(z)  dz  (db) ,  (105) 

g  o 

where  a(z)  is  in  nepers  and  z  is  in  meters.  Stogryn  gives  values 
for  the  integral  in  (105).  Table  XI  gives  these  values  of  L (z) 

O 

expressed  in  decibels.  The  majority  of  the  loss  is  incurred  in 

the  lower  five  kilometers  of  the  atmosphere.  This  emphasizes  the 

dependence  of  a m  and  aw  on  pressure.  The  total  loss  throughout 

the  entire  vertical  path  is  0.346  and  0.326  decibels  at  19.4  GHz 

and  35  GHz,  respectively.  This  is  approximately  one  order  of 

magnitude  higher  than  the  loss  at  frequencies  less  than  10  GHz. 

Calculations  of  L  have  also  been  made  by  Peake  (1967)  and 
g 

Harandino  (1967) .  These  publications  did  not  include  information 
about  the  exact  atmospheric  model  used.  Wulfsberg  (1967)  measured 
the  total  attenuation  of  solar  microwave  energy  by  the  atmosphere 
over  a  six-raor.th  period  at  15  GHz  and  35  GHz. 

Hogg  (1959)  computed  values  of  based  on  the  distribution 
of  temperature,  pressure  and  water-vapor  given  by  the  International 


Table  XI.  Zenith  attenuation  (L  (z),  db)  in  a  1960  ARDC  model  atmosphere  with  a  density  of  water 
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Standard  Atmosphere,  that  is,  he  assumed  that: 

T(z)  =  255  +5.5  (10  -  z)  (°K),  (106) 

p  (z)  =  2  (5  -  z)  (g/m3),  (107) 

w 

for  (0  <_  z  <_  5  km), 

p  (z)  =  0  (g/m3)  for  (z  >_  5  km),  (108) 

w 

and 

p(z)  -1039.91  (-0.014z)  (mb),  (109) 

where  z  is  the  altitude  in  kilometers. 

The  results  of  these  calculations  and  observations  of  L 

g 

are  presented  in  Taole  XII. 

The  calculated  values  in  Table  XII  agree  among  themselves  quite 

well  and  agree  closely  to  the  mean  total  attenuation  observed  by 

Wulfsbsrg  (1967).  His  observations  also  show  that  L  is  very  sensi- 

g 

tive  to  the  distribution  of  pw  for  frequencies  near  the  resonance 
absorption  line  for  water  vapor.  The  values  of  for  frequencies 
less  than  10  GHz  are  quite  low  and  indicate  that  radiation  at  these 
frequencies  would  suffer  a  loss  of  about  one  percent  when  traveling 
either  outward  from  the  surface  or  vice  versa. 

Sky  temperature 

If  o(z)  and  T(z)  are  known,  then  TSj^,  the  polarized  sky 
temperature,  may  be  calculated  from  (74).  Several  investigators 
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have  calculated  TSp^  under  specific  conditions  of  assumed  model 
atmospheres.  Several  other  investigators  hnv«  measured  the  sky 
radiation  with  ground-based  microwave  radiometers.  The  results 
of  these  calculations  and  observations  will  bo  reviewed  row  i.n 
order  to  understand  better  the  magnitude  and  variability  of  this 
parameter . 

Hogg  (1959)  calculated  values  of  TSp^,  under  tKe  conditions 
given  in  (106)  -  (109)  for  a  stratified  atmsopliere 

Weger  (1960)  calculated  values  of  TSp'  under  .rightly  dif¬ 
ferent  assumptions.  He  assumed  that  the  atmosphere  was  horizontally 
uniform  and  that  the  surface  of  Earth  was  flat.  He  considered 
absorption  above  the  tropopause  to  be  negligible.  Using  Van 

Vleck's  calculations  of  ag  for  a  pressure  of  1013.25  mb,  for  a 

3 

temperature  of  293°  K  and  for  equal  to  7.5  g/m  ,  he  formulated 

expressions  for  a  and  a  as  follows: 
ax  w 

a  =  a  (1013.25  mb,  293°  K,v)  exp  (-0.183z)  (km  *)  (110) 

ax  ox 


and 


aw  =  aw  (1013.25  mb,  293°  K,  7.5  g/m  ,v) 


exp  (-0.5z) 


(km  \  (111) 


where  z  is  the  altitude  in  kilometers.  The  surface  values  of  a 

ox 

and  a  may  be  taken  from  Van  Vleck  (1951a). 

slcv 

Croom  (1965)  calculated  values  of  T  by  assuming  an  I.C.A.Q, 
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Standard  Atmosphere  with: 

pw  =  10  exp  (-0.58z)  (g/m3),  (112) 

where  z  is  the  altitude  in  kilometers  (0  <  z  <  20  km) .  In  addition 
to  this  distribution  of  water  vapor,  he  assumed  that  a  layer  of 
water  vapor  existed  in  the  stratosphere  which  was  5  km  thick  and 
which  contained  0.01  grams  of  water  \  apor  per  cubic  meter.  The 
result  of  stratospheric  water  is  most  pronounced  at  22.235  GHz. 

Stogryn  (1967)  calculated  values  of  at  19. A  GHz  and 

-  p 

35  GHz  based  on  his  model  of  the  atmosphere  previously  give... 

gl  y 

Hogg  and  Semplak  (1963a)  have  calculated  values  of  t 
based  on  "dry,  average  and  humid”  weather.  These  values  are 
given  by  Peake  (1967)  without  explanation  of  the  precise  model 
employed  by  Hogg  and  Semplak.  Marandino  (1967)  also  gives  values 
of  TF,Cy  for  15.8  GHz,  22.2  GHz,  and  3A  GHz. 

p 

Wulfsberg  (1964)  has  approached  the  problem  of  calculating 
TSp7  from  an  engineering  and  empirical  point  of  view.  If  one 
considers  the  atmosphere  as  a  whole,  then  i ^  the  transmissivity,  is 
given  by: 

t  -  exp  (-seco;  fm  a(z)  dz)  (dimensionless).  (113) 

O  0 

If  one  assumes  furthermore  that  the  atmosphere  is  non¬ 
reflecting  then  the  absorptivity,  aQ,  is  given  as: 


(dimensionless) .  (HA) 
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By  Kirchhoff's  Laws,  the  emissivity  is  equal  to  the  absorp¬ 
tivity  and,  thus, the  sky  temperature  is  given  by 

TSpY  =  Tm  [1  "  exp  <°K),  0i5) 

where 

oo 

a  =  /  a(z)dz  (dimensionless)  (116) 

o  o 

and  =  effective  temperature  (°K)  of  the  atmosphere.  Without 

proof,  Wulfsberg  (1964)  and  Peake  (1967)  have  given  the  following 

empirical  expression  for  T  : 

m 

T  =  1.12  T  -  50  (°K) ,  (117) 

m  g 

where  T  =  the  temo  rrature  of  the  air  at  the  surface  (°K) . 
g 

Kreiss  (1968)  assumed  that  T  could  be  given  by  the  hypso- 
-  m 

sky 

metric  equation  for  a  layer,  and  he  calculated  values  of  T  ^ 

by  dividing  the  atmosphere  into  fourteen  layers  and  applying  an 

equation  similar  to  (115)  to  each  layer  and  summing  over  all  the 

layers  accordingly.  Wulfsberg  (1967)  stated  that  (117)  could  be 

used  in  place  of  a  more  rigorous  expression  to  give  correct 

values  of  T  to  within  2°  K. 
m 

Conway  and  Sakamoto  (1965),  Cummings  and  Hull  (1966),  and 
Decker  and  Dutton  (1968)  have  measured  Ts3c^  of  the  real  atmo- 

-  -  p 

sphere  at  various  frequencies. 


Table  XIII  presents  a  summary  of  some  of  the  calculated  and 
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measured  values  of  polarized  sky  temperature  determined  by  the 
aforementioned  investigators  for  various  zenith  angles. 

It  will  be  shown  later  that  approximately  two-thirds  of  the 
microwave  radiation  incident  upon  the  sea  surface  is  reflected 
or  scattered  back  into  the  atmosphere.  Also,  the  intensities 
measured  by  a  downward  looking  radiometer  are  typically  200°  K 
in  terms  of  brightness  temperature.  From  Table  XIII,  it  is 

slvV 

seen  that  T  •  is  usually  less  than  30°  K  for  frequencies  below 
about  16  GHz.  Thus,  approximately  20°  K  or  less  of  the  total 
brightness  temperature  is  due  to  reflected  sky  radiation  for 
frequencies  less  than  16  GHz.  In  most  cases  this  contribution 
to  the  total  radiation  is  quite  small  for  these  low  frequencies. 
Also,  the  effect  of  natural  variations  in  water  vapor  is  very- 
small  for  these  frequencies.  However,  the  amount  of  reflected 
sk--  radiation  is  quite  large  for  frequencies  above  16  GHz.  The 
calculations  by  Kreiss  (1968)  illustrate  the  pronounced  effect 
of  water  vapor  on  sky  temperature  at  frequencies  near  the  absorp¬ 
tion  line  of  water-vapor.  The  natural  variability  of  water  vapor 
in  the  atmosphere  can  cause  variations  in  the  microwave  radiance 
of  the  sea  of  10  to  30°  K  for  frequencies  ranging  from  19  GHz  to 
25  GHz. 

If  the  atmosphere  is  horizontally  uniform  over  any  particular 
area,  then  one  might  be  able  to  make  meaningful  measurements  of 
microwave  radiance  by  usj.ng  relative  measurements. 


Table  XIII.  Summary  of  calculated  and  observed  values  of  polarized  sky  temperature  (T 
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**  measurement 


Hydrometeors 


General 

Droplets  of  liquid  water  in  clouds  and  rain  absorb,  scatter, 
and  emit  microwave  radiation.  In  some  cases,  the  scattering  of 
microwaves  can  be  neglected  when  compared  to  processes  of  absorp¬ 
tion  and  thermal  emission.  Several  investigators  have  calculated 
the  polarized  sky  temperature  under  various  conditions  of  rain 
and  clouds  (Servey ,  1964;  Edison,  1966).  The  consideration  of 
this  general  problem  as  it  concerns  applied  microwave  radiometry 
has  just  begun.  Recently,  Kreiss  (1968)  calculated  theoretical 
values  of  microwave  radiance  over  marine  clouds  by  assuming  that 
scattering  is  negligible  and  that  ac,  the  absorption  coefficient 
of  cloud  droplets,  is  proportional  to  M,  the  liquid  water  content 
of  the  cloud  (grams  per  cubic  meter) ,  and  the  square  of  the  elec¬ 
tromagnetic  frequency.  Comparison  of  these  theoretically  derived 
values  to  aircraft  measurements  of  microwave  radiance  over  clouds 
in  the  Gulf  of  Mexico  shows  that  scattering  can  be  neglected  in 
all  cases  except  where  the  maximum  drop  sizes  are  large. 

Due  to  the  freezing  of  the  rotational  motion  of  water  in 
ice,  absorption  of  microwaves  by  ice  Is  quite  low.  This  leads  to 
a  correspondingly  low  thermal  emission  for  ice  clouds.  Flights 
of  microwave  radiometers  in  aircraft  have  shown  that  ice  clouds 
are  essentially  transparent  to  microwave  radiation  and  do  not 
significantly  emit  microwave  radiation. 


Mie  theory 


Mie  (1908)  developed  the  rigorous  solution  to  the  problem  of 
predicting  the  absorption  and  scattering  of  electromagnetic 
harmonic  waves  by  a  homogenous  dielectric  sphere  of  radius  a  and 
complex  index  of  refraction  m. 

Since  the  absorption  and  scattering  of  microwaves  by  clouds 
of  liquid  water  are  caused  by  a  finite  number  of  individual 
droplets,  it  is  convenient  to  introduce  the  concept  of  effective 
cross-sections  as  used  in  radar  meteorology.  Suppose  the  is 
the  power  density  (watts  per  square  meter)  of  radiation  incident 
upon  a  particular  droplet  of  liquid  water.  The  rate  of  extinction 
of  energy  by  scattering  (Pg ,  watts)  or  by  absorption  (P^ .watts) 
can  be  given  in  terms  of  an  effective  scattering  cross  section 
(Qg,  square  meters)or  an  effective  absorption  cross  section 
(Q  ,  square  meters),  respectively.  The  effective  cross  section 

A 

may  bear  little  relationship  to  the  actual  cross-sectional  area 
of  the  water  droplet.  The  energy  passing  through  these  effective 
cross  sections  is  the  energy  absorbed  or  scattered,  that  is: 

Qa  =  PA/W][  (m2),  (118) 

and 

(^  =  Pg/Wj.  (m2).  (119) 


The  rate  at  which  energy  is  removed  from  a  pencil  of 
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radiation  (P^,  watts)  can  be  given  in  terms  of  the  total 

2 

extinction  cross-section  (m  X  where: 


Or  -  %  +  qa 


(m2) .  (120) 


The  backscattering  of  microwave  energy  by  water  droplets  is 
of  importance  to  radar  meteorology  and  may  be  important  to  micro- 
wave  radiometry.  According  to  Bat tan  (1959),  Qg,  the  backscat¬ 
tering  cross  section,  is  "the  area  intercepting  that  amount  of 
power,  which,  if  scattered  isotropically,  would  return  to  the 
receiver  an  amount  of  power  equal  to  that  actually  received." 
Backscattering  of  solar  microwave  radiation  by  clouds  and  rains 
may  be  an  important  part  of  the  microwave  radiance  of  clouds  as 
viewed  from  above  during  daytime. 

0^,  Qa,  and  are  functions  of  a,  the  drop  radius,  v,  the 
frequency  of  the  radiation,  and  T,  the  temperature  of  the  drop. 

Aden  (1952)  gives  the  following  expressions  for  Qj»  Qg  aiul 
Q  from  the  Mie  theory: 

D 


*  -(*a2/62)  l  (An  +  2)  Re  (aR  +  bR) 
n=l 


(m2). 


(121) 


Q„  -  (*a2/52)  l  (4n  +  2)(|a  |2  +  |b  |2) 
*  n=l 


(m2). 


(122) 


and 


QB  -  (ira2/62) 


l  <2n  +  l)(an  -  bn) 

n*l 


(n2). 


2 


(123) 


S3 


where: 

<5  =  2ua/X  =  2irav/c  (124) 

and  X  is  the  wavelength  (meters). 

a^  and  are  the  complex  scattering-amplitude  coefficients 

and  are  expressed  in  terms  of  spherical  Bessel  functions  of  the 

first  and  third  kind  of  order  n.  See  Stephens  (1960)  for  exact 

expressions  for  a  and  b  . 

n  n 

The  absorption  or  scattering  per  meter  is  dependent  upon 
n(a),  the  number  of  water  droplets  per  unit  volume  of  radius  a, 
that  is: 

ed(v)  =  n(a)  QT  (a,v,T)da  (nepers),  (125) 

ou(v)  =  I*™  n(a)  Q.  (a,v,T)da  (nepers),  (126) 

a  0  A 

and 

sd(v)  =  /3m  n(a)  Qs  (a,v,T)da  (nepers),  (127) 

where  am  is  the  maximum  droplet  radius  found  in  the  cloud  and 
where  the  subscript  d  denotes  phenomena  associated  with  water 
droplets  only.  The  expressions  (125)  -  (127)  might  be  under¬ 
stood  best  as  the  cross  section  (m^)  per  unit  volume  (m^) . 

The  Mie  expressions  for  Q  ,  Q  and  Q„  are  quite  complicated, 
and  only  a  few  calculations  of  these  quantities  have  been  made  and 
reported  in  the  literature. 

Approximations  to  the  Mie  theory 


Ryde  (1946)  was  able  to  compute  values  of  attenuation  of 
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microwave  energy  by  clouds  and  rain  based  on  accurate  determina¬ 
tions  of  complex  dielectric  constant  of  water  by  Saxton.  He 
assumed  that  the  distribution  of  drop  sizes  was  that  given  by 
Laws  and  Parsons  (1943).  Unfortunately*  he  only  discussed 
and  Qg. 

In  the  case  of  small  water  droplets  (f  0.05),  Ryde  (1946) 
found  that: 

e  =  a  =  0.000832  C. Mv/cp T  (nepers) 

C  C  .  1  L 

=  4.093  (^MVc  l  (db/km),  (128) 

where  M  is  the  mass  concentration  of  liquid  water  (grams  of 
liquid  water  per  cubic  meter  of  air),  is  the  water  density 
(grams  per  cubic  centimeter) , c  is  the  speed  of  light  (centimeter 
per  second),  v  is  the  frequency  (Hz)  and  is  dependent  upon 
«-he  temperature  and  frequency.  Note  that  scattering  is  negli¬ 
gible  in  this  case  and  that  is  proportional  to  M.  Values  of 
ac/M  for  temperatures  of  0°  C  and  20*  C  were  computed  by  Ryde 
and  are  given  in  Table  XIV.  The  subscript  c  refers  to  cloud 
droplets . 

Table  XIV.  The  attenuation  coefficient  of  non-raining  clouds  and 
for  (ac,  db/km)  per  unit  concentration  of  liquid 
water  (M,  g/m3)  for  temperatures  of  0°  C  and  20°  C 
and  for  frequencies  of  3,  6,  10,  15,  and  30  GHz  (after 
Ryde,  1946) . 


Frequency  (GHz) 


3 


6 


10 


15  30 


Goldstein  (1951)  reasoned  that  o-ie  could  neglect  quadrupole 
and  higher  oscillations  in  the  water  droplet  and  consider  only 
the  scattering  and  absorption  that  results  from  dipole  inter¬ 
actions.  Thus,  all  scattering  amplitude  coefficients  except  a^, 
b^,  and  b^  could  be  assumed  to  be  zero.  For  small  drop  sizes 
(6  <  0.1),  one  can  express  the  Bessel  functions  as  infinite  series 
in  terms  of  ascending  powers  of  6.  Truncating  terms  involving 
6>  or  higher,  he  found  that: 

Qt  -  (c263/2ttv2)  (C2  +  C362  +  C463)  (m2),  (129) 

where  c  is  the  speed  of  light  (m/sec)  and  v  is  the  frequency  (Hz). 
The  constants  C^,  Cj  and  are  given  by: 

C2  =  6E'7[(e'  +  2)  +  e"2)  (dimensionless),  (130) 

C3  =  (e'715) 

+  25/ [(2 e'  +  4e"2]  +  l]  (dimensionless),  (131) 


3(7e'2  +  4 e’  -  20  +  It"2) 
[(e'  +  D1  +  e"2] 2 


r  -  1  r(e,-l)2(e,+2)2  +  e'^me'-lHs'-^-gHe"4' 
°4  “  3  [  [ (e+2)2  +  e"2)2  . 


(dimensionless) .  (132) 


If  e”  =  0,  C2  and  are  equal  to  zero, 

c  -  *  (e'  -  l)-2 

4  3  (e*  +  2) 


dimensionless) , 


(133) 
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and 

Qt  =  Qs  -  (26266/3^v2)  | (e ’  -  l)/(e’  +  2)|2  (m2),  (134) 

which  is  the  expression  for  the  scattering  cross-section  according 
to  Battan  (1959) . 

If  6<<1,  then  the  terms  of  (129)  involving  C0  and  may  be 
neglected  and 

Qt  »  Qa  =  4tt2vC2  a3/c  (m2) .  (135) 

Equation  (135)  is  known  as  the  Rayleigh  Approximation  for 
the  absorption  cross-section. 

The  volume  of  a  spherical  droplet  of  radius  a  is  (4ira i  13) 
cubic  meters.  Regardless  of  the  distribution  of  drop  sizes  in 
one  cubic  meter  of  atmosphere,  if  there  are  M  grams  of  liquid 
water  per  cubic  meter,  then  it  follows  that: 

QA  =  4.093  MClV/c  (m2).  (136) 

According  to  Goldstein  (1951),  is  proportional  to  frequency 

and 

a  =  0.000487  Mv2  (db/krn) ,  (137) 

c 

where  v  is  the  frequency  (GHz)  and  M  is  the  liquid  water  content 
(g/m3) . 

Values  of  ac/M  (db/km  per  g/m3)  calculated  by  Goldstein  are 
presented  in  Table  XV  for  18°  C. 
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Table  XV.  Absorption  coefficient  of  non-raining  clouds  and  fog 
(ac,  db/km)  per  unit  concentration  of  liquid  water 
(M,  g/m^)  for  a  temperature  of  18°  C  (after  Goldstein, 
1951). 

O 

Frequency,  v  (GHz)  otc/M  (db/km  per  g/m  ) 


3 

0.0045 

6 

0.0178 

10 

0.05 

15 

0.112 

24 

0.28 

30 

0.*38 

42.9 

0.876 

60 

1.65 

150 

7.14 

The  liquid  water  content  of  fog  may  be  as  high  as  one  gram 

per  cubic  meter.  Comparison  of  values  cf  in  Table  XV  to  values 

of  a  the  absorption  coefficient  of  molecular  oxygen  and  water 
8 

vapor  given  in  Table  X,  shows  that  absorption  and  emission  by 
clouds  and  fog  are  higher  than  absorption  and  emission  of  atmo¬ 
spheric  gases  for  frequencies  of  10  GHz  or  higher  and  for  surface 
pressure  and  a  temperature  of  18°  C. 

Goldstein  (1952)  gives  a  set  of  correction  factors  for  ac 
for  cases  where  the  temperature  is  not  equal  to  18°  C.  Let  F. 

JL 

(T)  be  such  that: 

ct  (T)  =  a  (18°C)  F.  (T) 
c  c  l 

Values  of  F^  (T)  are  given  in  Table  XVI. 


(db/km) ,  (138) 


oo 


ww 


Table  XVI.  Correction  factor,  F..(T)  (dimensionless),  for  a 
in  Table  XV  (after  Goldstein,  1952).  C 


Frequency,  v(GHz) 

8°  C 

o 

0 

O 

00 

O 

O 

20°  C 

30°  C 

-O 

o 

o 

3 

2.0 

1.25 

1.0 

0.95 

0.67 

0.59 

9.4 

1.98 

1.30 

1.0 

0.95 

0.7 

0.56 

24 

1.93 

1.29 

1.0 

0.95 

0.73 

0.57 

60 

1.59 

1.2 

1.0 

0.95 

0.73 

0.59 

Gcldstein  discusses  the  attenuation  bv  raindrops  and  calcu¬ 
lates  values  of  e^,  the  extinction  coefficient  for  raindrops,  for 
various  rates  of  precipitation  (R,  millimeter  per  hour)  based  on 
the  distribution  of  d  "ip  sizes  given  by  Laws  and  Parsons.  Some 
of  his  values  ?~e  listed  in  Table  XVII. 

From  Table  XVII,  one  can  see  that  very  substantial  attenu¬ 
ation  takes  place  for  heavy  rains  and  for  high  frequencies. 
Scattering  is  the  dominant  extinction  process  for  raindrops  having 
large  diameters  and  for  high  frequencies  (Goldstein,  1951) . 

Gunn  and  East  (1954)  have  made  calculations  of  eR  based  on 
distributions  of  drop  sizes  other  than  that  given  by  Laws  and 
Parsons  (1943).  Again,  no  distinction  was  made  between  scattering 
and  true  absorption.  Other  investigators  who  have  calculated 

Q_  oi  e„  using  the  Mie  theory  are  Tmai  (1957)  and  Stephens  (1961). 
T  R  1  " 

Emission  of  clouds  and  rain 


According  to  Deirmendj ian  (1963) , 
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a£  =  1.11  x  10  ^  (nepers),  (139) 

where  v  Is  the  frequency  (GHz)  and  M  is  the  liquid  water  content 
of  the  cloud  (g/m  ) .  Equation  (139)  is  valid  only  for  non¬ 
precipitating  clouds.  This  relationship  is  similar  to  (137). 

Following  a  suggestion  by  Buettner  (1963,  1965),  Kreiss 
(1968)  used  (139)  to  compute  values  of  polarized  sky  temperature 
under  various  conditions  of  cloud  thicknesses,  cloud  heights 
and  concentrations  of  liquid  water.  Kreiss  considered  forty- 
two  cases  which  are  listed  in  Table  XVIII.  See  his  paper  for 
mor-.  aetails.  Space  does  not  permit  the  presentation  of  his 
complete  tabulations,  but  some  values  are  presented  in  Table 

XIX.  His  curves  give  TS^  versus  v  for  (15  GHz  <  v  <  30  GHz). 

P  ~  ” 

The  values  of  polarized  sky  temperature  in  Table  XIX  show 
how  drastically  TS^  is  affected  by  various  meteorological 
conditions.  Since  the  sea  reflects  approximately  two-thirds 
of  the  microwave  radiation  incident  upon  its  surface  at  the  nadir 
angle,  it  is  easy  to  see  that  much  of  the  upvelling  radiation 
at  the  frequencies  shown  in  Table  XIX  originates  from  thermal 
emission  by  atmospheric  constituents.  Also,  tall  cumulus  clouds 
act  nearly  as  black  bodies  for  even  moderate  values  of  M.  One 
should  notice  that  the  polarized  sky  temperature  is  not  linearly 
proportional  to  M, the  liquid  water  content,  although  a£  is  pro¬ 
portional  to  M. 

Decker  and  Dutton  (1968)  were  able  to  map  the  liquid  water 


r 


91 


m 

#» 

in 

u 

#k 

in  * 

CM 

»  CM 

«> 

CM 

*» 

d) 

O 

m 

CM  CM 

fH 

m  »H 

m 

rH 

in 

u 

• 

• 

• 

• 

•  • 

• 

• 

• 

cfl 

, — . 

o 

O  * 

o 

0-4-0 

o 

>3-  O 

o 

►> 

rH 

1 

•u 

•» 

#> 

r. 

r»  r>  *k 

• 

•k  #k 

r* 

T3 

C 

m 

in 

m  » 

in 

m  cm  «n 

in 

cm  in 

m 

<r 

*H 

(Den 

CM 

cm 

CM  1/5 

CM 

CM  CM 

CM 

CM 

CM 

3 

n  E 

vO 

MO 

sO  • 

vr> 

•  #»  \£> 

• 

-  sO 

• 

#» 

c r 

C  \ 

O 

O  rH 

o  o  <r 

o 

O  rH  O 

o 

tH  O 

o 

rH 

*w 

o  oe 

-K 

* 

• 

• 

• 

• 

• 

• 

►j 

O  MX 

ZZZZZZZZZZ 

o 

o 

o 

o 

o 

o 

CO 

cn 

a 

c. 

X 

u 


03 

o 

rH 

u 


X  •  . 
H  O 


E 


:z:  z  z:  :z:  z  ^ 


O'  iP  o  ^ 

cm  to  co  cm 
2  2  X  O  M  C 


rH 

<r 

r- 

CM 

00 

m 

CO 

CO 

m 

iH 

CM 

s£> 

05 

-H 

O 

U 

X 

>> 

-C 

T3 

O 

U 

3) 

T3 


0) 

CZ 

U 

u-» 

O 


05 

t-t 

►J 


> 

X 


tfl 

H 


O  T3 

3 

CL  O 
°  -*  _ 
H  rJ  £ 


in 

o 


U-4 

o 

a)  r 
o.  o 
— « 
H  o 


o 

u 

3 

i-» 

C3 

L* 

o 

CL 

£ 

o 


05 

u 

o 

O  X} 
0)  B 
CT3  3 
U  S 


S2 


2:  2:  2:  2 


r^  m 

CO  O  O' 

<r  r-  CC  N 

2  2  r-  ,v  M  O 


r-* 

m 


4:  * 

zz  zz  tz  zz  zz  zz  ^  zz  ?: 


OJ  OJ  OJ  3  2  z.  U 


CO  CO  CO  CO 


O  H  N  c 

HMro<riOvONXC'«^^^  — 


CM 

CM 

I 

cm 


c 

cn 


H 

co 


H 

C/5 

H 


I 

m 

CM 


vC 

r> 


c 

co 


m 

oo 


so 

Os 

a  3 

05  O 

CC 

00 

vj 

Cx 

CC 

00 

CC 

00 

fH 

co  rH 

CM 

CM 

f- , 

U- 

CM 

CM 

CM 

CM 

NX 

os  u  £ 

2  Z  2  2  2  2  Z  2  2 

vO 

C/5 

00 

X 

sO 

vO 

sO 

O 

H 


05 

•H 

3 

CJ 

0 

0 

U 

0 

cj 

0 

c 

♦H 

»— < 

• — « 

tf 

t-* 

h-* 

M 

c 

•H 

CO 

c/5 

C/5 

»r 

CO 

CO 

CO 

CO 

o 

E 

3 

' — ' 

NX 

vx 

NX 

>~x 

'-X 

^x 

V— X 

05 

Q> 

3: 

O  05  2  C  X  2  C  x  D:  x 

C/5 

CC 

C/5 

CO 

CO 

CO 

CO 

CM 

<r 

i 


92 


93 


in 

<n 

H 

CM 

3 

• 

m 

3 

CM 

O 

o 

<r 

<r 

cn 

*H 

CM 

H 

m 

rH 

vO 

r- 

00 

o 

rH 

O 

/-N 

>4 

CM 

o 

z:  in 

<u 

O  co 

4J 

V-/  • 

3 

CN 

cn 

CM 

ON 

r>. 

O' 

£ 

H 

<r 

CM 

CM 

CJ 

>N 

C/5 

o 

3 

c 

O 

05 

-H 

3 

Vi 

croc 

3 

o  • 

> 

u  in 

j  cc 

O 

CM 

cn 

[L  rH 

rH 

H 

«— • 

T— 1 

Vl 

05 

*3 

c 

3 

/— \ 

X 

£ 

£ 

o 

CO 

* 

to 

Ln 

m 

0) 

<r 

m 

^  CL 

rH 

cc 

o 

0) 

•3 

1 

1 

3 

cc 

cc 

N—r 

O 

CM 

CM 

*H 

vC 

vC 

3 

1 

a 

Vi 

| 

E 

«= 

3 

* 

0 

5 

4J 

u 

i- 

CO 

4-1 

».  - 

14- 

Vi 

•H 

3 

•3 

•3 

CL 

•H 

3 

3 

E 

E 

4-1 

4-J 

0)  • 

J 

3 

cj 

4-1 

4  ■ 

Vi 

00 

0) 

XnO 

CO 

"3 

4-1 

4J 

^  CN 

a 

rj 

3 

3 

CO  »~t 

o 

i5 

to 

X 

*3  - 

a 

CO 

- 

r 

0)  01 

!  o 

CO 

V4 

CO 

X 

X 

N  05 

CO 

o 

CO 

X 

•H  *H 

rj 

a 

4-1 

*— t 

a 

3 

Vi  C 

o 

rH 

3 

*3 

. — • 

r— 

CO  V4 

*H 

*3 

J- 

”3 

*3 

*— 1  ^ 

4J 

-4 

3 

c 

-4 

o 

c- 

0 

r> 

o 

C 

a-  E 

— 1 

r4 

E 

3 

«— 1 

— • 

o 

Vi 

u 

a 

3 

3 

<44  i4 

c 

4J 

r 

C  <m 

to 

♦ 

> 

•> 

r 

' — ' 

0/ 

4-1 

"3 

u 

*3 

CO 

o 

*3 

U 

L 

3  co 

l> 

3 

3 

3 

3  C 

VH 

r 

*3 

*3 

•-«  O 

3 

;  •* 

C 

c 

3 

3  *H 

*H 

4-1 

3 

r— 

7Z 

3 

>  -W 

u 

J  _£ 

44 

3 

4_> 

«—> 

■*  1 

CO 

CO 

CO 

CO 

0*  *3 
£  C 

o  o 

co  3 


m 


e'¬ 

en 


it 


u 


O' 

CM 


cc 

c 


3 


co 


m 

cc 


oo 

<r 


v£> 

cm 


cn 

£ 


II 


O' 

cm 


u 

3 


CO 


r- 

r^ 


O 

<r 


t's 

O 

II 


£ 

r- 

in 

<r 

l 

/v\ 

CM 

vC 


C 

1~ 


X 

3 

4-> 

u 


C/I 


3 


4-» 

CO 


o 

ON 


o 

NO 


U) 

in 

CM 


in  O 

co 

CM  CM 


CM 

m 


00 

o 

m 

O 

«H 

cn 

vO 

CO 

t-H 

CM 

CM 

CM 

Vi 

• 

a 

3 

X 

.£5 

rH 

►H 

g 

X5 

X 

3 

E 

C-I 

o 

*— 1 

H 

3 

E 

*H 

r% 

X 

O 

3 

3 

Vi 

> 

H 

U 

L- 

X 

ON  o  cm  O  f-H  CM 

— •  cm  cm  <r  <r  <r 


94 


content  of  thunderstorms  over  Colorado  using  a  microwave  radiometer 
operating  at  10.7  GHz.  They  assumed  that  scattering  was  negligible 
and  thus  computed  values  of  M  times  column  length  that  are  unreal¬ 
istically  high. 

Cummings  and  Hull  (1966)  measured  sky  radiance  at  23  GHz 
and  35  GHz  under  various  meteorological  conditions.  These  values 
are  shown  in  Table  XX  using  their  word  descriptions  of  the  meteoro¬ 
logical  conditions  existing  at  the  time  of  the  observations. 


Table  XX.  Average  values  of^neasured  vertically  polarized  sky 
temperatures  (Tvs  *  ,  °K)  at  Columbus,  Ohio,  for  fre¬ 
quencies  of  23  GHz  and  35  GHz  taken  in  the  summers  of 
1964  and  1965  (after  Cummings  and  Hull,  1966). 

23  GHz  35  GHz 


Meteorological  condition 

**  =  0° 

*  =  45° 

It 

o 

o 

<!>  =  45° 

Clear  sky 

30 

47 

30 

40 

Fog 

— 

— 

40 

42 

Light  cloud  cover 

40 

60 

35 

50 

Moderate  dark  cloud  cover 

45 

100 

40 

65 

Dark  cloud  cover 

140 

160 

50 

70 

Light  rain 

200 

200 

90 

90 

Heavy  rain 

255 

255 

120 

120 

*<{/  is  the  zenith  angle. 


Comparing  measured  values  in  Table  XX  to  theoretical  values 
in  Table  XIX,  one  sees  that  the  atmosphere  was  slightly  cooler 
and  dryer  than  the  1962  Standard  Atmosphere  assumed  by  Kreiss. 

The  values  of  measured  sky  temperature  agree  very  well  with  cal¬ 
culated  values. 

Weger  (1960)  also  has  calculated  values  of  sky  temperature 
for  various  conditions  of  cloud  cover  and  rain  using  a  simplified 
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engineering  approach.  Table  XXI  shows  the  results  of  his  calcula¬ 
tions.  He  assumed  that  the  ARDC  Model  Atmosphere  is  valid  and 
that  relations  given  in  (IjlO)  and  (111)  are  valid. 

Table  XXI.  Values  of  vertically  polarized  sky  temperature 

(V*y,  °K)  by  Weger  (1960)  for  various  meteorological 
conditions . 


Meteorological  Conditions 

10 

Frequency, 

16.7 

.)  (GHz) 
24 

34.9 

Clear  sky 

4 

8 

27 

19 

Un^orm  cloud,  900-1800  m 

M  =  0.3  g/m^ 

4 

9 

32 

29 

Uniform  rain,  0-900  m 

R  =  4  mm/hr 

8 

20 

57 

77 

Weger's  values  differ  greatly  from  Kreiss'  values  in  Table 
XIX.  This  is  probably  a  reflection  of  the  simplified  approach 
used  bv  Weger. 

Wulfsberg  (1964)  measured  sky  temperatures  under  various 
conditions  of  cloudiness  and  rain.  His  results  are  given  in 
Table  XXII. 

Wulfsberg' s  measured  values  agree  closely  with  Weger's 
theoretical  values  at  near  17  GHz,  but,  they  differ  at  35  GHz. 

Summary 


From  the  preceeding  tables,  one  can  see  that  the  presence 
of  liquid  water  in  the  atmosphere  can  greatly  affect  the  intensity 
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Table  XXII.  Measured  values  of  polarized  zenith  sky  temperature  by 
Wulfsberg  (1964)  for  various  meteorological  conditions. 


Meteorological 

Condition 

15 

Frequency, v  fGHz) 
17 

35 

Cold,  dry,  clear  atmosphere 

3 

4 

10 

Hot,  wet,  clear  atmosphere 

8 

12 

29 

Overcast  clouds 

13 

— 

57 

Light  rain  (1  mm/hr) 

14 

20 

41 

Moderate  rain  (4  mm/hr) 

31 

26 

112 

of  downward  radiation  from  the  atmosphere.  All  of  the  quantitative 
calculations  have  been  made  with  the  assumption  that  clouds  do 
not  scatter  significantly  and  that  the  atmosphere  is  horizontally 
uniform  in  a  state  of  local  thermodynamic  equilibrium.  The  extinc¬ 
tion  in  rain  is  largely  due  to  scattering  (Deirmend j ian ,  1963) 
for  centimeter  wavelengths  (3-30  GHz).  One  is  prompted  to  ask  the 
question:  How  does  scattering  affect  the  relationship  between  the 
environment  and  the  intensity  of  microwave  radiation? 

Solar  microwave  radiation  is  a  major  source  of  extra-ter¬ 
restrial  microwave  radiation.  Forward  scattering  and  back- 
scattering  of  solar  radiation  could  result  in  a  significant 
addition  to  the  thermally  emitted  sky  radiation.  Side  scattering 
of  microwave  radiation  impinging  upon  a  cloud  from  horizontal 
directions  is  another  possible  source  of  emission  by  scattering. 

On  the  other  band,  one  must  consider  the  greatly  increased 
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extinction  that  results  in  clouds  due  to  scattering.  If  multiple 
scattering  is  dominant,  cne  would  expect  the  scattered  radiation, 
to  be  isotropic  and  equal  to  the  average  flux  per  unit  solid 
angle  of  the  radiation  intercepted  by  the  cloud.  Before  this  matter 
can  be  pursued  further,  one  must  determine  the  polarized  intensity 
of  microwave  radiation  emerging  from  the  surface.  This  is  the 
subject  of  Chapter  V. 


i 
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CHAPTER  V 

EXTINCTION  AND  EMISSION  IN  THE  HYDROSPHERE 

General 

When  microwave  radiation  encounters  the  air-sea  interface, 
it  is  partially  reflected  or  scattered  back  into  the  atmosphere. 
The  remaining  radiation  enters  the  hydrosnhere  and  is  quickly 
absorbed  since  sea  water  is  a  very  absorptive  material.  The 
hydrosphere  also  emits  microwave  radiation  into  the  atmosphere. 

One  may  use  certain  models  of  th^  hydrosphere  to  predict  its 
emissive  and  reflective  characteristics.  As  a  first  approxi¬ 
mation,  one  may  assume  that  the  air-sea  interface  is  perfectly 
flat  and  that  the  active  surface  layer  is  homogenous  in  its 
temperature  and  dielectric  properties.  The  active  surface  layer 
of  the  hydrosphere  is  that  part  of  the  hydrosphere  in  which  the 
microwave  energy  fiom  the  atmosphere  is  absorbed  and  the  ther¬ 
mally  emitted  microwave  radiation  originates.  The  active  surface 
layer  is  usually  no  thicker  than  a  few  multiples  of  the  wave¬ 
length  of  the  microwave  radiation  under  consideration.  In  the 
case  of  two  homogenous  dielectric  materials  separated  by  a  flat 
boundary,  one  may  use  Fresnel's  Laws  of  Reflection  to  predict 
the  polarized  reflectivities  of  the  interface.  Through  thermo¬ 
dynamic  considerations,  one  may  relate  the  polarized  emissivities 
of  the  hydrosphere  to  its  polarized  reflectivities  by  use  of 
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(58)  and  (59). 

The  actual  air-sea  interface  is  not  flat;  indeed,  it  is 
quite  rough.  In  fact,  it  is  sometimes  impossible  to  define  a 
single  air-sea  interface  under  conditions  of  high  shearing  stress 
when  bubbles  and  foam  form  on  the  sea  surface.  Blown  salt  spray 
may  also  complicate  the  problem. 

One  may  add  a  certain  degree  of  reality  to  the  case  of  a 
perfectly  flat  sea  surface  by  assuming  that  the  air-sea  interface 
has  an  irregular  shape  which  is  smooth.  According  to  Beckman 
and  Spizzichino  (1963) ,  a  surface  may  be  treated  as  a  smooth 
surface  if  its  radius  of  curvature  is  many  times  larger  than  the 
wavelength  of  the  radiation  under  consideration.  These  assump¬ 
tions  hold  very  well  for  the  actual  sea  surface  when  infrared 
and  visible  radiation  are  under  consideration.  In  the  case  of 
microwave  radiation,  only  gravity  waves  have  radii  of  curvature 
suitable  to  the  condition  of  smoothness  required.  At  any  rate, 
one  may  assume  that  Fresnel's  law  of  reflection  is  valid  at 
a  particular  point  on  the  surface,  a..  ..e  may  compute  mean 

polarized  emissivities  by  integrating  over  the  slope  distribution 
(Buettner  et  al.,  1968). 

Finally  one  may  treat  the  sea  surface  as  a  scattering  sur¬ 
face.  A  knowledge  of  y,  the  bistatic  scattering  cross  section 
per  unit  area,  is  needed  in  this  case.  Scientists  have  been 
unable  to  determine  the  form  of  y  from  purely  theoretical  grounds. 
Thus,  they  resort  to  making  educated  guesses  as  to  the  form  of  y 
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for  the  sea  and  other  natural  surfaces. 

Peake  et  al. ,  (1966)  have  suggested  several  forms  for  y 
for  land  surfaces.  Stogryn  (1967)  and  Janza  (1968)  have  pro¬ 
posed  models  for  y  for  the  sea  surface. 

The  purpose  of  this  chapter  is  to  review  the  problem  of 
predicting  the  polarized  reflective  and  emissive  temperatures  of 
the  sea  for  the  three  cases  of  a  flat  surface,  an  irregular 
smooth  surface,  and  a  rough  surface  of  the  sea. 

Flat  Sea  Surface 


Fresnel's  Laws  of  Reflection 


The  polarized  electric  reflection  coefficients  for  two 
homogenous  dielectrics  separated  by  a  flat  boundary  are  given 
by  Jordon  (1950)  as: 

2 

Pv  =  (e2/ei)cos  ij)-  /(e2/d)  -  sin  (dimensionless) ,  (140) 

(E2^Gl)cos  ~  sin'^ 

and 

o 

cosiii  ~  /( to/ci)  ~  sin~i><  (dimensionless) ,  (141) 

COSlJj  +  A£2/cX)  -  sin^ 

where: 

p  =  vertically  polarized  reflection  coefficients  for  the 
electric  field  (dimensionless), 

=  horizontally  polarized  reflection  coefficients  for 
the  electric  field,  (dimensionless). 
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e^  =  the  complex  dielectric  constant  of  the  first  medium 
(dimensionless) , 

e 2  =  the  complex  dielectric  constant  of  the  second  medium 
(dimensionless) , 

and 

ifc  =  angle  of  incidence. 


The  radiation  is  assumed  to  be  plane  polarized  radiation 
coming  from  medium  'one'  toward  medium  'two*.  In  the  case  of 
a  homogenous  layer  of  sea  water  overlain  by  air,  e^,  is  equal 
to  unity,  and  is  equal  to  e’-  je"  as  given  in  Chapter  II  for 
sea  water. 

The  polarized  reflectivity  is  equal  to  tne  square  of  the 
polarized  electric  reflectivity.  Thus: 


r 


v 


(e*  -  je")cos-J>  -  Ac*  -  je")  -  singly 
(e'  -  je")cosi!>  +  /fe'  -  je")  -  sin^P 


2 


(dimensionless) , 


(142) 


and 


rh  = 


cosil1  -  /(e*  -  je")  -  sin^ft 


cos<|<  +  v/(e'  -  je")  -  sin^O 


(dimensionless) . 


(143) 


Reflective  temperature  of  a  flat  sea  suxrace 

Values  of  e'  and  e"  were  computed  in  Chapter  II  and  were 
given  for  specific  values  of  temperature,  ralinity,  and  frequency 
in  Table  VII  as  obtained  from  equations  (33),  (34),  (36),  (37), 
and  (38).  Using  these  equations,  one  may  compute  r  as  a  function 


102 


of  v,  S,  and  T  for  any  given  and  polarization.  If  TS' ^  is  also 

P 

known,  then: 

Tj  (v,S,T,*>  =  rp  (v,S,T,*)  T*£y  (v,*)  (°K)  (1AA) 

Emissive  temperature  of  a  flat  iea  surface 

When  a  flat  sea  surface  is  under  conditions  of  local  thermo¬ 
dynamic  equilibrium,  the  polar „ed  emissivities  of  that  surface 
are  given  as: 


£  -  1  -  r 

P  P 


(dimensionless) .  (1A5) 


From  (58)  and  (1A5)  it  follows  that: 


T®  «  (1  -  r  )  T 
p  p  sea 


(°K) , 


(1A6) 


where  Tp  is  the  polarized  emissive  temperature  of  the  sea. 

Brightness  temperature  of  a  flat  sea  surface 

The  term  polarized  brightness  temperature  (Tp,  °K)  is  now 
defined  as  the  effective  blackbody  temperature  associated  with 
the  total  polar xzed  upwelling  microwave  radiaticnal  intensity. 
Tp  is  a  function  of  height  and  other  variables.  Now,  it  may  be 
stated  that: 


T  (z  =  0)  =  r  Tsky  +  (1  -  r  )T 
p  P  P  P  sea 


(°K). 


For  purposes  of  illustration,  values  of  T  (0)  have  been 

P 


(1A7) 


i  no 


computed  for  an  incidence  angle  of  0°.  Table  XXIII  shows 
e  r 

values  of  r  ,  ?  ,  T  ,  T  ,  and  T  for  frequencies  of  1,  6, 

P  P  P  P  P 

9  3,  15.8,  19.35,  22.235,  ?:;d  34  GHz,  for  temperatures  of  10°  C 
and  20°  C,  and  for  salinities  of  10  and  35  o/oo.  The  assumed 
values  of  Tsj^  were  taken  from  Table  XXITI  for  average  meteoro¬ 
logical  conditions. 

Measurements  of  the  microwave  radiance  of  the  sea  are 
available  for  incidence  angles  near  0°  and  for  frequencies  of 
9.3,  15.8,  19.35,  22.235,  and  34  GHz.  From  Table  XXIII,  it  is 
seen  that  does  not  change  more  nan  3°  K  for  a  change  in  sea 
temperature  from  10°  C  to  20°  C  at  these  frequencies.  This  is 
due,  in  part,  to  the  fact  that  the  polarized  emissivity  of  the 
sea  decreases  with  increasing  temperature;  consequently,  the 
polarized  emissive  temperatures  of  the  sea  remain  fairly  constant. 
In  fact,  for  a  frequency  of  34  GHz,  the  sea  emits  less  radiation 
when  the  temperature  is  20J  C  than  when  it  is  10°  C! 

The  frequency  domain  below  8  GHz  is  of  interest  because  of 
the  low  attenuation  by  the  atmosphere  and  the  small  values  of 
polarized  sky  temperature.  The  values  in  Table  XXIII  for  fre¬ 
quencies  of  1  GHz  and  6  GHz  show  many  interesting  features 
about  t’.e  emission  and  reflection  of  the  sea  surface  for  these 
low  microwave  frequencies.  The  effect  of  changes  in  salinity 
on  Tp  at  1  GHz  is  quite  large.  For  a  water  temperature  of  20°  C, 
a  change  in  salinity  from  10  o/oo  to  35  o/oo  results  in  a  decrease 


of  Tp  by  18.9°  K.  There  is  also  a  pronounced  effect  on  due  to 
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temp  . mature  at  this  frequency. 

0 

Figure  8  shows  a  plot  of  T  versus  T  for  frequencies  of 

p  sea 

1,  5.4,  9.3,  15.8,  19.35,  22.235,  and  34  GHz  and  for  various 
values  of  salinity.  The  effect  of  salinity  is  negligible  for 
frequencies  higher  than  5.4  GHz.  The  curves  in  Figure  8  show 
that  the  same  emissive  temperature  may  be  the  result  of  two 
different  values  of  sea  temperature  for  frequencies  of  15.8, 

19.35  and  22.234  GHz.  At  5.4  GHz,  there  is  no  salinity  effect, 
and  it  appears  that  the  emissivity  is  practically  constant  with 
temperature.  The  pronounced  effect  of  salinity  at  1  GHz  is 
obvious  from  Figure  8.  Also,  it  should  be  recalled  from  Chapter 
IV  that  the  polarized  sky  temperature  for  these  frequencies  is 
small  and  fairly  constant  since  it  is  mostly  the  result  of  pressure 
broadening  oxygen  emission.  The  natural  variability  of  water 
vapor  has  little  effect  at  these  frequencies. 

It  is  proposed,  therefore,  that  the  temperature  and  salinity 
of  sea  water  could  be  measured  remotely  bv  use  of  a  dual-channel 
microwave  radiometer  operating  near  5.4  GHz  and  1  GHz.  The  first 
step  would  be  to  determine  the  temperature  of  the  water  through 
measurements  at  5.4  GHz,  and,  having  determined  the  temperature, 
one  could  infer  the  salinity  of  the  water  from  the  measurements 
at  1  GHz.  The  accuracy  of  these  environmental  measurements  is 
of  course,  dependent  upon  the  accuracy  of  the  remote  measurements. 
Also,  the  existence  of  foam  on  the  sea  surface  would  invalidate 


these  measurements. 
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Any  further  discussion  of  these  possibilities  will  be  post¬ 
poned  until  roughness  effects  and  instrument  effects  have  been 
considered . 

Marandino  (1967)  studied  the  brightness  temperature  of  the 

sea  as  a  function  of  sea  temperature  and  frequency.  For  ocean 

0 

water,  he  determined  AT^,  the  change  in  polarized  emissive  tem¬ 
perature  (°K),  that  results  when  the  temperature  of  the  sea  sur¬ 
face  is  changed  from  273°  K  to  300°  K.  His  results  are  shown  in 
Figure  9. 

He  found  that  the  polarized  emissive  temperature  is  most 
sensitive  to  changes  in  temperature  at  5  to  6  GHz  and  is  least 
sensitive  to  temperature  at  1.4  GHz. 

Sirounian  (1968)  studied  the  effect  of  salinity  on  at 
a  frequency  of  3  GHz.  He  found  that  the  effect  of  salinity  was 
largest  for  incidence  angles  near  zero  and  that  a  salinity 
change  from  6.7  o/oo  to  32  o/oo  resulted  in  a  decrease  of  ? 
of  3  to  4°  K  for  a  sea  temperature  of  25°  C  and  for  an  angle 
of  incidence  of  zero.  These  results  agree  with  the  values  ob¬ 
tained  by  Paris  (1968). 

Smooth  Irregular  Sea-Surface 

Buettner,  et.  al.  (1968)  have  considered  the  effect  of 
roughness  on  Tp  for  centimeter  wavelengths.  They  assumed  that 
the  sea  surface  acts  as  a  specular  reflector  for  any  given  facet 
of  the  surface  and  determined  the  effective  emissivitv  by 
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integrating  over  the  distribution  of  sea  slopes  given  by  Cox 
and  Munk  (1954a,  1954b).  Their  estimates  show  that  the  emissivjrv 
increases  with  increasing  roughness;  but,  the  increase  is  largest 
at  large  incidence  angles  of  observation.  In  fact,  the  emissivitv 
does  not  appear  to  be  a  function  of  roughness  for  incidence  angles 
near  zero.  This  is  probably  due  to  the  fact  that  the  reflectivity 
of  a  flat  sea  surface  is  fairly  constant  with  angle  for  incidence 
angles  near  zero  (Saxton  and  Lane ,  1952) 

In  addition  to  the  change  in  the  effective  emissivitv  of 
the  sea  surface  with  roughness,  one  should  consider  the  effect 
of  roughness  on  the  reflected  radiation.  T\«>  effects  are  present 
in  this  regard.  First,  the  effective  reflectivity  or  albedo 
changes  with  the  distribution  of  slopes  to  a  slight  extent. 

Second,  the  source  of  the  reflected  radiation  changes  since  sky 
radiation  at  all  incidence  angles  is  reflected,  to  some  extent, 
in  a  given  direction  by  the  irregular  surface.  The  sky  temperature 
is  highly  dependent  upon  incidence  angle  for  any  given  atmospheric 
condition  since  the  path  length  increases  with  increasing  inci¬ 
dence  angle.  From  Table  XIII,  one  can  see  that  the  sky  tempera¬ 
ture  may  be  60°  K  higher  for  incidence  angles  of  60°  than 
for  incidence  angles  near  zero  for  frequencies  near  the  absorp¬ 
tion  line  of  water  vapor.  This  effect  is  not  so  severe  for  fre¬ 
quencies  less  than  10  GHz.  In  addition  to  other  effects,  one 
effect  of  roughness  is  to  increase  the  amount  of  reflected 
radiation  for  frequencies  near  22.235  GHz. 


Rough  Sea  Surface 


General 


The  actual  sea  surface  is  rough  and  scatters  microwave 
energy.  To  determine  the  value  of  the  polarized  brightness 
temperature  of  a  rough  surface,  one  must  know  its  bistatic  scat¬ 
tering  coefficients  per  uni.-  area,  as  defined  by  (75)  through  (78) 
skv 

and  T  7  ,  the  polarized  sky  temperatures  as  a  function  of  inci¬ 
dence  angle.  Peake  (1959)  shows  that: 

Th  *  Yhh  +  if?  Yhv]  d!>  CK),  (148) 

and 

Tv  -  <»4”){[Tvky  +  Thky  \„]  d!>  0*).  C1M) 

where  the  integrations  indicated  in  (148)  and  (149)  are  performed 
over  the  upper  hemisphere.  From  (58),  (80),  (81),  (86),  the 
polarized  emissive  temperatures,  and  T^,  are  expressed  as: 

Th  '  Tsea  U  '  f4hl,  +  <*»•  <150> 

and 

Te  =  T  {l  -  (1/4*)]  [y  +  y  , ]dQ}  (°K) ,  (151) 

v  sea  J  w  vh  ’ 

where  the  integration  is  performed  over  the  upper  hemisphere. 

The  definitions  of  bistatic  cross  section  per  unit  area 


U3 


give  no  indication  cf  the  physics  of  scattering  as  such.  The 
problem  of  scattering  of  electromagnetic  waves  by  an  irregular 
surface  cannot  be  approached  in  a  rigorous  fashion.  Some  approx¬ 
imations  have  been  developed. 

Kfrchhoff  approximation 

Stogryn  (1967)  reasored  that  the  Kirchhoff  ->roximation  can 
be  used  to  give  the  scattering  coefficients  per  unt  area  for  a 
very  rough  surface.  He  assumed  that  the  distribution  of  wave 
slopes  given  by  Cox  and  Hunk  (1954a,  1954b)  is  valid  and  that  the 
root-mean-squared  wave  slopes  (g^  and  g^)  depend  on  the  wind 
speed  (w,  m/sec)  at  41  feet  above  the  mean  sea-surface,  tha-  is: 

g^  =  0.003  +  1.92  x  10"3  v,  (152) 

and 

gj  «  3.16  x  10~3  w,  (153) 

where  the  x-axis  lies  cross  wind  and  the  y-axis  lies  parallel 
to  the  wind. 

The  exact  expressions  for  y^,  ''Vv’  ^hv*  anc*  rvh  are  ^u*te 
complex.  The  reader  is  referred  to  the  compact  expression  of 
these  quantities  in  Stogryn's  p.-.per  since  it  ±s  impossible  to 
obtain  any  analytic  information  from  these  expressions  explicitly. 
The  results  of  his  studies  agree,  in  general,  with  the  results  of 
the  previous  case  of  a  smooth  irregular  surface.  In  general. 
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the  horizontally  polarized  brightness  temperature  increases  with 
increasing  wind  speed  for  incidence  angles  between  zero  and  70° 
at  all  frequencies  in  the  microwave  domain.  This  increase  is 
largest  for  incidence  angles  near  50° . 

Qualitatively,  it  appears  that  the  horizontally  polarized 
brightness  temperature  may  be  increased  by  as  much  as  20°  K  due 
to  increased  roughness.  The  vertically  polarized  brightness  temp¬ 
erature  of  the  sea  is  not  noticeably  affected  by  roughness 
(Stogrvn,  1967). 

Any  further  discussion  of  the  roughness  effect  will  be  post¬ 
poned  until  Chapter  VI  where  the  extinction  and  emission  of  the 
intervening  atmosphere  between  the  aircraft  or  spacecraft  and 
the  sea  will  be  considered. 

Foam,  Bubbles,  and  Ice 


General 

In  each  of  the  three  cases  treated  previously,  it  was  assumed 
that  a  single  liquid-gas  interface  existed  between  these  two  homo¬ 
genous  dielectrics.  This  surface  was  allowed  to  be  quite  irregular 
and  rough.  In  the  case  of  the  real  air-sea  boundary,  large  surface 
shearing  stresses  often  result  in  the  formation  of  bubbles  and 
foam  on  the  sea  surface.  In  this  case,  one  cannot  use  equations 
developed  for  a  single  boundary  between  two  homogenous  dielectrics. 
There  is  no  physical  theory  that  is  adequate  to  describe  the 
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effective  dielectric  constant  of  this  domain  consisting  of  multiple 
layers  of  water  films  and  air  spaces  of  random  orientation  and 
thickness. 

Also,  in  some  parts  of  the  World  Ocean,  the  surface  of  the 
hydrosphere  is  covered  to  some  extent  by  ice  of  various  thick¬ 
ness  and  shape.  The  great  difference  between  the  dielectric 
constant  of  ice  and  of  sea  water  has  been  noted  in  Chapter  II. 

This  difference  should  result  in  a  large  difference  between  the 
polarized  brightness  temperature  of  sea  water  and  that  of  sea  ice 
on  the  water. 

Foam  and  bubbles 

The  only  information  in  the  literature  on  the  effects  of 
foam  on  the  polarized  brightness  temperature  of  the  sea  is  the 
experiment  performed  by  Wi 11 lams  (1968).  Using  the  microwave 
radiometers  on  the  Convair  240A  operated  by  NASA,  he  made 
measurements  of  the  antenra  temperature^  of  water  in  a  large 
portable  swimming  pool  placed  in  front  of  the  aircraft  w’hile  it 
was  parked  on  the  runway  of  the  Miami  International  Airport  in 
June  1967.  Measurements  were  made  of  the  reflected  sky  radiation, 
a  smooth  water  surface,  a  surface  of  bubbles,  and  a  surface  of 
foam.  His  values  are  shown  In  Table  XXTV. 

7 

This  quantity  is  discussed  in  detail  in  Chapter  VII. 
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Table  XXIV.  Results  of  ground  measurements  of  the  emission  of 

the  sky,  smooth  water,  bubbles,  and  foam  by  Williams 
(1968). 


Antenna  temperature,  Ta  (°K) 


Condition 

9.4  GHz 

15.8  GHz 

22.235  GHz 

34  GHz 

Reflected  Sky 

90 

120 

150 

170 

Smooth  Water 

170 

190 

218 

224 

Bubbly  water 

185 

210 

240 

260 

Foamy  Water 

295 

280 

285 

310 

Serious  doubt  must  be  cast  on  the  absolute  accuracy  of  these 
measurements  due  to  the  fact  that  the  measurements  were  made  in 
the  near  field  of  the  radiometers.  Note  the  unusually  high  values 
of  reflected  sky  temperature. 

It  is  evident  from  values  in  Table  XXIX  that  the  emissivity 
of  foam  approaches  unity.  Th’is,  as  Williams  proposed,  one  should 
be  able  to  measure  the  amount  of  foam  coverage  on  the  sea  sur¬ 
face  by  observing  the  increase  in  the  antenna  temperature  that 
results  when  a  large  area  is  viewed  by  a  microwave  radiometer  in 
which  there  exists  some  foam  covered  areas  and  some  open  areas. 

The  percentage  of  coverage  by  foam  is  related  to  the  surface  wind 
speed  (Munk,  1947;  Blanchard,  1963).  Foam  appears  to  raise  the 
apparent  brightness  temperature  of  the  sea  surface  significantly 
for  all  frequencies,  polarizations,  and  incidence  angles. 
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Ice 


There  have  been  many  excellent  papers  written  on  the  micro- 
wave  emission  of  ice.  The  paper  by  Pascalar  and  Sakamoto  (1965) 
is  quite  representative.  Two  cases  will  be  considered.  First, 
predictions  are  made  of  the  polarized  brightness  temperature  of  a 
large  ice  domain  such  as  that  of  an  iceberg.  Secondly,  the 
brightness  temperature  of  horizontal  sheets  of  ice  overlying 
water  is  studied. 

In  the  first  case,  one  may  consider  the  problem  of  determining 
the  reflectivity  of  a  homogenous  dielectric  having  a  flat  surface. 
The  dielectric  constant  of  ice  can  be  taken  to  be  3. 165- j 0.003 
(see  Chapter  II).  Using  (143),  one  finds  that  the  reflectivity 
of  ice  is  0.08  for  ij;  equal  to  zero.  Therefore,  the  emissivity 
of  ice  for  incidence  angles  near  0°  is  about  0.92.  If  one  consi¬ 
ders  the  difference  between  the  emission  of  ice  at  a  temperature 
of  273°  K  and  the  combined  emission  and  reflection  of  sea  water 
at  a  temperature  of  273°  K  and  for  a  frequency  of  9.2  GHz,  for 
example,  then  one  finds  that  the  polarized  brightness  temperatures 
are  251°  K  and  107°  K,  respectively.  The  difference  between  these 
two  values  is  144°  K  and  is  easily  measured  by  microwave  radio¬ 
meters  in  use  today.  In  fact,  scanning  microwave  radiometers 
operating  near  13  GHz  are  being  used  by  the  U.S.  Coast  Guard  to 
survey  icebergs  in  the  North  Atlantic  Oceans  (Roeder ,  1967). 


In  the  field,  the  microwave  radiometers  do  not  ’see'  as  large  a 
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brightness  temperature  difference  as  indicated  above  since  the 
iceberg  does  not  occupy  the  entire  field  of  view  of  the  sensor. 
The  size  of  the  increase  in  brightness  temperacure  is  used, 
qualitatively,  to  distinguish  between  large  and  small  icebergs. 

Since  the  imaginary  part  of  the  dielectric  constant  of  ice 
is  small,  one  would  expect  ice  to  be  nearly  lossless.  This 
property  of  ice  can  be  used  to  measure  remotely  the  thickness 
of  horizontally  uniform  ice  sheets  overlying  water.  The  reflec¬ 
tivity,  r*,  of  an  ice  sheet  of  thickness  2.  over  water  for  normal 
incident  is  easily  deduced  from  standard  text  books  on  electro¬ 
magnetic  theory  (for  example,  Ramo,  et.al. ,  1965)  and  may  be 
expressed  as: 


|<2L-21)/(Zl+Zl)|2 


(dimensionless),  (154) 


where  Z^  is  the  impedence  of  air  (377  ohms)  and  where: 


7  -  z  ^w  cos(2r£/lj)  4-  jZj  sin  ( 2tt £./ a -j ) 

L  iZ  .  cos(2it£/X.)  +  jZ  sin  (2tt£/X.) 
i  i  w  l 


(ohms) , 


Z 

w 


/p  /e  (e  '  -  je  ") 
0  0  w  w 


(ohms) , 


(155) 

(156) 


Z 


i 


e . 

l 


(ohms),  (157) 


X 


i 


=  c/v 


(meters),  (158) 


e  *  -  ie  "  is  the  complex  dielectric  constant  of  water,  and  £ ,  is 
w  J  w  1 

the  dielectric  constant  of  ice  (about  3.165). 


no 


One  should  note  that  the  arguments  of  the  trigonometric  func¬ 
tions  in  (155)  involve  the  thickness  of  the  ice  and  the  wave¬ 
length  of  the  electromagnetic  wave  in  the  ice  (X^,  meters). 

Since  these  functions  are  periodic,  one  would  expect  that  r*  is 
periodic  in  Thus,  the  polarized  brightness  temperature  is  not 
single  valued  with  respect  to  the  ice  thickness.  It  does  not 
appear  that  measurements  of  ice  thickness  could  be  made  by  a 
single  channel  microwave  radiometer.  This  could  be  done,  however, 
with  a  multi-frequency  radiometer  since  X^  would  change  with 
frequency. 
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CHAPTER  VI 

INTENSITY  OF  UPWELLING  MICROWAVE  RADIATION 

General 

The  purposes  of  this  chapter  are  to  illustrate  the  intensity 
field  of  upwelling  microwave  radiation  for  several  general  cases 
as  a  function  of  appropriate  physical  parameters  of  the  hydro¬ 
sphere  and  atmosphere  and  to  discuss  the  relative  advantages  and 
disadvantages  of  using  various  polarizations,  incidence  angles, 
and  microwave  frequencies  to  study,  remotely,  certain  aspects  of 
this  marine  environment.  One  should  recall  at  this  point  that 
intensity  is  a  quantity  that  is,  strictly  speaking,  impossible 
to  determine  with  practical  sensors.  Regardless,  the  concept  of 
intensity  must  be  retained  until  the  next  chapter  where  the 
effects  of  the  antenna  characteristics  of  the  microwave  radio¬ 
meter  used  to  measure  microwave  radiation  will  be  considered. 
Then,  the  intensity  field  may  be  weighted  by  the  gain  pattern  of 
the  antenna  to  predict  theoretical  values  of  actual  measurements 
by  microwave  radiometers. 

This  chapter  is  concerned  only  with  the  upwelling  microwave 
radiation  that  is  passing  through  some  point  above  the  surface  of 
Earth  in  the  atmosphere  where  aircraft  may  operate  or  in’free' 
space  where  an  Earth  satellite  may  operate. 

Let  T  (h  ,<J>)  be  the  polarized  brightness  temperature  of  the 
P  o 
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upwelling  microwave  radiation  at  height,  h  (meters),  and  for  an 

o 

incidence  angle  of  i|j.  In  Chapter  V,  (0,\J»)  was  discussed. 

In  the  case  of  a  non-scattering,  horizontally  uniform  atmosphere, 
it  follows  from  (87)  that: 

T  (h  ,<!/)  =  T  (0,«J.)  exp[-sec<|>  /h°  a(z)  dz] 

pop  o 

+  sect// /^o  a(z)  T(z)  exp[-seciii/ho  a(z)  dz]dz 

(°K).  (159) 

The  first  term  of  (159)  refers  to  the  portion  of  the  emitted 
and  reflected  microwave  radiation  that  reaches  height,  hQ,  after 
partial  extinction.  The  second  term  of  (159)  refers  to  the  direct 
contribution  from  the  atmosphere  between  hQ  and  the  surface. 

One  may  immediately  describe  two  special  cases.  In  the  first 
case,  hQ  is  above  20  kilometers.  In  this  case,  the  radiometer  is 
above  most  of  the  radiating  part  of  the  atmosphere.  Due  to  the 
high  transmissivity  of  the  atmosphere,  one  may  assume  that  the 
atmosphere  emits  upward  at  the  same  intensity  that  it  emits  down¬ 
ward  for  the  same  incidence  angle.  The  author  has  computed  the 
brightness  temperature  of  upwelling  microwave  radiation  at  19.4 
GHz  using  Stogryn's  (1967)  values  for  fractional  absorption  and 
found  that  it  differs  by  only  a  few  nvmdreths  of  a  degree  from  the 
sky-temperature  calculated  for  the  same  incidence  angles.  Thus, 
cor  a  clear  atmosphere,  it  is  assumed  that: 

TaJr  (v,*,hQ  1  20  km)  =  TSJy  (v,*)  (°K)  , 


(160) 


Ill 


sir 

where  T  represents  the  second  term  of  (159/  in  terms  of  an 

r 

aoparent  brightness  temperature  (°K),  and  wl.t  re  TS^  is  the 

p 

polarized  sky  temperature  as  used  in  Chapter  IV.  The  angle, 
ip,  represents  both  the  incidence  angle  and  angle  of  observation 
since  they  are  equal  interior  angles. 

In  the  second  case,  hQ  is  very  small  (less  than  one  kilo¬ 
meter).  In  this  case,  one  may  treat  the  intervening  atmosphere 
as  an  isothermal,  homogenous  layer  and  compute  its  emission  by: 


7  W,h  <  1  km)  T  [1  -  exp  (-ah  ]  (°K),  (161) 

p  o  —  o 

where  T  is  the  mean  temperature  of  the  layer  (°K)  and  where  a  is 
the  mean  value  of  the  absorption  coefficient  of  the  layer. 

Five  special  cases  will  be  considered  in  the  next  sections. 
These  cases  are  listed  in  Table  XXV.  The  purpose  in  presenting 
these  cases  is  to  illustrate  the  effect  of  salinity,  temperature, 
roughness,  incidence  angle,  polarization,  frequency,  humidity, 
clouds,  and  altitude  on  the  field  of  brightness  temperature  at 
some  point  above  Earth.  The  effect  of  temperature  as  a  variable 
has  been  considered  before  in  Chapter  V  (see  Figure  8) . 


Salinity  Effect 


One  of  the  most  important  parameters  of  the  hydrosphere  is  its 
salinity.  The  salinity  of  the  surface  layer  of  the  sea  is  usually 
34 to  37  o/oo.  Near  coastlines,  the  salinity  may  vary  between  5  o/oo 
and  37  o/oo  due  to  the  influence  of  coastal  runoff  of  precipitation 


Table  XXV.  List  of  special  cases  considered  in  Chapter  VI. 
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or  river  outflow.  Temperatures  in  the  ocean  are  usually  between 
•1°  C  and  30°  C.  To  illustrate  the  effect  of  changes  in  salinity 
on  the  polarized  brightness  temperature,  the  following  environ¬ 
mental  conditions  have  been  assumed: 

1.  Let  the  temperature  of  sea  water  be  held  at  10  °C. 

2.  Assume  that  the  sea  surface  is  flat,  smooth,  and  non- 
f oamy . 

3.  Let  the  incidence  angle,  i|>,  be  zero  so  that  the  maximum 
effect  of  salinity  may  be  ascertained  (Sirounian,  1968) .  The 
previous  studies  of  the  effect  of  roughness  show  that  the  eiiis- 
sivities  of  a  rough  and  smooth  surface  for  near  zero  degrees  are 
practically  identical  as  long  as  there  is  no  foam  on  the  surface. 

4.  Assume  that  there  are  no  clouds  present  in  the  atmosphere 

and  that  the  atmosphere  is  in  a  standard  state  as  far  as  its 

temperature,  humidity,  and  pressure  dependence  or  height  are 

concerned . 

5.  Let  the  height  of  the  observer  be  one  kilometer. 

6.  Let  the  salinity,  S,  vary  from  0  o/oo  to  36  o/oo. 

Figure  10  shows  plots  of  T  (v,h  =  1,  <Ji  *  0,  S,  T  *=  10  C°) 

p  o 

versus  S  for  frequencies  of  1,  5.4,  9.3,  15.8,  19.35,  22.235,  and 
34  GHz.  The  polarized  sky  temperatures  have  been  taken  from 
Table  XIII.  The  specific  values  used  in  the  calculations  are 
given  in  Table  XXVI.  The  mean  temperature  of  the  layer  of  atmo¬ 
sphere  between  the  surface  and  one  kilometer  is  taken  to  be  283°  K. 
The  absorption  coefficient  for  molecular  oxygen  and  water  vapor 


12  f» 


are  taken  from  Table  X  for  a  density  of  water  vapor  of  7.5  grams 
per  cubic  meter.  These  values  are  also  shown  in  Table  XXVI. 

The  values  of  the  transmissivity  and  enissivity  of  the  layer  of 
atmosphere  from  the  surface  to  a  height  of  kilometer  a.'*e  also  shown 
in  Table  XXVI  as  they  are  derived  Trom  the  mean  absorption  coef¬ 
ficient. 

It  can  be  seen  from  Figure  10  that  salinity  affects  the 
polarized  brightness  temperature  significantly  for  a  frequency 
of  1  GHz  and  insignificantly  at  the  other  frequencies  shown. 

Calculations  were  made  of  che  nolarized  brightness  temperature 
under  the  same  conditions  as  above  with  two  exceptions:  the 
temperature  of  the  sea  water  was  assumed  to  be  30°  C,  and  the  mean 
temperature  of  the  layer  from  the  surface  to  one  kilometer  was 
assumed  to  be  294°  K.  Tne  values  of  sky  temperature  also  were 
increased  to  more  representative  values  for  humid,  summertime 
conditions  over  the  tropical  ocean.  The  results  of  the  calcu¬ 
lations  for  this  case  are  shown  in  Figure  11. 

One  can  see  that  the  salinity  effect  at  1  GHz  is  much  greater 
for  a  water  temperature  of  30°  C  than  for  a  temperature  of  10°  C. 
Again,  the  polarized  brightness  temperatures  for  the  other  fre¬ 
quencies  shown  in  Figure  11  are  practically  constant  with  salinity. 
Thus,  it  appears  that  one  may  be  able  to  determine  the  temperature 
and  salinity  of  the  sea  surface  to  some  degree  of  accuracy  by 
measuring  the  intensity  of  microwave  radiation  at  5.4  GHz  and  1.0 


Table  XXVI.  Assumed  atmospheric  parameters  and  derived  quantities  used  in  calculating  values 
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GHz.  The  existence  of  significant  foam  cover  and  rain  would 
negate  this  possibility. 


Effect  of  Water  Vapor 


The  contribution  of  thermal  emission  by  water  vapor  in  the 
atmosphere  to  the  overall  polarized  brightness  temperature  is 
quite  large  for  frequencies  near  the  resonance  line  of  water 
vapor.  If  one  assumes  a  certain  form  for  the  distribution  of 
water  vapor  with  height,  it  is  possible  to  compute  values  of 
polarized  brightness  temperature  for  an  incidence  angle  of  zero 
degrees  for  various  conditions  of  humidity. 

Suppose  one  wishes  to  measure  the  total  amount  of  precipi- 
table  water  (W,cm)  in  the  marine  atmosphere.  It  has  been  shown 
in  Chapter  V  that  the  emissive  temperature  of  sea  water  is 
practically  constant  for  the  usual  ranges  of  temperature  and  sa¬ 
linity  for  ;  requencies  of  19.35  GHz  and  22.235  GHz.  It  is  assumed 
that  is  120°  K  and  122°  K,  respectively.  The  reflectivities 
for  these  frequencies  are  approximately  0.60  and  0.59,  respectively. 
Kreiss  (1968)  shows  that  the  zenith  sky  temperature  ranges  from 

9  to  43°  K  for  19.35  GHz  and  from  14  to  106°  K  for  22.235  GHz  under 
various  meteorological  conditions.  Thus,  the  overall  brightness 
temperature  of  upwelling  microwave  radiation  may  change  by  as 
much  as  147°  K  due  to  changes  in  the  emission  of  water  vapor  in 
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0°  for  several  conditions  of  temperature  and  humidity  in  the 
atmosphere  (see  Table  XIX) .  Using  his  values  of  precipitable 
water  for  these  conditions  of  temperature  and  humidity,  one  may 
plot  Tp  (v,hQ  L  20^an»  4*  =  0)  versus  W  for  frequencies  of  19.35 
GHz  and  22.235  GHz  as  shown  in  Figure  12  (circles).  An  indepen¬ 
dent  calculation  by  Staelin  (personal  communication,  1967)  is 
also  included  in  Figure  12  (triangles). 

There  appears  to  be  a  simple  correlation  between  these 
polarized  brightness  temperatures  and  the  total  amount  of 
precipitable  water  in  the  atmosphere.  The  values  given  by 
Kreiss  were  computed  for  nine  conditions  of  temperature  and 
humidity  in  the  atmosphere;  however,  the  effect  of  temperature 
does  not  appear  to  be  significant.  This  fortunate  circumstance 
probably  is  due  to  the  fact  that  the  amount  of  moisture  present 
in  the  atmosphere  is  highly  dependent  upon  its  temperature  pro¬ 
file.  The  values  shown  in  Figure  12  are  for  hQ  above  20  kilo¬ 
meters.  As  the  observer  decends  to  lower  altitudes,  the  polarized 
brightness  temperature  decreases  proportionately.  Even  at  the 
surface,  one  still  would  be  able  to  distinguish  between  various 
conditions  of  precipitable  water  due  to  the  facet  that  the  sur¬ 
face  of  the  water  reflects  about 
upon  it. 


Effect  of  Roughness 


Stogryn  (1S67)  has  presented  a  very  complete  treatment  of  the 


POLARIZED  BRIGHTNESS  TEMPERATURE  OF  SEA  WATER,  TD  (°K) 
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effects  of  a  distribution  of  wave  slopes  on  the  polarized 
brightness  temperature.  Figure  13  shows  the  plots  of  polarized 
brightness  temperature  versus  incidence  angle  for  horizontally 
and  vertically  polarized  radiation  at  19.4  GHz  for  the  following 
conditions: 

1.  The  temperature  of  the  sea  water  is  held  constant  at 

17°  C,  and  the  salinity  of  the  sea  water  is  assumed  to  be  36o/oo. 
One  may  say  that  these  conditions  are  representative  of 
the  average  surface  state  of  the  ocean  in  mid-latitudes.  It  can 
be  seen  in  Figure  8  that  the  emissive  temperature  of  the  sea  is 
fairly  constant  with  temperature  and  salinity  at  19.4  GHz. 

2.  The  temperature  profile  of  the  atmosphere  was  assumed 
to  be  that  of  the  ARDC  standard  atmosphere. 

3.  The  density  of  water  vapor  was  assumed  to  vary  with 
altitude  as  given  by  (104) . 

4.  The  Kirchhoff  Approximation  was  assumed  to  be  valid, 
and  the  slope  statistics  developed  by  Cox  and  Munk  (1954a,  1954b) 
were  adopted. 

Values  of  and  are  plotted  in  Figure  13  versus  iji  for 
hQ  equal  to  one  kilometer  and  for  w  equal  to  0,  4,  8,  and  14 
meters  per  second. 

Figure  13  shows  that  the  roughness  effect  is  greatest  for 
horizonttl  polarization  and  for  incidence  angles  near  50°.  The 
main  cause  of  this  phenomenon  is  the  change  in  the  source  of  the 
reflected  sky  radiation  that  results  when  the  surface  is  roughened 


POLARIZED  BRIGHTNESS  TEMPERATURE  OF  SEA  WATER 


INCIDENCE  ANGLE,  degrees 

Figure  13.  The  polarized  brightness  temperatures  of  sea  water  versus 
incidence  angle  for  horizontal  and  vertical  polarization, 
for  fully  developed  sea  driven  by  surface  winds  of  0,  A, 
and  1A  m/sec,  for  an  altitude  of  1  km,  and  for  a  frequency 
of  19. A  GHz  ‘after  Stogryn,  1967). 
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(Stogryn,  1967). 

It  has  been  noted  before  that  the  root-mean-squared  slopes. 

and  g^,  of  the  sea  surface  are  related  to  wind  speed  by  (152) 

and  (153).  Saunders  and  Wilkens  (1968)  say  that  if  the  sun  is 

o. erhead,  the  angular  radii  (radians)  of  the  sun's  glitter  pattern 

as  seen  in  photographs  are  approximately  2g  and  2g  along  direc- 

x  y 

tions  cross  wind  and  upwind,  spectively.  Thus,  one  may  use 
photographs  of  the  sun's  glitter  pattern  to  provide  'ground  truth' 
about  the  roughness  of  the  sea  surface  (Cox  and  Munk,  1954a,  1954b). 

Table  XXVIJ  shows  the  difference  between  horizontally  polar¬ 
ized  brightness  temperatures  for  incidence  angles  of  0°  and  45° . 

These  theoretical  values  of  T^  as  a  function  of  roughness 
are  valid  only  if  there  is  no  foam  on  the  sea  surface.  For  surface 
winds  above  seven  meters  per  second,  one  would  expect  some  finite 
area  of  the  sea  surface  to  be  covered  by  foam  and  bubbles. 


Table  XXVII.  Difference  between  the  horizontally  polarized 

brightness  temperature  of  sea  water  at  an  incidence 
angle  of  0°  (T^  (0°))  ind  that  of  sea  water  at  an 
incidence  angle  of  45°  (Tp  (45°))  for  a  frequency 
of  19.4  GHz  versus  the  speed  of  the  wind  at  41  feet 
above  sea  level  (w,  m/sec)  for  3  fully  developed  sea. 


Wind  speed,  w,  (m/sec) 


0 

4 

8 

14 


Th  (*  =  0)  -  Th  (*  =  45°) 


19°  K 
14°  K 
8°  K 
1°  K 
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Effect  of  Liquid  Water 

Kreiss  (1968)  has  computed  values  of  the  polarized  brightness 

of  upwelling  microwave  radiation  over  various  types  of  clouds 

3 

containing  various  concentrations  of  liquid  water  (M,  g/m  ). 

These  computations  are  presented  in  his  paper  for  specific  cases 

and  values  of  M  as  T  (v,h  >30  km,  i b  =0)  versus  v  for  v 

p  o 

between  15  GHz  and  30  GHz.  It  may  be  more  convenient  to  present 
his  data  in  a  slightly  different  manner. 

Consider  the  following  case: 

1.  The  frequency  is  19.35  GHz. 

2.  Standard  temperature  and  humidity  profiles  as  given  by 
Kreiss  (1968)  are  assumed  to  be  valid.  In  addition,  the  relative 
humidity  is  assumed  to  be  100%  in  clouds. 

3.  Equation  (139)  is  assumed  to  be  valid. 

4.  The  cloud  bases  are  assumed  to  be  688  meters  above  the 
sea  surface. 

5.  Clouo  tops  are  assumed  to  be  at  1949  meters  for  'stratus', 
3012  meters  for  'thick  stratus'  and  7185  meters  for  'tall  cumulus' 
clouds. 

Figure  14  shows  plots  of  T  (19.35  GHz,  h  >20  km,  =  0) 

p  o 

versus  M  for  'stratus',  'thick  stratus',  and  'tall  cumulus’  clouds. 

Figure  15  shows  a  plot  of  ?  versus  the  product  of  M  times 
the  thickness  of  the  cloud.  This  latter  product  m3v  be  called 
the  concentration  of  liquid  water  per  unit  horizontal  area. 


POLARIZED  BRIGHTNESS  TEMPERATURE  OF  SEA  WATER,  T  (°K) 
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of  0°. 


T 


5  10  15  20  25 

TOTAL  VERTICAL  CONTENT  OF  LIQUID  WATER,  M5hQ  (1000  g/m2) 
Figure  15.  The  polarized  brightness  temperature  of  sea  water  versus 
the  product  of  the  thickness  of  intervening  clouds  and 
their  concentration  of  liquid  water  for  a  frequency  of 
19.35  GHz  and  for  an  incidence  angle  of  0°. 
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The  dashed  lines  in  Figure  34  indicate  the  temperature  of  the 
top  layers  of  the  stratus  (A),  thick  stratus  (B) ,  and  tall  cumulus 
(C).  It  is  apparent  from  Figure  14  that  the  polarized  brightness 
temperature  over  clouds  for  a  frequency  of  19.35  GHz  depends  on  the 
type  of  cloud  and  its  concentration  of  liquid  water.  The  plots 
in  Figure  15  indicate,  furthermore,  that  the  temperature  of  the 
cloud  is  very  important  in  determining  the  microwave  radiance 
over  the  cloud.  In  this  figure,  the  effects  of  total  vertical 
water  content  are  normalized  by  using  M  times  AhQ  as  the  indepen¬ 
dent  variable. 

Of  course,  if  one  knew  the  thickness  and  height  of  the  clouds 
being  viewed,  one  could  use  curves  in  Figure  14  to  ascertain  its 
mean  value  of  H  from  measurements  of  the  brightness  temperature. 

It  appears  that  one  may  be  able  to  determine  the  hei^-  .t  of  the 
topmost  clouds  in  the  atmosphere  by  measuring  the  infrared  radia- 
tional  intensity  near  wavelengths  of  10  microns  as  by  measurements 
of  the  amount  of  attenuation  of  near  visible  radiation  due  to 
carbon  dioxide  (Hanel,  1961). 

Kreiss  (1968)  assumed  that  the  sea  surface  acts  as  a  smooth, 
flat  surface  —  even  under  raining  clouds.  Rain  tailing  on  the 
sea  surface  creates  bubbles  and  foam.  This  would  tend  to  raise 
the  radiance  of  the  cloud  at  microwav.  frequencies. 
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CHAPTER  VII 
MICROWAVE  RADIOMETRY 

General 

The  historical  develonment  of  microwave  radiometry  from 
radio  astronomy  to  a  discipline  of  terrestrial  remote  sensing 
was  given  in  Chapter  I.  The  three  basic  components  of  a  micro- 
wave  radiometer  are  its  antenna  system,  its  receiver,  and  its 
output  recorder . 

The  antenna  system  intercepts  a  portion  of  the  total  electro¬ 
magnetic  radiation  incident  upon  it  and  transfers  this  power  to 
the  input  terminals  of  the  receiver  along  connecting  transmission 
lines  or  wave  guides.  Antennas  used  in  microwave  radiometry  are 
highly  directive,  that  is,  they  receive  energy  mostly  from  a  cer¬ 
tain  set  of  directions  relative  to  the  axis  of  the  antenna. 

The  power  supplied  to  the  input  terminals  of  the  receiver 
is  greatly  amplified  and  filtered  so  that  only  power  having  fre¬ 
quencies  from  v  -  A v/2  to  v  +  Av/2  is  passed  to  its  output  ter¬ 
minals.  Av  is  the  band  pass  of  the  radiometer,  and  v  is  the  cen¬ 
tral  frequency  of  the  radiometer. 

The  output  from  the  amplifier  is  converted  to  direct  current 
and  may  be  recorded  on  a  strip  chart  or  on  magnetic  tape. 

Since  microwave  radiation  is  weak,  the  contribution  of 
.uermal  emission  in  the  antenna,  transmission  lines,  and  first 
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stages  of  amplification  may  be  comparable  to  the  radiational 
power  intercepted  by  the  antenna. 

The  purpose  of  this  chapter  is  to  present  certain  concepts 
necessary  to  allow  one  to  relate  the  field  of  polarized  brightness 
temperatures  present  at  some  location  in  space  to  measurements 
made  by  typical  microwave  radiometers.  The  concepts  presented 
in  the  next  sections  are  discussed  in  great  detail  by  Bracewell 
(1962)  in  his  paper  on  techniques  of  radio  astronomy. 

Microwave  Receivers 

The  functions  of  the  receiver  portion  of  a  microwave  radio¬ 
meter  are  to  amplify  the  power  supplied  to  its  input  terminals  by  the 
antenna  system,  to  convert  this  power  to  more  suitable  frequencies, 
and  to  tune  the  output  to  a  certain  band  of  frequencies.  It  is 
important  to  remember  that  the  radiometer  is  not  monochromatic. 

A  typical  bandwidth  of  microwave  radiometers  used  today  is 
0.1  GHz. 

A  block  diagram  of  a  typical  microwave  receiver  of  a  Dicke 
radiometer  is  shown  in  Figure  16. 

The  antenna  system  supplies  a  small  amount  of  power  to  the 

input  terminals  of  the  microwave  receiver.  This  power  is  typically 
—20 

10  watts  per  Hertz  and  is  carried  by  a  wide  range  of  frequen¬ 

cies.  The  antenna  system  displays  some  degree  of  selectivity  and 
does  not  deliver  all  of  the  power  incident  upon  its  physical  are.’ 
to  the  input  terminals.  This  tuning  effect  is  small,  however. 
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compared  to  tuning  that  occurs  in  the  detectors  of  the  microwave 
receiver.  The  antenna  system  does  exhibit  one  important  aspect  of 
selectivity,  viz.,  it  chooses  a  certain  polarization.  If  the  inout 
power  is  essentially  constant  with  frequency,  one  defines  the 
bandwidth  (Av,  Hertz)  as  "the  difference  between  the  frequencies 
at  which  the  output  voltage  (or  other  output  indication)  falls  to 
the  value  corresponding  to  a  halved  power  input  at  the  mid¬ 
frequency,  v "  (Bracewell,  1962)  .  The  shape  of  the  band  pass  is 
usually  assumed  to  be  Gaussian. 

Let  be  the  power  supplied  by  the  antenna  system  to  the 
input  terminals  and  carried  by  frequencies  from  v  -  Av/2  to 
v  +  Av/2.  Formally,  one  may  say  that  the  actual  power  passed  by 
the  reciver  is  times  a  constant  factor  since  power  is  passed 
outside  Av  and  since  not  all  the  power  existing  in  Av  is  passed. 

On  the  other  hand,  calibrations  of  the  radiometer  consist  of 
comparing  this  power  to  the  power  produced  by  a  matched  resistor 
or  produced  when  the  antenna  is  surrounded  by  a  blackbody.  One 
would  expect  the  same  biasing  effect  on  the  calibration  power 
as  on  the  power  supplied  by  the  antenna  system.  It  is  further 
assumed  that  P^  is  constant  over  Av. 

Microwave  frequencies  are  much  too  high  to  be  handled  by 
ordinary  components  of  radio  amplifiers.  If  v  is  less  than  10  GHz, 
however,  Pj  may  be  araolified  before  the  signal  is  converted 
to  lower  frequencies  by  mixing  the  original  signal  with  a  lower 
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frequency  signal. 

The  video  amplifier  then  boosts  the  power  to  a  level  that  is 
recordable.  The  signal  is  rectified,  integrated,  and  fed  to  the 
output  recorder. 

The  most  serious  limitation  placed  on  early  microwave  receiv¬ 
ers  was  drift  in  amplification.  Dicke  (1946)  proposed  a  method 
whereby  the  effect  of  relatively  slow  fluctuations  in  amplification 
could  be  lessened.  The  power  from  a  fixed,  known,  internal  source 
is  fed  into  a  switch  which  alternately  passes  power  from  the 
antenna  and  the  internal  source  at  a  rapid  pace.  Any  drift  in 
amplification  affects  the  amplification  of  the  antenna  power  and 
the  internal  power  in  the  same  amount;  therefore,  the  difference 
between  these  signals  should  not  be  affected  by  amplification 
fluctuations.  The  signal  from  the  reference  generator  which 
drives  the  Dicke  switch  also  is  supplied  to  the  final  detector 
stage.  The  antenna  power  is  then  recovered. 

Since  a  continuous  spectrum  of  frequencies  is  passed  by  the 
microwave  receiver,  modulations  are  evident  on  the  wave  form  of 
che  amplified  signal.  Bracewell  (1962)  shows  that  the  root- 
mean-square  deviation  in  the  output  voltage  is  proportional  to 
the  mean  value  of  the  output  voltage  divided  by  the  square  root 
of  the  integration  time  multiplied  by  the  bandwidth. 


The  power 


wafts)  produced  by  a  matched  resistor  at 


temperature,  T  (°K),  over  the  frequency  range  from  v  -  flv/2  to 
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V  +  Av/2  is  called  Johnson  noise  and  is  given  by  Bracewell  (1962) 
as : 


Pj  =  kT  Av  (°K) ,  (162) 

where  k  is  Boltzmann's  constant. 

Using  (162)  it  is  possible  to  express  in  terms  of  the 
temperature  of  a  matched  resistor  (T . ,  °K)  which  produces  a  power 
P^.  Thus, 

T  =  P±  /  kAv  (°K).  (163) 

is  called  the  apparent  input  temperature. 

The  value  of  T^  is  determined  by  the  total  amount  of  power 
received  by  the  antenna  system  of  the  radiometer  and  supplied  to 
the  input  terminals.  A  significant  portion  of  may  originate 
in  the  antenna  system  and  connecting  transmission  lines  since  these 
devices  also  emit  Johnson  noise. 

The  output  of  the  radiometer  may  be  expressed  in  terms  of 
apparent  input  temperature.  From  the  preceeding  discussion  of 
voltage  accuracies,  one  may  say  that  the  root-mean-squarc  deviation 
of  the  measured  apparent  input  temperature  AT^  is  given  by: 

AT (T  +  T.)/(t.  Av)1/2  (°K),  (164) 

1  Nil 

where  is  the  contribution  to  the  overall  signal  due  to  noise 
emission  in  the  radiometric  components  (°K) ,  and  t^  is  the 


145 


Integration  time  (seconds) . 

is  called  the  apparent  temperature  of  the  receiver  noise 
and  may  be  found  in  the  literature  in  terms  of  the  noise  figure 
(NF)  where: 

T„  =  (NF  -  1)290  (°K).  (165) 

Currently  used  microwave  receivers  have  noise  figures  close 
to  unity. 

Typical  values  of  AT^  are  5°  K  in  airborne  microwave  radio¬ 
meters.  The  next  generation  of  airborne  microwave  radiometers 
should  have  values  of  AT^  near  0.5°  K.  Some  radio  telescopes 
have  values  of  AT^  near  0.01°  K. 

Now  it  is  necessary  to  relate  to  the  field  of  brightness 
temperature  and  certain  antenna  parameters. 

Antenna  Effects 

It  is  convenient  to  define  a  new  set  of  direction  angles 
relative  to  the  axis  of  a  high  gain  antenna.  Let  (9  ,<{>)  be  this 
set  as  shown  in  Figure  17.  The  angle,  9,  lies  in  the  vertical 
plane  of  polarization,  and  4>  lies  in  the  horizontal  plane  of 
polarization.  The  axis  of  the  antenna  may  be  inclined  to  the 
vertical  by  an  angle  \|i  which  is  called  the  angle  of  incidence  or 
angle  of  viewing.  If  is  zero,  the  vertical  and  horizontal 
planes  lose  their  significance.  In  this  case  the  'vertical  plane' 
is  defined  as  the  vertical  plane  that  contains  the  path  of  the 
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platform  carrying  the  microwave  radiometer.  In  most  cases,  the 
field  of  intensity  is  symmetrical  to  the  normal. 

The  characteristic  of  an  antenna  system  are  described  in 
terms  of  its  radiating  properties.  Let  P  be  the  total  power 
radiated  by  an  antenna.  It  follows,  then,  that  the  average 
power  (P,  watts/steradian)  radiated  per  unit  solid  angle  is  Pt/4u 
Antennas  are  usually  highly  directive,  that  is,  they  radiate  un¬ 
evenly  with  (0,<*>).  Let  P(9,<J>)  describe  the  distribution  of  the 
radiated  power  where: 

P  =  /  P(0,$)  dfi  (watts).  (166) 

C  4lT 

If  P^^  is  the  power  supplied  to  the  antenna  terminals,  then 
c.ue  defines  an  mtenna  efficiency  factor,  n  ,  such  that: 

cl 

n  =  P  /  P.  (dimensionless).  (167) 

a  t  x 

Let 

D{6 , d>)  *  P(0,O)  /  P  (dimensionless),  (168) 

and 

G(6,<J>)  =  P(0,<j>)  4tt  /  P^  (dimensionless).  (169) 

D(0,4>)  is  called  the  directive  gain  pattern  of  the  antenna, 

and  G(0,4>)  is  called  the  gain  pattern  of  the  antenna. 

From  (167),  (168),  and  (169)  it  follows  that: 
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G(0,<j>)  =  n  D(0,$)  (dimensionless).  (170) 

Let  (0,<(>)  be  the  distribution  of  the  polarized  intensity 
of  the  upwelling  radiation  with  respect  to  the  antenna  axis. 

The  total  amount  of  power  (P  ,  watts)  in  polarization  p  inter- 
cepted  by  the  antenna  is  given  from  (39)  by: 

Pa  =  4ir  £v  aiea  Tp  cose '  dA  dv  d^  (watts).  (171) 

It  is  assumed  tha  I  (0,4>)  is  constant  over  the  area  of  the 
antenna  and  over  the  bandwidth.  The  angle  6'  is  the  angle  be¬ 
tween  the  axis  of  the  antenna  and  the  direction,  (0  ,<{>).  The  ef¬ 
fective  area  of  the  antenna,  (Ag,  square  meters)  is  introduced 
in  (171)  and  represents  the  efficiency  of  the  antenna  in  gathering 
electromagnetic  energy  as  well  as  the  lessening  of  the  physical 
area  of  the  antenna  due  to  the  factor,  cos0'.  Thus,  (171) 
becomes : 

P  «  I  I  (9 ,<J>)  A(0,<f>)  Av  dn  (watts).  (172) 

3  nil  p  6 

According  to  Jordon  (1950) , 

Ae  (0,<fr)  -  (c2/4ttv2)  G  (0,*)  (ra2),  (173) 

where  c  is  the  speed  of  light  in  a  vacuum. 

Using  (54)  to  express  I  (0 ,4>)  in  terms  of  T  (0,<!>)  and 

P  P 

using  (173),  one  obtains  the  following  expression  for  P  : 

3 

P  »  (kAv/47r)  /  T  (9,4>)  G  (0,<j>)  dfi  (watts).  (174) 

3  *t7T  p 


1  /'rt 
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Using  (162),  one  may  express  P  in  terms  of  the  temperature 

3 

of  a  matched  resistor,  (T  ,  °K)  or  simply  antenna  temperature 

3 

where: 

P  ■  kAvT  (watts) .  (175) 

3  3 

Thus,  (174)  becomes: 

Ta  «  (1/4*)  ^  Tp  (6,0  G  (9,0  da  (°K).  (176) 

When  a  radome  is  used  to  protect  the  antenna,  one  must  con¬ 
sider  the  .effects  of  absorption  and  emission  of  energy  by  the 
radome.  If  is  the  transmissivity  of  the  radome,  then: 

Ti'V4*)  L\  (6,0  G  (6,0  dP:  +  TR  (1  -  1^) 

(°K),  (177) 

where  TR  is  the  effective  thermometric  temperature  of  the  radome 
and  is  the  apparent  input  temperature  as  measured  by  the  radio¬ 
meter.  In  practice,  and  TR  have  been  difficult  to  determine. 

The  shape  of  G(0,O  is  usually  assumed  to  be  Gaussian  for 
(8,0  near  (0,0''.  In  this  case,  one  may  describe  G(0,O  hy 
giving  G(0,0),  and  the  beamwidths,  0  and  The  beamwidths  are 
the  angles  between  values  of  6  and  where  the  gain  is  equal  to 
half  of  the  gain  along  the  axis.  In  most  cases,  the  gain  is 
given  in  terms  of  decibel  reduction  relative  to  G(0,0).  Thus, 
the  gain  is  3  decibels  below  G(0,0)  when  6  =  ±0  or  I  =  ±$ 

In  some  cases  the  angular  widths,  0.,  and  <{>  ,  of  the  solid 
angle  containing  the  main  lobe  of  the  gain  pattern  are  given. 
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The  power  delivered  to  the  input  terminals  of  the  receiver 
by  the  antenna  comes  from  all  portions  of  the  4m  steradians 
surrounding  a  location  is  space.  Thus,  one  cannot  assume,  sum¬ 
marily,  that  the  radiometer  is  'seeing'  one  particular  area  under 
view.  This  is  especially  important  when  Tp  varies  strongly  with 
0  and  6.  Over  the  ocean,  however,  it  has  been  shown  in  past 
chapters  that  Tp  varies  slowly  with  $  especially  for  ^  near 
zero  degrees. 

Consider  the  case  of  an  antenna  placed  in  a  constant  tempera¬ 
ture  enclosure.  In  this  case: 

Ta  *  (Tp/ Air)  ^  G  (e,*)  dfi  *  Tp  (°K),  (178) 

since 

G  (0 ,4>)  da  =  Air.  (179) 

In  a  similar  fashion,  one  may  break  the  integral  in  (176) 
into  parts  and  say  that: 

Ta-  <Tp(0,0)/4„)  «,*)  dSl 

lobe 

+[Tplde/4*|  side  g  «,♦>  da 

'  lobes 

+  T£ack/47rj  G  (9,4>)  dfi  (°K) .  (180) 

lobes 

In  each  term  of  (180),  the  integrals  have  the  dimensions  of 
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solid  angle  and  may  be  referred  to  as  the  equivalent  solid  angle 
of  the  main  lobe  (Q^),  the  side  lobes  (fig),  and  the  back  lobes 
(flg),  respectively. 

In  most  cases,  /  4tt  is  0.7  to  0.9.  It  will  be  assumed 
that  the  antenna  temperature  is  equal  to  the  brightness  temperature 
along  the  antenna  axis  as  a  first  approximation.  The  mean  bright- 

— g 

ness  temperature  for  side  lobes  (T  ,  °K)  is  usually  as  high  or 

higher  than  for  ij >  near  zero  degrees,  and  the  mean  brightness 

temperature  for  back  lobes  (T  ,  °K)  is  much  lower  than  T  (0,0). 

P  P 

An  excellent  discussion  of  these  antenna  parameters  is  given 
by  Ko  (1964)  for  several  specific  cases. 

In  summary,  one  should  note  that  is  measured  by  the  radio¬ 
meter.  If  one  knows  L  and  G  (6,<*>),  it  is  still  impossible  to 

S 

obtain  (0,<>)  without  making  certain  simplifying  assumptions. 

In  the  case  of  measurements  over  the  ocean,  it  is  assumed  that 

does  relate  only  to  the  value  of  at  (e,$)  equal  to  zero,  i.e.. 


T  (♦)  -  [T,  -  TR  (1  -  V)  /  Lr 


<"K),  (182) 


where  is  the  calibrated  reading  of  the  radiometer. 
If  no  radcme  is  used,  it  is  assumed  that: 


T  W  =  T  (*)  (°K) .  (183) 

P  a 

In  any  event,  one  would  be  able  to  measure  relative  changes 
in  apparent  antenna  temperature  that  occur  along  a  flight  line 


%  r  A 


and  assume  that  the  corresponding  change  in  relative  polarized 
brightness  temperature  would  be  comparable. 


CHAPTER  VIII 


RESULTS  OF  MICROWAVE  MEASUREMENTS 

General 

Few  measurements  have  been  made  of  microwave  radiance  over 
sea  water  under  conditions  in  which  the  distributions  of  essential 
atmospheric  and  hydrospheric  parameters  were  known  adequately. 
Indeed,  it  may  be  impossible  to  determine  the  distribution  of 
liquid  water  in  a  cloud  field  by  in  situ  measurements.  Atmo¬ 
spheric  temperature  and  humidity  soundings  often  are  not  available 
over  the  ocean.  Surveys  of  the  ocean's  surface  by  research  ves¬ 
sels  do  not  result  in  truly  synoptic  measurements.  The  nature 
of  the  small-scale  roughness  of  the  sea  surface  is  difficult  to 
deduce  from  standard  wave  instruments. 

The  National  Aeronautics  and  Space  Administration  (NASA) 
has  been  actively  sponsoring  research  and  development  in  applied 
remote  sensing.  Under  NASA's  Earth  Resources  Survey  Program, 
many  airborne  tests  have  been  made  of  a  set  of  remote  sensors 
that  will  probably  be  placed  in  a  spacecraft  in  the  1970's. 
Applications  are  being  developed  under  this  program  for  all  of 
the  earth  sciences.  At  the  present  time,  the  Goddard  Space  Flight 
Center  and  the  Environmental  Sciences  and  Service  Administration 
are  actively  sponsoring  research  in  remote  sensing  for  meteoro¬ 
logy.  It  is  really  impossible  to  separate  meteorology  from  the 
other  earth  sciences  since  one  must  always  look  through  the 
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atmosphere  when  one  is  trying  to  study  some  aspect  of  the  land 
or  the  ocean  from  space.  Over  the  ocean,  one  must  always  be  aware 
of  the  mutual  influences  between  the  atmosphere  and  hydrosphere. 

The  Spacecraft  Oceanography  Project  has  been  sponsoring 
studies  at  Texas  A&M  University  of  possible  applications  of 
remote  sensing  over  the  Mississippi  Delta  as  a  part  of  NASA's 
Earth  Resources  Survey  Program.  One  of  the  objectives  of  this 
study  is  to  determine  whether  or  not  remote  sensors  can  be  used 
from  aircraft  and  spacecraft  to  study  the  distribution  of  the 
outflow  of  the  Mississippi  River  into  the  Gulf  of  Mexico.  The 
chief  indicators  of  outflow  distribution  are  the  temperature  and 
salinity  fields.  Sharp  boundaries  exist  between  masses  of  river, 
coastal,  and  ocean  water  in  the  mixing  area.  Microwave  radio¬ 
meters  have  been  flown  over  this  area  during  various  seasons  of 
the  year  at  altitudes  of  458  meters  and  3048  meters.  These 
radiometers  operate  at  the  central  frequencies  of  9.2,  15.8, 
22.235,  and  34.0  GHz.  Calibrated  apparent  input  temperatures 
have  been  supplied  to  the  author  by  the  Jet  Propulsion  Laborato¬ 
ries  and  the  Manned  Spacecraft  Center  from  aircraft  missions 
flown  in  February  1967  and  February  1968. 

The  Goddard  Space  Flight  Center  conducted  tests  of  several 
remote  sensors  over  portions  of  the  Gulf  of  Mexico  in  June  1967 
using  a  Convair  990.  Included  among  these  remote  sensors  was  a 
microwave  radiometer  which  scans,  electronically,  up  to  50° 
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to  either  side  of  the  flight  track  and  which  operates  at  a  central 
frequency  of  19.35  GHz.  Dr.  William  Nordberg  and  Dr.  Jack  Conaway 
of  the  Goddard  Space  Flight  Center  have  courteously  supplied  the 
author  with  calibrated  temperatures  obtained  by  this  radiometer 
along  with  certain  supporting  photographs  and  infrared  data. 

It  has  been  possible  to  obtain  some  surface  measurements 
made  during  these  aircraft  flights.  The  purpose  of  this  chapter 
is  to  describe  the  microwave  radiometers  flown  on  the  Convair 
240A  and  Convair  990,  to  present  their  measurements  in  a  concise 
form,  and  to  discuss  their  scientific  value  for  oceanography  and 
meteorology . 


Microwave  Radiometers  in  the  NASA  926 

System  description 

Blinn  (1967a)  has  written  a  fairly  complete  description  of 
the  microwave  radiometers  installed  in  a  Convair  240A  which  has 
been  designated  as  the  NASA  926.  There  are  two  separate  micro- 
wave  radiometers  in  this  system  —  the  MR-62  and  the  MR-64.  The 
MR-62  is  practically  identical  tc  the  microwave  radiometers  flown 
on  the  Mariner  II  spacecraft  near  Venus  in  1962.  It  operates 
at  15.8  GHz  and  22.235  GHz  and  uses  a  common  antenna.  The  MR-64 
operates  at  9.2  GHz  and  34  GHz  and  uses  another  antenna.  The 
antennas  of  these  systems  are  shown  in  Figure  18.  The  parabolic 
rt:f lector  of  the  MR-62  is  shown  on  the  left  side  of 


iete: 
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Figure  18,  and  one  can  easily  see  the  series  of  concentric,  flat 
reflectors  which  reflect  visible  light  back  to  space  and  which 
concentrate  microwave  radiation  into  the  horn  in  front  of  the 
reflectors.  One  should  notice  the  size  and  orientation  of  the 
waveguides  which  transfer  the  microwave  signal  from  the  horn 
to  internal  parts  of  the  radiometers.  The  MR-62  is  cross  polarized 
and  the  MR-64  is  polarized  in  the  same  plane  for  both  channels. 

Both  of  these  antennas  may  be  rotated  through  90°  to  change  the 
polarization  planes.  These  antennas  are  mounted  in  the  nose  of 
the  NASA  926  (see  Figure  18  and  19)  and  may  be  tilted  to  zenith 
angles  between  0°  and  45°. 

A  radome,  consisting  of  eleven  layers  of  several  dielectric 
materials,  protects  the  antennas  from  the  airstream.  The  trans¬ 
missivities  (Lr)  of  this  radome  are  given  in  Table  XXVIII  along 
with  some  other  parameters  of  the  radiometers. 


Table  XXVIII.  Significant  parameters  of  the  MR-62  and  MR-64 
microwave  radiometers  on  board  the  NASA  926 
(after  Blinn,  1967a). 
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Figure  19.  Location  of  remote  sensors  in  the  NASA  926. 
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During  a  flight,  the  temperature  of  the  radome  is  recorded  by 
using  two  thermistors  embedded  in  the  radome.  It  was  not  possible, 
however,  to  calibrate  these  thermistors.  Blinn  (1967a)  suggests 
that  one  should  assume  that  the  rad(.  e  temperature  is  equal  to 
the  temperature  of  the  air  recorded  by  standard  aircraft  instru¬ 
ments  . 

Pre-flight  and  post-flight  calibrations  of  both  radiometers 
are  performed  on  the  ground  by  filling  the  horns  of  the  antenna 
with  highly  absorptive  material.  Readings  are  made  of  the  antenna 
temperature  of  this  material  which  is  assumed  to  act  as  a  black- 
body.  This  method  of  calibration  is  liab1*  to  produce  inaccurate 
calibration  points.  The  ’accuracy'  of  these  measurements  is 
suspect;  however,  the  relative  values  can  be  used  with  confidence. 

The  outp..  -  of  the  radiometers  are  recorded  on  magnetic  tape 
by  modulating  the  frequency  of  a  54  KHz  carrier  frequency.  Known 
increments  of  internal  noise  occasionally  are  added  to  the  power 
received  by  the  antenna  system  to  provide  scde  calibration. 

These  microwave  radiometers  have  been  used  to  make  measure¬ 
ments  of  the  tiux  of  upwelling  microwave  radiation  over  an  area 
south  of  the  main  mouths  of  the  Mississippi  River  along  certain 
radial  flight  lines  emanating  from  South  Pass.  Three  extensive 
aircraft  missions  have  been  conducted  over  this  area  in  the  past 
two  years  witv  these  radiometers.  Figure  20  shows  the  basic 
flight  lines  flown  cn  these  missions.  Data  taken  along  Lines  1, 

4,  and  3  will  be  discussed  in  this  chapter. 


Figure  20.  Uasie  flight  lines  flown  over  Mississippi  Delta  over  past  two  years. 
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Mission  41 

Mission  41  was  flown  on  February  24,  1967,  at  an  altitude  of 
458  meters.  Table  XXIX  summarizes  the  significant  flight  para¬ 
meters  and  radiometric  parameters  for  this  mission. 

The  sky  was  clear  during  this  mission.  A  strong  cold  front 
had  swept  over  the  area  in  the  last  24  hour  .  ofore  the  mission 
was  flown.  The  wind  was  about  7  meters  per  oecond  from  the  north 
at  flight  level.  No  surface  measurements  were  available.  The 
rate  of  river  outflow  was  above  normal. 

The  measurements  obtained  by  the  microwave  radiometers  during 
Mission  41  were  calibrated  and  converted  to  digital  form  by  Blinn 
(1967b).  These  data  show  that  the  microwave  radiance  was  fairly 
constant  over  the  entire  area  when  polarization  and  zenith  angles 
were  held  constant.  Occasionally,  the  data  would  indicate  large 
increases  or  decreases  on  all  four  channels.  These  changes  took 
place  over  a  period  of  three  to  five  seconds.  Further  investi¬ 
gations  revealed  that  these  changes  were  approximately  equal  to 
the  noise  levels  occasionally  introduced  into  the  radiometric 
system  to  allow  one  to  calibrate  the  scale  of  the  output.  Thus, 
it  appears  that  these  excursions  in  the  data  were  due  co  improper 
calibration  procedures.  Ignoring  these  excursions,  one  may  use 
the  calibrated  data  to  obtain  some  information  about  the  marine 
environment  on  this  particular  day.  It  will  be  shown  in  the  analy¬ 
sis  of  data  taken  on  Mission  66  that  the  temperature  and  salinity 


Table  XXIX.  Flight  parameters  and  radiometric  parameters  for  Mission  41 
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of  the  water  had  virtually  no  effect  on  the  microwave  radiance  of 
the  sea  for  these  frequencies. 

The  values  of  apparent  input  temperatures  furnished  to  the 
author  by  Blinn  (1967a)  are  averages  over  one-second.  From  these 
data,  ten-second  averages  were  calculated  for  certain  times  along 
the  flight  track. 

An  examination  of  these  values  showed  that  no  steady  trends 
in  antenua  temperature  were  evident  as  one  moved  from  river  water 
to  sea  water.  Average  values  of  apparent  input  temperature  were 
selected  for  each  of  the  flight  lines,  and  these  values  were  cor¬ 
rected  for  radome  effects  by  using  values  of  L  given  in  Table 

K 

XXIX  in  (182).  The  temperature  of  the  radome  was  assumed  to  be 
279°  K  which  was  the  reported  temperature  of  the  air  outside  the 
aircraft  during  the  flights.  The  resulting  values  of  estimated 
brightness  temperatures  are  shown  in  Table  XXX. 

The  slight  disagreement  between  values  of  estimated  brightness 
temperature  over  lines  4  and  1  underscores  the  lack  of  absolute 
accuracy  of  these  measurements.  It  is  apparent  from  the  detailed 
computations  involved  in  arriving  at  values  of  given  in  Table 

XXX  that  the  emissive  temperature  of  the  radome  is  quite  appreci¬ 
able.  The  air  temperature  at  flight  level  was  reported  to  be  279° 
K  throughout  the  daytime  and  nighttime  flights. 

The  main  difference  between  Flights  2  and  3  on  Line  l  is  the 
change  in  incidence  angle  from  10°  to  45°.  Marandino  (1967) 


Table  XXX. 


Mean  values  of  brightness  '■emperatures  for  Mission  41. 
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174 

has  calculated  the  distribution  of  with  incidence  angle  for  a 
flat  sea  surface  using  Fresnel's  Law  of  Reflection.  His  calcula¬ 
tions  show  an  increase  01"  34°  K,  37°  K,  and  36°  K  in  vertically 
polarized  brightness  temperature  when  ij-  changes  from  10°  to  43° 
for  frequencies  of  9.2,  15.8,  22.234,  and  34  GHz,  respectively. 
The  observed  values  of  these  increases  given  in  Table  XXX  are 
40°  K,  2°  K,  19°  K,  and  7°  K,  respectively.  The  discrepancies  in 
the  latter  three  values  are  probably  due  to  sea  roughness  and 
foam  since  some  white  caps  were  visible  in  the  photographs  taken 
on  this  mission.  i  -andino's  calculations  for  horizontal 
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polarization  show  a  decrease  in  the  brightness  temperature  of 
27°  K,  28°  K,  17°  K,  and  24°  K,  respectively,  for  these  same 
frequencies.  The  observed  values  of  these  decreases  given  in 
Table  XXX  are  approximatley  16°  K,  17°  K,  1°  K,  and  13°  K  if  one 
assumed  that  equals  T  at  ii>  equal  to  10°.  Again,  it  appears 
that  roughness  causes  the  horizontally  polarized  brightness 
temperature  to  be  higher  for  a  rough  surface  than  for  a  smooth  sur¬ 
face  for  'j'  near  45° .  Using  the  apparent  wind  speed  at  flight 
level,  one  would  estimate  the  surface  wind  to  be  about  7  meters 
per  second.  Stogrvn  (1967)  predicted  that  the  difference  between 
T^  (10°)  and  (45°)  for  a  frequency  of  35  GHz  and  for  average 
conditions  is  14°  K.  This  value  i.  remarkably  close  to  the  13°  K 
difference  observed.  Due  to  the  extremely  short  fetch  in  this 
case,  one  must  conclude  that  the  local  wind  is  the  controlling 
factor  in  determining  the  effect  of  roughness  on  microwave 
radiation  for  these  high  frequencies. 

The  only  other  information  that  might  be  inferred  from  the 
data  in  Table  XXX  is  about  the  humidity  of  the  atmosphere.  The 
curves  in  Figure  12  are  valid  if  one  is  viewing  the  atmosphere 
from  an  altitude  greater  than  20  kilometers.  Values  of  for  a 
frequency  of  22.234  GHz  were  competed  fcr  an  altitude  of  458  meters 
for  various  values  of  precipitable  water  (W,  cm)  using  Kreiss' 
values  of  sky  temperature  and  values  of  absorption  coefficients 
for  molecular  oxygen  and  water  vapor  listed  in  Table  X.  From 
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these  values  and  the  value  of  in  Table  XXX, it  appears  that  W 
was  approximately  2  centimeters.  This  value  appears  to  be  higher 
than  one  would  expect  on  a  cold  winter  day.  The  increase  in 
apparent  brightness  may  be  due  to  a  small  amount  of  breaking 
waves  on  the  surface. 

Mission  50 

Mission  50  was  flown  on  June  12,  1967  at  458  meters  above- 
sea  level.  The  measurements  made  by  the  microwave  radiometers 
during  this  mission  have  neither  been  calibrated  nor  reduced  to 
digital  form.  The  flight  lines  flown  on  this  mission  are  shown 
in  Figure  21  as  solid  lines.  The  dashed  line,  in  this  figure  are 
the  approximate  flight  lines  flown  by  the  Convair  990  on  June  6, 
1967,  sponsored  by  the  Goddard  Space  Flight  Center.  The  data 
obtained  by  this  flight  (Flight  13)  will  be  discussed  later  in 
this  study. 

Mission  66 

Mission  66  was  flown  on  February  5  end  10,  1968  at  an  altitude 
of  3048  meters  and  458  meters,  respectively.  The  raw  output  from 
the  MR-62  and  MR-64  microwave  radiometers  was  calibrated  and 
converted  to  digital  form  by  the  Manned  Spacecraft  Center  (1968). 
There  were  no  irregularities  in  these  data  as  there  were  in  the 


reduced  data  for  Mission  41. 


The  9.2  GHz  channel,  however,  did 


168 


not  function  properly  except  for  one  brief  segment. 

Flight  1  of  Mission  66  was  flown  on  February  5,  1968  from 

g 

19:23:30  2  *-o  19:29:00  Z  at  an  altitude  of  3048  meters  over  Line 

1  only.  Turbulence  and  fuel  problems  forced  a  cancellation  of 
the  rest  of  this  flight.  The  air  temperature  at  flight  level 
was  274°  K. 

Flight  11  of  Mission  66  was  flown  on  February  10,  1968,  from 
16:08:55  Z  to  16:50:52  Z  at  an  altitude  of  458  meters  and  covered 
Lines  1,  41,  and  4.  The  air  temperature  was  283°  K  at  flight  level 
during  Lines  1  and  41,  and  it  was  284°  K  during  Line  4. 

An  extensive  survey  of  the  test  site  was  conducted  by  R/V 
ALAMIN0S ,  the  departmental  research  vessel  of  Texas  A AM  University, 
from  February  9  to  February  13,  1968,  during  cruise  6S-A-1.  The 
track  chart  for  this  cruise  is  shown  i.  Figure  22.  Figu.es  23 
and  24  show’  the  distributions  of  surface  temperature  and  surface 
salinities  based  on  the  data  obtained  during  cruise  6S-A-1 . 

The  temperature  and  salinity  of  the  surface  layer  are  plotted 
versus  the  distance  from  South  Pass  (d ,  miles)  along  Flight  Line  1 
in  Figure  25. 

The  actual  measurements  made  by  the  microwave  radiometers  were 
furnished  to  the  author  as  one-seccnd  averages  of  apparent 
input  temperature  versus  the  time  (z).  Ten-second  averages  were 
computed  for  every  20  seconds  to  smooth  the  data.  The  apparent 
- g - 

Z  refers  to  Greenwich  Meridian  Time. 
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Figure  22.  Cruise  track  of  R/V  AL AMI NOS  during  Cruise  68-A-l. 


Figure  23.  Analysis  of  surface  temperature  for  68-A-l,  February  1968. 
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DISTANCE  FROM  SOUTH  PASS  ALCUG  LINE  1,  d(n.  «.) 

Figure  25.  The  distribution  of  tempera '-are  and  salinity  a’ong 
Line  1  fror:  analyses  in  Figures  23  and  24. 
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input  temperatures  were  corrected  for  radorae  effects  as  before. 

The  resulting  estimated  brightness  temperatures  are  plotted  versus 
d  in  Figures  26  and  27  along  with  the  theoretical  values  based  on 
the  distributions  of  tenperature  and  salinity  given  in  Figure  25. 
The  measured  values  are  shown  as  circles,  triangles,  or  squares 
on  Cue  following  figures.  The  theoretical  values  are  plotted 
as  solid  lines,  these  lines  are  adjusted  in  absolute  value  to 
allow  one  to  compare  the  remtive  changes  that  occurred  as  the 
flight  line  war  flown. 

The  measured  values  of  brightness  temj  arature  at  a  frequency 
of  9.2  GHz  show  that  this  channel  was  not  operating  properly. 

The  shape  of  the  theoretical  curve  for  15.8  GHz  fits  the  general 
shape  of  the  measured  profile-  There  appears  to  be  a  regular 
fluctuation  in  the  signal  measured  at  458  meters  at  the  frequencies 
of  15-8,  22.234,  and  34  GHz.  On  the  other  hand,  an  overall  drift 
is  apparent  in  the  readings  at  3048  meters. 

In  all  cases,  there  appears  to  be  no  information  in  these 
measurements  about  the  distribution  of  temperature  or  salinity. 

The  fact  that  the  readings  of  the  34  GHz  channel  are  much  higher 
for  low  levels  than  for  middle  levels  gives  one  reason  to  suspect 
the  readings  frt  ■  the  MR-64  on  both  of  these  flights. 

The  mean  value  of  the  readings  of  the  22.235  GHz  channel 
indicate  that  the  precipitable  water  was  approximately  2.5  cm. 
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lg  e  26.  Estimated  brightness  temperatures  measured  over  Line  1  of 

Mission  66  at  altitudes  of  3048  meters  (squares)  and  458 

l*iehtn(cirrlea)  co"esPondin8  theoretical  values  of 

Jf!!!  I**9?* atUreS  (80lid  lines)  versus  Stance  from 
South  Pass  for  frequencies  of  9.2  GHz  and  15.8  GHz. 
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DISTANCE  FROM  SOUTH  PASS  ALONG  LINE  1,  d  (n.  m.) 

Figure  27.  Estimated  brightness  temperatures  measured  over  Lire  i  of 
Mission  66  at  altitudes  of  3048  i..eters  (squares)  and  45P 
meters  (circles  &  triangles)  and  corresponding  theoretical 
values  of  brightness  temperatures  (solid  lines)  versus 
distance  from  South  Pass  for  frequencies  of  22.235  Gdz  and 
34  GHz. 


Implications  for  future  missions 


Using  the  distribution  of  temperature  and  salinity  shown  in 
Figure  25,  the  author  computed  theoretical  values  of  polarized 
brightness  temperature  for  each  nautical  mile  along  Line  1  for 
an  incidence  angle  of  0°  and  for  the  frequencies  of  5. A  GHz  and 
1-0  GHz-  The  effects  of  the  intervening  atmosphere  between  the 
observer  and  the  sea  surface  were  neglected  in  these  calculations. 
The  zenith  sky  temperatures  were  assumed  to  be  10°  K  at  1.0  GHz 
and  5°  K  at  5. A  GHz.  The  results  of  these  calculations  are  shown 
in  Figure  28.  It  has  been  shown  in  Chapters!'  and  VI  that  the 
polarized  brightness  temperature  of  the  sea  is  a  function  of  sea 
temperature  only  at  5. A  GHz.  The  polarized  brightness  temperature 
radiation  at  a  frequency  of  1.0  GHz  is  strongly  dependent  on  salin¬ 
ity  and  slightly  dependent  on  sea  temperature.  Therefore,  it  ap¬ 
pears  that  one  should  be  able  to  survey  the  distribution  of  temper¬ 
ature  and  salinity  along  coastal  areas  using  a  dual-frequency 
radiometer  which  operated  at  5.4  GHz  and  1.0  GHz. 

In  the  summer,  the  distribution  of  temperature  along  coastal 
regions  is  often  inconclusive;  however,  the  salinity  distribution 
does  act  as  an  indicator  of  coastal  mixing  and  other  coastal 
processes.  A  survey  was  made  of  the  area  south  of  the  Mississippi 
Delta  by  the  R/V  ALAMIN0S  in  late  August  1967.  Figures  29  and  30 
show  the  distributions  of  surface  temperature  and  surface  salinity, 
respectively,  derived  from  the  data  obcained  during  this  cruise 
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Figure  28.  Theoretical  oolarized  brightness  temperatures  versus 

distance  from  South  Pass  for  the  distribution  of  tempera¬ 
ture  and  salinity  given  in  Figure  24,  for  frequencies  of 
5.4  GHz  and  1.0  GHz,  and  for  an  incidence  angle  of  0°. 


Figure  30.  Analysis  of  surface  salinity  for  67-A-7,  August  1967 


(67-A-7).  These  figures  show  the  consistency  of  salinity  as  an 
indicator  of  coastal  and  river  phenomena. 

A  new  set  of  microwave  radiometers  is  scheduled  to  be  in¬ 
stalled  in  the  P3A  aircraft  sponsored  by  NASA  as  a  part  of  the 
Earth  Resources  Survey  Program.  One  channel  of  these  radiometers 
operates  at  a  frequency  of  1.4  GHz.  It  is  hoped  that  these 
radiometers  can  be  used  in  April  1969  to  test  the  hypothesis 
presented  in  this  report. 

Microwave  Radiometers  in  the  Convair  990 

Systems  description 

Under  the  sponsorship  of  the  Goddard  Space  Flight  Center, 
an  electronically  scanned  radiometer  operating  at  a  frequency 
of  19.35  GHz  with  a  bandwidth  of  400  MHz  has  been  flown  over  parts 
of  the  Gulf  of  Mexico  during  June  1967  (Catoe,  et  al_. .  1967). 

This  radiometer  uses  an  18  inch  by  18  inch  array  of  dipoles  as  an 
antenna.  The  phase  relationships  between  respective  dipoles  are 
controlled  in  such  a  manner  as  to  form  a  highly  directive  gain 
pattern  that  can  be  step  scanned  sidewavs  up  to  ±50°  from  the 
nadir.  The  ratio  of  the  effective  solid  angle  of  the  main  lobe 
of  this  gain  pattern  to  4tt  steradians  is  about  0.9  or  higher 
depending  on  the  looking  angle.  This  indicates  that  the  phased 
dipole  antenna  is  quite  directive  and  surpresses  sidelobes  ef¬ 
ficiently.  The  sensitivity  of  the  measured  antenna  temperatures 
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is  claimed  to  be  about  2°  K  (Catoe,  et  al. ,  1967).  The  polariza¬ 
tion  is  horizontal.  The  Goddard  Space  Flight  Center  has  calibrated 
and  digitized  the  readings  of  this  radiometer.  The  output  is 
displayed  as  apparent  antenna  temperature  versus  incidence  angle 
and  time.  The  entire  angular  span  of  100°  is  broken  into  39 
separate  angles,  and  the  readings  are  presented  for  each  of  the  scan 
lines  which  are  approximately  2  sec  apart  in  time. 

It  has  been  mentioned  in  Chapter  VI  that  there  appear  to 
be  three  applications  for  this  instrument.  Fir_t,  one  may  be 
able  to  determine  the  toughness  of  the  sea  surface  by  comparing 
the  measured  values  of  horizontally  polarized  radiation  for  ^ 
near  zero  degrees  to  those  at  < l>  near  50°.  Secondly,  the  amount 
of  precipitable  water  in  the  atmosphere  may  be  inferred  from  the 
absolute  value  of  horizontally  polarized  radiation  for  ^  near 
zero  degreees  under  conditions  of  clear  sky.  Thirdly,  it  may  be 
possible  to  determine  the  concentration  of  liquid  water  in  under¬ 
cast  clouds.  These  three  applications  will  be  treated  briefly 
in  the  next  sections. 

Roughness 

This  microwave  radiometer  was  flown  over  the  central  western 
Gulf  of  Mexico  on  June  5,  1967,  at  an  altitude  of  about  10,000 
meters.  Reports  from  surface  ships  indicated  that  the  sea  tempera¬ 
ture  was  approximately  18.0°  C  over  the  area  and  that  the  salinity 
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was  approximately  36.5  o/oo  (Arnold,  personal  communication, 

1968).  Photographs  taken  from  the  Convair  990  during  this  flight 
(Flight  12)  indicated  that  the  sky  was  practically  clear  and  that 
the  surface  was  practically  smooth  and  flat  for  the  portion  of 
this  flight  from  18:13:00  Z  to  18:13:30  Z.  The  readings  of  the 
radiometer  were  averaged  over  this  period  for  each  angular 
position.  These  averages  are  plotted  in  Figure  31  versus  scan 
angle  (ij))  and  are  connected  by  a  dashed  line.  The  theoretical 
distribution  of  horizontally  polarized  brightness  temperature 
with  incidence  angle  given  by  Stogryn  (1967)  for  a  flat  sea  sur¬ 
face  is  shown  in  Figure  31  as  a  solid  line.  The  measured  antenna 
temperatures  and  the  theoretical  brightness  temperatures  agree 
for  i})  near  0°  and  for  \p  near  ±50°.  The  disagreement  between  these 
values  for  ip  near  ±20°  is  probably  due  to  antenna  effects. 

During  the  time  near  18:04  Z  the  sea  surface  was  reported  by  the 

pilot  to  be  a  rough  sea  surface,  no  white  caps - sea  rough 

enough  for  Doppler  to  operate  (Tobin,  1967) .  During  this  period 
of  time,  the  horizontally  polarized  brightness  temperature  at  ±50° 
was  about  4°  K  lower  than  that  at  0°.  From  Table  XXVIII,  it 
appears  that  this  decrease  would  be  caused  by  a  surface  wind  of 
about  10  meters  per  second.  This  value  exceeds  the  critical  wind 
speed  at  which  waves  would  begin  to  break  on  the  surface.  Con¬ 
sidering  the  accuracy  of  these  measurements  and  effects  of  the 
antenna,  one  should  be  satisfied  that  the  emissivity  of  the  water 
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INCIDENCE  Ai:GLE,  (degrees/ 

Figure  31.  The  average  value  of  Che  horizontally  polarized  bright¬ 
ness  temperature  (circles  and  dashed  line)  measured 
during  Flight  12  from  18 : 13: 00 Z  to  18:13:33Z  versus 
incidence  angle  and  corresponding  theoretical  values 
(solid  line)  from  Stogrvn  (1967)  for  a  frequency  of 
19.35  GHe. 
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had  apparently  risen  for  41  near  50°  due  to  the  increase  in 
surface  roughness. 

Precipitable  water 

The  highest  altitude  flown  by  the  Convair  990  was  about 
11,000  meters  above  the  sea  surface.  Figure  32  shows  the  distri¬ 
bution  of  the  estimated  brightness  temperature  for  an  incidence 
angle  of  0°  versus  time  for  a  portion  of  Flight  12.  The  coast  of 
Texas  near  Houston,  Texas,  was  crossed  during  this  time.  One 
should  note,  incidentally,  the  large  decrease  that  occurs  in 
the  estimated  brightness  temperature  when  one  goes  from  land  to 
water.  One  may  even  distinguish  between  the  Intercoastal  Canal 
and  the  narrow  strip  of  land  that  separates  it  and  the  Gulf  of 
Mexico.  It  appears  that  the  average  value  of  the  brightness 
temperature  observed  over  the  land  was  about  148°  K.  The  curves 
in  Figure  12  may  be  used  to  estimate  the  amount  of  precipitable 
water  in  the  atmosphere  for  this  June  day  as  being  about  1.9 
centimeters.  This  value  appears  to  be  reasonable  for  the  season 
and  the  observed  clear  sky. 

Liquid  water 

On  June  6,  1967,  the  Convair  990  was  able  to  fly  over  various 
types  of  clouds  in  the  Gulf  of  Mexico.  One  portion  of  this  flight 
(Flight  13)  passed  directly  over  a  small  group  of  cumulus  clouds. 
The  calibrated  values  of  brightness  temperature  from  the  19.35  GHz 
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radiometer  are  plotted  in  Figure  33  versus  time  for  an  incidence 
angle  of  0°.  The  actual  values  are  shown  as  circles  in  this  figure 
and  the  smoothed  analysis  of  these  data  is  indicated  by  the  solid 
line.  One  should  note  that  a  slight  increase  occurred  in  the 
readings  when  the  visible  edge  of  the  cloud  was  passed.  When 
the  aircraft  was  over  the  core  of  the  cloud,  an  enormous  increase 
was  noted  in  the  brightness  temperature.  It  is  probable  that  this 
increase  occurs  over  portions  of  the  cloud  which  contain  water 
droplets  large  enough  and  dense  enough  to  fall  as  rain.  The  top 
of  the  cloud  was  estimated  to  be  at  an  altitude  of  4600  meters. 

The  cloud  bases  were  estimated  to  be  about  600  meters  above  the 
sea  surface.  The  thickness  of  this  cloud  seems  to  correspond 
most  closely  to  the  case  of  a  'thick  stratus'  cloud  considered 
before  in  Chapter  VI.  If  one  assumed  that  the  curve  for  thick 
stratus  in  Figure  14  is  valid  for  the  cloud  overflown  in  Figure 
32,  one  may  compute  a  horizontal  profile  of  the  mean  concentra¬ 
tion  of  liquid  water  in  this  cloud.  The  computed  values  for  the 
curve  in  Figure  15  were  reduced  by  a  factor  of  0.6  to  allow  for  the 

difference  in  cloud  thickness.  The  resulting  distribution  of  the 

3 

mean  concentration  of  liquid  water  (M,  g/m  )  is  plotted  versus 
time  (Z)  in  Figure  34.  The  increase  in  brightness  temperature 
over  the  brightness  temperature  observed  over  the  clear  ocean  was 
•:sed  in  the  computations.  The  actual  distribution  of  liquid 
water  in  this  particular  cloud  was  unknown.  The  raindrops  probably 
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produced  bubbles  on  the  sea  surface.  In  this  case,  it  is  inadvis¬ 
able  to  use  Kreiss'  curves  for  obtaining  M. 

The  microwave  radiometer  on  board  the  Convair  990  scans 
sideways  to  the  path  of  the  aircraft,  and  it  is  possible  to  obtain 
an  image  or  map  radiometric  values  along  a  strip.  It  is  hoped 
that  this  radiometer  could  be  flown  on  one  of  the  Nimbus  space¬ 
craft.  Using  this  one  radiometer  alone,  it  should  be  possible 
to  survey  rainbands  at  sea,  determine  the  moisture  content  of 
the  atmosphere  over  the  oceans  and  possibly  determine  the  roughness 
of  the  sea  surface.  Land-water  boundaries  and  water-ice 
boundaries  could  be  surveyed  in  remote  areas  of  the  world. 


CHAPTER  IX 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

Microwave  radiometry  has  a  great  potential  as  an  observational 
tool  in  meteorology  and  oceanography.  The  greatest  disadvantage 
is  that  there  are  the  large  number  of  variables  that  significantly 
affect  the  intensity  field  of  microwave  radiation  at  some  location 
above  the  surface  of  the  sea.  Nevertheless,  further  research  in 
microwave  radiometry  is  fully  justified. 

There  appear  to  be  no  serious  gaps  in  the  state  of  knowledge 
concerning  the  various  aspects  of  microwave  radiometry;  however, 
all  aspects  must  be  investigated  more  fully  to  enable  one  to  use 
this  tool  with  confidence.  There  is  a  lack  of  accurate  measure¬ 
ments  of  microwave  radiance  over  the  sea  under  conditions  in  which 
appropriate  environmental  support  measurements  are  available. 

Most  of  the  theoretical  predictions  that  have  been  made  of  micro- 
wave  radiance  have  yet  to  be  confirmed  by  actual  measurements 
over  the  ocean. 

Dielectric  properties  of  water 

The  ability  to  predict  the  emission  and  extinction  of  sea 
water  is  dependent  upon  one's  knowledge  of  its  complex  dielectric 
constant  as  a  function  of  salinity,  temperature,  and  frequency 
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of  observation.  The* fundamental  basis  of  these  investigations, 
at  the  present  tjima,  are  the  observations  by  Lane  and  Saxton 
of  the  electric  behavior  of  an  aqueous  solution  of  sodium  chloride 
made  in  the  1940’s.  Although  there  is  no  reason  to  expect 
that  the  electric  properties  of  such  a  solution  and  that  of 
actual  sea  water  differ  significantly,  it  is  desirable  to  determine 
as  accurately  as  possible  the  electric  properties  of  sea  water 
under  various  conditions  of  salinity  and  temperature  throughout 
the  microwave  region  of  the  electromagnetic  spectrum. 

Radiative  transfer 

It  is  possible  to  specify  the  intensity  of  microwave  radiation 
in  a  horizontally  homogenous,  non-scattering  atmosphere  overlying 
a  homogenous  water  surface  by  solution  of  the  equation  of  radiative 
transfer  if  one  can  specify  the  distributions  of  temperature  and 
absorption  coefficient  with  height  and  depth.  This  approximation 
to  the  real  atmosphere  and  hydrosphere  is  valid  as  long  as  there 
are  no  clouds  of  large  water  droplets  in  the  atmosphere  and  as  long 
as  the  sea  surface  is  relatively  undisturbed.  Otherwise,  scatter¬ 
ing  cannot  be  neglected.  If  one  is  to  consider  scattering,  the 
problem  of  solving  the  equation  of  radiative  transfer  becomes 
enormously  complicated  since  the  entire  field  of  intensity  must 
be  specified  at  each  point  in  the  atmosphere.-  In  many  cases, 
the  real  atmosphere  is  not  horizontally  uniform,  for  example,  in 
the  vicinity  of  thunderstorms  and  fronts.  These  situations,  too. 
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must  be  considered  in  the  future. 

Thermodynamic  equilibrium 

The  source  function  for  thermal  emission  by  atmospheric  and 
hydrospheric  constituents  is  assumed  to  be  Planckian  for  all  cases 
considered.  This  form  is  valid  only  if  a  condition  of  thermo¬ 
dynamic  equilibrium  exists  throughout  the  entire  environment.  It 
has  been  noted  that  the  amount  of  thermal  energy  transferred  as 
microwave  radiation  is  much  smaller  than  that  transferred  as  in¬ 
frared  radiation.  Suppose  that  a  radiational  imbalance  existed  in 
some  location  in  the  atmosphere  as  far  as  a  particular  band  of 
microwave  frequencies  is  concerned.  One  would  expect  the  result¬ 
ant  rate  of  heating  or  cooling  to  be  insignificantly  small. 

Thus,  the  idea  that  this  imbalance  would  be  lessened  greatly  through 
changes  in  internal  energy  appears  to  be  invalid  in  this  case. 

A  final  proof  for  any  theory,  of  course,  lies  in  its  ability  to 
predict  results  that  agree  with  measurements  under  controlled 
conditions  or  well-known  conditions  in  the  environment. 

Water  vapor  and  molecular  oxygen 

The  theory  of  absorption  and  emission  by  water  vapor  and 
molecular  oxygen  in  the  atmosphere  was  formulated  by  Van  Vleck 
in  the  1940's.  The  distributions  of  these  atmospheric  constitu¬ 
ents  are  easily  specified  in  the  case  of  oxygen,  but  the  distri¬ 
bution  of  water  vapor  is  usually  unknown.  It  appears  that  the 


r 


193 


natural  variability  of  water  vapor  in  the  atmosphere  has  little 
effect  on  the  amount  of  microwave  energy  absorbed  or  emitted  by 
the  atmosphere  when  the  frequency  is  less  than  8  GHz.  Most  studies 
of  the  microv/ave  emission  of  water  vapor  are  made  by  assuming  a 
simple  expcnential  decay  of  its  density  with  height.  Real  distri¬ 
butions  should  be  considered  to  determine  fully  if  this  variable 
can  be  surveyed  remotely  through  measurements  of  appropriate 
microwave  radiances. 

It  appears  that  the  several  engineering  approaches  to  the 
problem  of  emission  of  water  vapor  in  the  atmosphere  yield  only 
approximations.  If  one  is  able  to  use  relative  measurements 
of  microwave  radiance  as  opposed  to  absolute  measurements,  it  is 
not  necessary  to  be  able  to  predict  accurately  the  contribution 
of  the  overall  radiation  that  comes  from  the  atmosphere. 

Hydrometers 

The  effect  of  extinction  or  emission  by  ice  crystals  in  the 
atmosphere  can  be  neglected  in  most  cases.  The  extinction  of 
liquid  water  is  linearly  proportional  to  the  concentration  of 
liquid  water  in  a  cloud  of  wate*-  droplets  if  the  maximum  size  or 
the  droplets  is  small,  for  example,  as  in  the  case  of  droplets 
in  fog  and  non-raining  clouds.  Extinction  in  raining  clouds 
becomes  more  complicated  and  scattering  must  be  considered.  Since 

c 

scattering  and  absorption  of  clouds  of  large  droplets  must  be 
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treated  by  the  precise  theory  developed  by  Mie,  the  distribution 
of  drop  sizes  in  the  cloud  becomes  a  factor.  The  investigations 
in  this  area  are  elementary  at  present  and  much  work  is  needed. 

One  needs  to  ascertain  the  effect  of  bubbles  and  foam  that  are 
formed  on  the  sea  surface  where  rains  fall  into  the  sea. 

Air-sea  interface 

Since  the  ocean  returns  approximately  two-thirds  of  the  micro- 
wave  radiation  incident  upon  its  surface  to  the  atmosphere  as 
scattered  or  reflected  radiation,  it  is  necessary  to  determine  the 
effect  of  various  configurations  of  sea  surfaces.  Emipirical 
methods  must  be  used  to  specify  the  scattering  of  the  rough  sea 
as  a  function  of  incidence  angle,  polarization,  sea  temperature, 
sea  salinity,  electromagnetic  frequency,  and  roughness.  Based  on 
meager  evidence  .j  date,  it  appears  that  the  sea  surface  behaves 
as  a  specular  reflector  for  viewing  angles  near  normal  incidence 
(Stogryn,  1967).  If  foam  exists  on  the  sea  surface  as  a  result 
of  large  surface  shearing  stresses,  one  must  expect  roughness 
to  be  a  factor  at  all  angles. 

Sea  temperature 

Due  to  the  relaxation  of  water  molecules,  the  emissivity  of 
sea  water  decreases  with  increasing  temperature  for  frequencies 
above  about  8  GHz.  Current  microwave  measurements  lie  in  this 
range  of  frequencies,  therefore,  it  is  not  possible  to  determine 
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the  •  'anperature  of  the  sea  surface  by  available  radiometers. 

It  appears  that  the  emissivity  of  sea  water  remains  fairly 
constant  or  increases  slightly  with  temperature  for  frequencies 
near  5.4  GHz. 

Salinity 

The  effect  of  changes  in  the  salinity  of  sea  water  on  the 
brightness  temperature  of  sea  water  is  negligible  for  frequencies 
above  4  GHz.  Investigations  of  this  effect  for  frequencies  above 
1  GHz  revealed  that  the  effect  of  salinity  is  dominant  for  fre¬ 
quencies  between  1  GHz  ana  3  GHz  and  that  the  effect  of  salinity 
is  largest  at  1.0  GHz  and  for  viewing  angles  near  normal  incidence. 
It  appears  that  microwave  radiometers  operating  at  1.0  GHz  and 
5.4  GHz  could  be  used  to  survey  remotely  temperature  and  salinity 
along  coastal  and  river  regions  of  the  world.  Imaging  radiometers 
operating  at  these  frequencies  might  be  able  to  map  the  location 
of  major  water  mass  boundaries. 

Foam,  bubbles,  ice,  and  land 

Microwave  measurements  and  theory  indicate  that  one  should  be 
able  to  determine  the  amount  of  foam  or  ice  covering  the  ocean 
surface  at  some  location  and  to  survey  land-sea  and  ice-sea 
boundaries  remotely  by  using  scanning  microwave  radiometers.  Ac¬ 
curate  measurements  are  needed  to  determine  the  effect  of  foam 
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The  current  generation  of  microwave  radiometers  should  be 
accurate  enough  to  allow  scientists  to  make  microwave  measurements 
that  can  be  used,  practically,  for  a  large  number  of  applications. 
Sensitivities  of  0.5°  K  are  needed  in  most  applications  to  ocean- 
graphy.  Presently  available  microwave  radiometers  are  sufficiently 
accurate  for  meteorological  applications. 

Applications 

In  summary,  the  following  oceanographic  measurements  appear 
feasible  from  airborne  or  spaceborne  microwave  radiometers: 

(1)  determination  of  surface  roughness,  (2)  determination  of 
surface  temperature,  (3)  determination  of  surface  salinity,  (A) 
detection  of  icebergs,  (5)  survey  of  ice  fields,  (6)  determination 
of  foam  cover,  and  (7)  survey  of  precipitation. 

The  following  meteorological  measurements  seem  feasible: 

(1)  measurement  of  moisture  content  of  atmosphere  from  spacecraft, 

(2)  survey  of  rainbands,  (3)  determination  of  concentration  of 
liquid  water  in  non-raining  clouds,  and  (A)  sounding  of  temperature 

c 

(not  covered  in  detail  in  this  study). 

Recommendations 

The  following  suggestions  are  offered  for  future  research  and 


action: 
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1.  Numerous  measurements  of  the  microwave  radiance  over  the 
ocean  are  needed  under  known  environmental  conditions  so  that 
existing  microwave  theory  may  be  tested  throughly  and  so  that  cer¬ 
tain  empirical  determinations  can  he  made,  for  example,  the  nature 
of  the  bistatic  scattering  cross-section  for  a  real  ocean  sur¬ 
face.  Atmospheric  soundings  of  temperature,  humidity,  and  pres¬ 
sure  are  needed  to  support  any  measurement  program  for  meteorolog¬ 
ical  applications.  An  appropriate  set  of  microwave  radiometers 
should  be  mounted  on  a  research  vessel  so  that  accurate  surface 
measurements  can  be  made  of  the  brightness  temperature  of  the  sky 
and  the  sea.  To  avoid  the  effects  of  being  in  the  near  field, 
one  would  mount  antennas  on  a  high  platform. 

2.  Flights  should  he  made  of  the  microwave  radiometers 
scheduled  to  be  placed  on  the  P3A  aircraft  sponsored  by  NASA  over 
the  outflow  region  of  the  Mississippi  Delta  to  test  the  feasibility 
of  measuring  salinity  with  the  1.4  GHz  channel. 

3.  Extensive  studies  of  the  electric  properties  of  sea  water 
as  a  function  of  its  temperature,  its  salinity,  and  its  dissolved 
gases  should  be  made  for  frequencies  throughout  the  microwave 
regions.  Controlled  studies  of  the  emission  of  sea  foam  and 
bubbles  should  be  included  in  this  study. 

4.  The  effects  of  liquid  water  in  the  atmosphere  should  be 
studied  rigorously  using  the  Mie  equations  to  specify  scattering 
and  using  realistic  models  of  atmospheric  clouds.  This  study 
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would  constitute  a  major  study  into  practically  all  area  of  concern 
to  microwave  physics.  Computer  programs  would  have  to  be  developed. 
Simultaneous  surveys  of  the  atmosphere  by  radiometer,  radar,  and 
in  situ  instruments  would  need  to  be  made  under  various  meteorolog¬ 
ical  conditions. 
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